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Abstract 

In  rarefied  gas  flow  problems  there  are  two  effects  which  influence 
the  flow;  (1)  collisions  of  gas  molecules  with  other  gas  molecules,  and  (2) 
collisions  of  g£.s  molecules  with  solid  surfaces  (the  gas-surface  interaction). 
This  study  deals  with  free  molecular  flow  in  which  the  effect  of  gas-gas 
collision  can  be  neglected  and  the  gas-surfacc  interaction  has  the  dominate 
influence  on  the  flow.  The  results  of  this  study  have  application  to  satel¬ 
lites  since  free  molecular  flow  conditions  occur  at  orbital  altitudes  above 
120  km. 

Knowledge  of  the  gas-surface  interaction  is  required  in  order  to 
determine  the  aerodynamic  properties  of  satellites.  At  satellite  velocities 
(7-8  km/sec)  the  interaction  of  neutral  atmospheric  gas  molecules  with  the 
satellite  surface  occurs  at  energies  in  the  1  to  10  ev  range.  It  is  just 
this  energy  range  which  has  not  been  satisfactorily  duplicated  in  the  labora¬ 
tory;  therefore,  at  present,  laboratory  gas-surface  interaction  data  can  not 
be  applied  directly  to  the  determination  of  the  aerodynamic  properties  of 
satellites.  It  is  proposed  in  this  study  that  satellite  experiments  be 
performed  to  obtain  the  needed  information  from  measurements  of  the  aero¬ 
dynamic  properties  of  satellites.  In  order  to  interpret  the  satellite  data, 


a  generalized  gas-surface  interaction  model  was  developed  and  used  in  the 
analysis  of  this  study. 

Gas-surface  interaction  models  such  as  those  of  Maxwell,  Schamberg, 
and  Nocilla,  contain  two  or  more  parameters  which  may  be  adjusted  to  cover  a 
certain  range  of  possible  gas- surface  interactions.  Although  such  specific 
models  may  be  used  to  develop  the  aerodynamic  equations  of  satellites,  the 
validity  of  these  models  in  this  application  has  not  been  determined.  The 
results  of  this  study  show  that  the  proposed  generalized  model  is  necessary 
in  the  interpretation  of  measured  satellite  aerodynamic  properties. 

In  the  past,  the  interpretation  of  measured  satellite  aerodynamic 
properties  to  obtain  information  on  the  gas-surface  interaction  and  orbital 
gas  density  has  not  been  successful  for  two  reasons;  (1)  the  uncertainty  in 
the  validity  of  gas-surface  interaction  models,  and  (2)  insufficient  data  to 
allow  a  determination  of  the  orbital  gas  density  and  at  least  two  gas-surface 
interaction  parameters.  The  results  of  this  study  illustrate  strongly  the 
feasibility  of  performing  a  satellite  experiment  in  which  accurate  data  could 
be  obtained  on  the  gas  density  and  gas-surface  interaction  parameters  by 
measuring  the  drag,  3pin  rate  slowdown  and  spin  axis  precession  rate  of  a 
spinning  convex  satellite. 

The  results  of  the  study  on  the  aerodynamic  properties  of  spinning 
convex  bodies  have  exhibited  a  number  of  interesting  effects  associated  with 
the  spin  of  the  body  and  the  gas-surface  interaction.  For  example,  the  drag 
and  lift  of  a  spinning  body  was  found  to  be  greater  than  that  cl  a  non-spin¬ 
ning  body.  It  was  also  found  that  there  exists  a  spin  induced  lateral  lift 
force  which  is  analogous  to  the  Magnus  effect  but  is  opposite  in  direction. 


In  addition,  spin  induced  aerodynamic  torques,  perpendicular  to  the  spin 
axis,  are  significant  on  bodies  at  angles  of  attack  to  the  flow. 

The  gas-surface  interaction  was  found  to  have  a  strong  influence 
in  determining  the.  aerodynamic  properties  of  both  spinning  and  non-spinning 
bodies.  Both  analytical  and  numerical  results  were  obtained  for  the  aero¬ 
dynamic  properties  of  four  basic  body  shapes  (disk,  cylinder,  cone,  and 
sphere)  to  study  the  effects  of  spin,  angle  of  attack,  and  the  gas-surface 


interaction. 


ACKNOWLEDGMENTS 


in 


The  author  thanks  Professor  Shee-Mang  Yen  who,  as  the  principal 
advisor,  has  generously  contributed  his  time  and  wisdom  in  the  prepara¬ 
tion  of  this  dissertation  and  Professor  Howard  W.  Knoebel  who,  with 
Professor  Yen,  was  directly  responsible  for  the  able  direction  of  this 
study. 

Since  June  of  1968,  this  study  has  been  performed  under  a  con¬ 
tract  with  the  NASA  George  Marshall  Space  Flight  Center.  The  author  is 
deeply  indebted  to  Mr.  James  0.  Ballance  of  NASA  for  his  active  support 
of  this  work. 

The  author  thanks  those  who  have  provided  encouragement  and 
guidance  through  discussions;  Dr.  Kenneth  Moe,  Dr.  James  L.  Myers,  Jr., 
Dr.  Franklin  M.  Propst,  and  Dr.  Wallace  W.  Youngblood. 

The  author  thanks  Richard  Flood  who  did  the  computer  program¬ 
ming  in  the  numerical  studies,  Arnold  Z.  Cohen  who  gathered  information 
on  past  satellites  and  the  atmpspheric  environment,  and  Nicki  Stillings 
who  carefully  typed  this  manuscript. 

The  author  takes  pleasure  in  thanking  his  wife,  Laurel,  for  her 
encouragement  and  her  many  sacrifices. 


4  ",  s  „  *'■?  >V  *^7“^}  V?' <^3i" ’  * 


iv 

TABLE  OF  CONTENTS 

Page 

1.  INTRODUCTION .  1 

2.  A  GENERALIZED  GAS -SURFACE  INTERACTION  MODEl .  7 

2.1.  A  Generalized  Model  for  the  Interaction .  7 

2.1.1.  Reflected  Velocity .  9 

2.1.2.  Angle  of  Reflection . 11 

2.1.3.  Reflected  Number  Flux..... .  14 

2.1.4.  Parameters  of  the  Generalized  Model . 15 

2.2.  Subclasses  of  the  Generalized  Model .  15 

2.2.1.  Reduction  to  Maxwell  Model .  16 

2.2.2.  Reduction  to  Schamberg  Model .  17 

2.2.3.  Reduction  to  Nocilla  Model . 23 

2.3.  Incorporation  of  Experimental  Results . .  29 

2.3.1.  Current  Status  of  Molecular  Beam  Experiments...'...  29 

2.3.2.  Incorporation  of  Intensity  Distribution .  30 

2.3.3.  Incorporation  of  Force  Measurements . . .  31 

2.4.  Significance  of  the  Proposed  Model .  34 

3 .  BASIC  AERODYNAMICS  EQUATIONS .  37 

3.1.  Basic  Equations  of  Force  and  Torque . . .  37 

3.2.  Aerodynamics  Equations  for  a  Class  of  Spinning  Body  Shapes  39 

3.2.1.  Coordinate  Systems . "45 

3.2.2.  Spinning  Disk  at  Angles  of  Attack.. .  48 

3.2.3.  Spinning  Cylinder  at  Angles  of  Attack .  52 

3.2.4.  Spinning  Cone  at  Angles  of  Attack . 55 

3.2.5.  Spinning  Sphere  at  Angles  of  Attack .  60 

3.3.  Applications . 65 

4.  AERODYNAMIC  PROPERTIES  OF  SELECTED  SHAPES .  66 

4.1.  Special  Case  of  the  Maxwell  Model .  66 

4.1.1.  Spinning  Disk  Properties  in  Terms  of  the 

Maxwell  Model  Parameters .  70 


V 


Page 

4.1.2.  Spinning  Cylinder  Properties  in  Terms  of  Maxwell 

Mode1.  Parameters . 71 

4.1.3.  Spinning  Cone  Properties  in  Terms  of  Maxwell 

Model  Parameters . 73 

4.1.4.  Spinning  Sphere  Properties  in  Terms  of  the 

Maxwell  Model  Parameters .  76 

4.2.  Solutions  in  Terms  of  the  Generalized  Model  Parameters...  77 

4.2.1.  Analytical  Results  (Zero  Spin) .  77 

4. 2. 1.1.  Drag  and  Lift  of  Flat  Plate  at  Angles 

of  Attack . 78 

4. 2. 1.2.  Drag  of  Cylinder  with  Axis  Perpen- 

t  dicular  to  Flow .  79 

4. 2. 1.3.  Drag  of  Cone  with  Axis  Parallel  to  Flow.  80 

4. 2. 1.4.  Sphere  Drag .  81 

4.2.2.  Numerical  Methods .  81 

4. 2. 2.1.  Single  Variable  of  Integration .  82 

4. 2. 2. 2.  Double  Integration .  82 

4. 2. 2. 3.  Graphical  Display  of  Results .  83 

4.3.  Discussion  of  Results .  84 

4.4.  Conclusions .  95 

.  .  < 

5.  AERODYNAMIC  PROPERTIES  OF  SATELLITES .  98 

5.1.  Coordinate  Transformation .  98 

5.2.  Instantaneous  Aerodynamic  Properties  of  Satellite .  101 

5.2.1.  Components  of  Force .  101 

5.2.2.  Components  of  Torque  Acting  on  Satellites .  104 

5.3.  Average  Aerodynamic  Properties  of  Satellites..., .  106 

5.3.1.  Spinning  Disk . . .  107 

5. 3. 1.1.  Spinning  Disk  with  One  Side  Exposed 

to  the  Flow . r .  108 

5. 3. 1.2.  Spinning  Disk  with  Both  Sides  Exposed 

to  the  Flow. . . .  108 

5.3.2.  Spinning  Spherical  Satellite .  109 

5.3.3.  Spinning  Cone-Disk  Composite  Satellite .  110 

5.4.  Average  Aerodynamic  Properties  of  Tumbling  Bodies .  Ill 


vi 

Page 

5.4.1.  Average  Drag  Coefficient  of  a  Random-Tumbling.- 

Disk .  113 

5.5.  Discussion  of  Results . . .  114 

6.  FEASIBILITY  OF  A  SATELLITE  EXPERIMENT  TO  DETERMINE  THE  GAS- 
SURFACE  INTERACTION  PARAMETERS  AND  THE  ATMOSPHERIC  DENSITY. ...  116 

6.1.  Gas-Surface  Interaction  Experiments  Utilizing  Satellite. 

Aerodynamic  Properties .  >117 

6.1.1.  Method  1:  Utilization  of  Data  From  Satellites 

of  Different  Shapes .  119 

6.1.2.  Method  2:  Utilization  of  Data  From  a  Single 

Satellite .  121 

6.2.  Accuracy  of  a  Satellite  Experiment..., .  123 

6.3.  Magnitudes  of  Measurable  Quantities  in  a  Satellite 

Experiment . 128 

6.4.  Assessment  of  Results  Obtained  From  a  Satellite 

Experiment . . .  136 

6.5.  Significance  of  Satellite  Experiments  to  Determine 

Gas"Surface  Interaction  Parameters  and  the  Atmospheric 
Density .  141 

6.5.1.  Consideration  of  Surface  Conditions .  141 

6.5.2.  Consideration  of  the  Composition  of  the 

Atmosphere .  143 

6.5.3.  Variation  of  Atmospheric  Density  Models .  145 

7.  SUMMARY  AND  CONCLUSIONS .  146 


REFERENCES 


150 


! 

I 

Q 


0 

0 

0 

0 

0 

0 


1 


1 .  INTRODUCTION 

In  the  free  molecular  flow  regime,  intermolecular  collisions  may 
be  neglected  and  the  gas-surface  interaction  then  becomes  the  dominating 
influence  in  this  flow  regime.  For  earth  satellites,  free  molecular  con¬ 
ditions  exist  at  all  altitudes  above  100  mi  (161  km.).^  It  is  therefore 
essential  that  the  effect  of  the  gas-surface  interaction  be  considered  in 
the  determination  of  the  aerodynamic  properties  of  satellites. 

At  satellite  velocities,  the  gas  molecules  of  the  atmosphere  im¬ 
pinge  on  the  satellite  surface  at  velocities  in  the  order  of  7  to  8  km/sec- 
Taking  into  account  the  molecular  weight  of  the  molecules  composing  the 
atmosphere,  the  interaction  energy  associated  with  satellite  velocities  is  in 

the  range  of  1  to  10  ev.  Laboratory  experiments  using  molecular  beam  tech- 
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niques  have  not  been  successful  in  duplicating  these  interaction  energies. 
Therefore,  the  character  of  the  gas-surface  interaction  at  satellite  velo¬ 
cities  and  thus  the  effect  of  the  interaction  on  satellite  aerodynamic  pro¬ 
perties  has  not  been  determined.  It  has  been  possible  to  construct  models  of 
the  gas-surface  interaction  from  physical  principles  and  experimental  results 
obtained  for  interaction  energies  less  than  1  ev.  Such  models  contain  two 
or  more  parameters  which  may  be  adjusted  to  include  a  certain  range  of  pos¬ 
sible  interaction.  It  has  not  been  determined,  however,  how  well  these 
models  approximate  the  actual  gas-surface  interaction  that  occurs  in  the 
satellite  environment. 

The  fact  that  aerodynamic  properties  of  bodies  in  a  free  molecular 
flow  depend  on  the  gas-surface  interaction  suggests  that  measurements  of 
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satellite  aerodynamic  properties  could  yield  information  on  the  gas-surface 
interaction.  However,  measurements  of  satellite  aerodynamic  drag  have  been 
inconclusive  in  determining  information  on  the  gas-surface  interaction  for 
two  reasons;  1)  uncertainties  in  the  satellite  environment;  primarily,  un¬ 
certainty  in  the  atmospheric  density  and  2)  uncertainties  in  the  aerodynamic 
properties  of  satellites  and,  thus,  in  the  interpretation  of  the  measure¬ 
ments.  Drag  measurements  really  only  determine  the  product  of  density  and 
drag  coefficients  since  neither  is  known  separately. 

The  problems  with  interpreting  drag  measurements  suggest  that  an 

additional  aerodynamic  property  of  a  satellite  should  be  measured,  such  as 

the  slowdown  rate  of  a  spinning  satellite,  which,  when  combined  with  the 

drag  measurement  of  that  satellite,  would  provide  a  means  of  separating  the 

effects  of  density  and  drag  coefficient  or  gas-surface  interaction.  Analyses 

of  this  type  have  been  performed  on  drag  and  spin  rate  decay  data  for  paddle- 

wheel  shaped  satellites,  from  which  estimates  of  the  density  and  a  gas-sur- 
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face  interaction  parameter  were  obtained.  *  These  results  are,  however, 
subject  to  uncertainties  which  are  much  the  same  as  those  associated  with 
drag  measurements.  First,  since  the  measurements  of  drag  and  slowdown  rate 
were  a  function  of  at  least  three  unknowns  (orbital  gas  density  and  two  or 
more  gas-surface  interaction  parameters),  a  value  for  at  least  one  of  the 
unknowns  had  to  be  assumed  in  order  to  obtain  estimates  of  the  other  two. 
Secondly,  since  the  validity  of  any  particular  model  of  the  gas-surface  inter¬ 
action  has  not  been  established,  the  interpretation  of  measurements  may  be 
made,  as  in  References  2  and  3,  using  a  number  of  different  models. 


Since  each  model  used  produces  a  different  estimate  of  density,  the  experi¬ 
ment  is  able  to  determine  only  a  possible  range  on  the  density. 

The  first  point  made  above  suggests  that  additional  aerodynamic 
properties  should  be  measured  in  order  to  remove  the  necessity  to  make 
assumptions  which  can  introduce  errors  in  the  interpretation  of  the  measure¬ 
ments.  There  are  six  aerodynamic  properties  to  be  considered  corresponding 
to  three  components  of  force  and  three  components  of  torque  acting  on  a 
satellite;  however,  the  properties  must  net  only  be  measurable  but  must  also 
be  independent  functions  of  the  quantities  to  be  determined.  The  free  mole¬ 
cular  aerodynamic  properties  (drag,  lift,  and  torque)  of  non-spinning  bodies 
are  known  to  depend  on  the  angle  of  attack  of  the  body  and  the  gas-surface 
interaction  (see  for  example  References  5  and  6),  For  non-spinning  satel¬ 
lites,  however,  the  cts  of  lift  and  torque  properties  cannot  be  easily 
assessed  because  the  orientation  of  the  satellite  with  respect  to  the  flow 
is  usually  unknown  and  probably  random.  Spinning  satellites,  on  the  other 
hand,  maintain  relatively  fixed  orientations  in  space.  If  the  orientation 
of  the  spi|i  axis  of  the  satellite  is  known,  1l  is  possible  to  determine  the 
angle  of  attack  on  the.  satellite  with  respect  to  the  flow  at  any  position 
in  the  orbit.  This  suggests,  then,  that  the  aerodynamic  properties  of  spin¬ 
ning  satellites  may  provide  the  measurables  needed  for  determining  the  gas 
density  and  gas-surface  interaction  parameters. 

Past  studies  of  the  free  molecular  aerodynamic  properties  of 
spinning  bodies  do  not,  however,  provide  a  sufficient  basis  for  proposing 
a  satellite  experiment  such  as  suggested  above.  The  analysis  of  aerodynamic 
torque  on  spinning  satellites  is  usually  made  on  an  approximate  basis 
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considering  only  the  moments  of  drag  forces  about  the  center  of  mass  of  the 
satellite  (see  for  example  References  7,  8)  which,  in  general,  do  not  give  a 
complete  understanding  of  the  influence  of  the  gas-surface  interaction. 

More  exact  analyses  of  the  aerodynamic  torques  have  been  made  in  References 
9,  10,  and  11  for  the  case  of  a  spinning  spherical  satellite  in  which  it  was 
found  that  the  aerodynamic  torque  properties  are  strongly  dependent  upon  a 
single  parameter  of  a  specific  gas-surface  interaction  model. 

The  objective  of  this  study  is  to  analyze  more  fully  the  influence 
of  the  gas-surface  interaction  on  the  aerodynamic  properties  of  spinning 
bodies  and  to  propose  satellite  experiments  to  accurately  determine  the  gas 
density  and  the  gas- surface  interaction. 

In  order  to  remove  uncertainties  introduced  by  a  variety  of  pos¬ 
sible  gas-surface  interaction  models  that  can  be  used  in  such  a  study  and 
satellite  experiments  (second  point  made  aoove)  ,  a  generalized  gas-surface 
interaction  model  is  developed  which  is  designed  to  cover  a  wider  range  of 
possible  gas-surface  interactions  than  models  currently  being  used.  The 
generalized  model  contains  currently  accepted  models  as  subclasses  and  has 
the  additional  advantage  of  being  able  to  incorporate  laboratory  results 
and  models  which  may  be  suggested  in  the  future.  The  description  of  this 
generalized  model  is  given  in  chapter  2. 

In  chapter  3,  the  generalized  model  is  used  to  develop  the  equa¬ 
tions  expressing  the  aerodynamic  properties  of  spinning  and  non-spinning 
bodies  in  a  free  molecular  flow.  The  results  obtained  may  be  interpreted 
in  terms  of  any  of  the  gas-surface  interaction  models  contained  as  sub¬ 
classes  in  the  general  model.  The  aerodynamic  equations  are  developed  in 
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a  general  manner  which  makes  them  applicable  to  bodies  of  various  shapes. 

In  chapter  4,  results  are  obtained  for  a  disk  (or  flat  plate), 
cylinder,  cone,  and  sphere  for  arbitrary  angles  of  attack  and  for  both  the 
spinning  and  non-spinning  cases.  These  results  reveal  the  strong  influence 
of  the  gas -surface  interaction  on  the  aerodynamic  properties  of  spinning 
bodies . 

In  chapter  5,  the  aerodynamic  properties  of  spinning  satellites 
is  studied  to  determine  the  importance  of  the  gas-surface  interaction  on  the 
average  aerodynamic  properties  of  satellites.  These  results  suggest  pos¬ 
sibilities  for  performing  satellite  experiments.  The  random  tumbling  pro¬ 
blem  is  also  studied  in  this  chapter. 

In  chapter  6,  satellite  experiments  are  proposed  and  the  feasi¬ 
bility  of  performing  these  experiments  is  investigated  by  assessing  the 
possible  accuracy  and  the  magnitude  of  measurable  quantities  needed  to 
determine  the  unknowns  of  atmospheric  density  and  gas-surface  interaction 
parameters . 

The  feasibility  of  the  proposed  satellite  experiment  is  enhanced 
by  results  obtained  in  a  study  performed  by  the  Coordinated  Science  Lab¬ 
oratory  pertaining  to  the  measuremer.c  of  satellite  precession  rates  which 
could  be  caused  by  a  general  relativity  effect.  This  study  determined  that 
extremely  accurate  measurements  of  even  small  precession  rates  are  possible 

by  using  a  completely  passive  optical  readout  techniqie  utilizing  observa- 
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tions  of  sunlights  reflected  by  the  satellite  surface.  ’  On  this  basis 
then,  it  is  proposed  that  for  certain  satellite  shapes  there  are  at  least 
three  measurable  aerodynamic  properties  (drag,  slowdown  torque,  and 


precession  torque)  which  can  be  utilized  to  determine  more  precisely  the 
atmospheric  density  and  the  character  of  the  gas-surface  interaction  at 
satellite  velocities. 
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2.  A  GENERALIZED  GAS-SURFACE  INTERACTION  MODEL 

Insufficient  information  is  available  at  present  to  warrant  choos¬ 
ing  a  specific  gas-surface  interaction  model  to  represent  the  reflection  of 
molecules  which  impinge  a  surface  at  satellite  velocities.  For  this  reason, 
a  generalized  gas-surface  interaction  model  is  developed  such  that  it  con¬ 
tains  various  possible  gas-surface  interaction  models  or  subclasses  includ- 
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ing  the  models  of  Maxwell,  Nocilla  ,  and  Schamberg 

The  generalized  model  will  be  used  in  the  development  of  equations 
to  express  the  aerodynamic  properties  of  spinning  bodies  in  subsequent 
chapters.  The  resulting  equations  have  the  advantage  that  they  can  then  be 
interpreted  in  terms  of  any  of  the  gas-surface  interaction  models  contained 
as  subclasses  of  the  general  model. 

2.1.  A  Generalized  Model  for  the  Interaction 

Consider  a  stream  of  mono-energetic,  -mi-directional  neutral  gas 
molecules  impinging  upon  a  solid  surface  at  an  angle  of  0  with  respect  to 
the  plane  of  the  surface  (see  Figure  2.1).  Also  consider  that  the  molecules 
are  all  reflected  in  a  beam  which  is  axial  symmetric  about  an  axis  which 
makes  an  angle  9.  with  respect  to  surface,  in  the  plane  formed  by  the  imping- 

J 

ing  molecules  and  the  surface  normal  (see  Figure  2.1).  The  subscript,  j, 
on  8^  may  take  on  values  of  l,2,3---to  represent  cases  in  which  the  reflec¬ 
tion  can  be  modeled  as  being  composed  of  two  or  more  beams  which  are  axial 
symmetric  about  axes  which  make  angles  with  respect  to  the  surface  of  9^, 

0.,  0,  .  respectively.  The  purpose  for  adding  the  versatility  of 

4  j 
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Figure  2.1.  Notation  used  in  the  generalized  model. 
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of  using  more  than  one  reflected  beam  will  be  illustrated  later  in  this 
chapter.  For  the  present,  however,  consider  the  beam  in  the  0^  direction. 

In  order  to  describe  the  gas-surface  interaction,  the  following 
three  quantities  must  be  determined  about  the  reflected  beam. 

1.  The  velocity. 

2.  The  angle  of  reflection. 

3.  The  number  flux. 

The  relations  of  these  three  quantities  with  the  incident  beam  properties 
are  described  in  the  following  three  sections. 

2.1.1.  Reflected  Velocity 

In  general,  the  velocity  of  individual  molecules  reflected  from 

the  surface  will  be  distributed  in  some  arbitrary  manner.  The  distribution 

of  velocities  of  a  large  number  of  reflected  molecules  could  be,  for  example, 

Maxwellian,  or  constant  (no  distribution),  or  any  one  of  any  numerous 

possible  distributions.  For  purposes  of  calculating  the  force  on  the  surface 

in  free  molecular  flow,  however,  the  specific  distribution  of  velocities  is 

not  important  since  once  the  molecules  leave  the  surface  they  do  not  again 

hit  the  surface,  and  they  do  not  collide  with  the  impinging  molecules. 

Only  the  average  velocity  of  the  reflected  beam  is  needed  in 

determining  the  momentum  of  reflected  molecules  and  then  the  force  on  the 

surface.  Therefore,  a  vector  velocity  U.  is  defined  to  represent  the 

average  velocity  of  the  beam  of  molecules  reflected  in  the  direction  0^. 

The  velocity  U.  is  also  in  the  direction  of  0.  since  the  beam  is  assumed  to 
J  j 

be  symmetrical  in  velocity  ~..stribution  about  the  axis  at  angle  0.. 
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with  T^,  the  temperature  of  the  reflected  molecules.  The  operation  in 
going  from  temperature  to  velocity  in  Equation  2.4  depends  on  the  distribu¬ 
tion  function  of  reflected  molecules  and  also  on  the  definition  of  temper¬ 
ature.  For  the  purposes  of  this  study,  Equation  2.2  will  be  used  without 
reference  to  temperature. 

2.1.2.  Angle  of  Reflection 

The  reflection  of  a  beam  of  molecules  from  a  solid  surface  was 
first  considered  by  Maxwell  to  be  analogous  to  the  reflection  of  light 
from  a  surface.  He  postulated  that  molecules  could  be  reflected  elastically 
or  specularly  much  as  light  from  a  perfect  mirror,  or  the  molecules  could 
reflect  diffusively  as  light  does  from  a  rough  surface  (see  Figure  2.2). 

Even  though  this  treatment  of  the  angle  of  reflection  may  be  elementary, 
the  Maxwell  model  has  found  wide  applications. 

Jn  1959,  R.  Schamberg^  proposed  a  gas-surface  interaction  model 
which  allowed  for  reflections  at  angles  between  the  limits  of  specular  and 
diffuse.  Schamberg  postulated  that  the  angle  of  reflection  should  be 
related  to  the  angle  of  incidence  of  the  molecular  beam..-  As  an  example  of 
such  a  relationship,  Schamberg  introduced  a  parameter,  v,  defined  by 

cos  0r  =  (cos  0i)V;  v  >  1  2.5 

where  9^  is  the  angle  of  the  reflected  beam  of  molecules  and  6^  is  the  angle 
of  incidence.  In  the  limits  of  u  =  l  and  v  -*  the  Schamberg  model  reduces 
to  the  cases  of  specular  and  diffuse  angles  of  reflection  respectively. 


RS-AZl 


Figure  2.2.  NotP'cion  used  in  the  Maxwell  model 
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Although  recent  experimental  results  using  molecular  beam  tech¬ 
niques  indicate  that  Equation  2.5  is  not  in  general  correct,  experimental 
results  do  indicate  that  the  angle  of  reflection  is  a  function  of  the 
angle  of  incidence. (the  experimental  results  will  be  discussed  later  in 
this  chapter).  In  order  to  include  an  angle-oi:-ref  lection  law  in  the 
generalized  model,  consider  the  functional  relationship  between  the  angle  of 


reflection,  9  ,  and  the  angle  of  incidence,  0,  to  be  in  the  form  of  a 
general  polynomial  of  degree  N.  That  is,  let 


N 

8,  =  a.  +  2  (b.)  9n  2.6 

j  J  i  J  n 

j  j  n=i  j 


where  a^  and  (b^)^  n  =  1,2, ...N  are  constant  coefficients. 

As  appropriate  experimental  results  become  available,  the  constants 

a.  and  (b  )  can  be  found  by  fitting  a  polynomial  to  the  experimental  data. 

J  ^ 

Since  appropriate  experimental  data  is  not  now  available,  assume,  as  a  first 
approximation,  that  the  functional  relationship  for  the  angle  of  reflection 
is  linear  in  0.  That  is,  let 


9.  =  a .  +  b  .0 
J  3  J 


2.7 


which  contains  the  two  unknown  constants  a.  and  b..  Unless  there  is  a 

J  J 

systematic  irregularity  in  the  surface,  the  reflection  of  a  molecule  beam 
which  is  incident  normal  to  the  surface  (8  =  should  also  be  normal  to 
the  surface  (9.=  j) .  Using  this  reasoning,  one  of  the  unknown  constants  in 
Equation  2,8  can  be  eliminated  to  obtain  a  functional  relationship  dependent 


upon  only  one  unknown  constant,  P^  ,  where 

TT 


9.  =  f  p.  +  ci-Pj)  e 


2.8 
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This  form  of  the  angle-of-ref lection  law  will  be  used  in  later  chapters  to 
illustrate  the  effect  of  angle  of  reflection  on  the  aerodynamic  properties 
of  convex  bodies  in  free  molecular  flow. 


2.1.3.  Reflected  Number  Flux 

If  a  solid  surface  is  neither  a  source  or  sink  for  molecules,  the 
number  flux  of  reflected  molecules  must  equal  the  number  flux  of  incident 
molecules.  The  number  flux  of  incident  molecules,  N,  is  defined  as 

N  =  £  u-n  2.9 

m 

— ■> 

where  p  is  the  density  of  the  incident  gas,  n  is  the  unit  normal  to  the 
surface,  and  U  is  the  incident  velocity  with  respect  to  the  surface.  If 
all  the  molecules  were  reflected  in  a  single  beam,  ®refiecte(j  ~  N  =mU-n. 
For  the  generalized  model,  a  parameter  a ^ ,  is  introduced  which  relates  the 
number  flux  reflected  in  the  0^.  beam,  ,  with  the  incident  number  flux. 


N.  =  <7.  N  =  a.  m  U-n  2.10 

J  J  J  m 

For  cases  when  all  the  incident  molecules  are  reflected  in  the  0^  direction, 
cjj  =  1.  For  the  more  general  case  when  the  reflection  is  composed  of  J 
symmetric  beams  having  direction  0^,  . 0j, 


A*’1 

where  J  is  the  number  of  reflected  beams. 


vvr  ~  *~izr-  -1 


n 
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2.1.4.  Parameters  of  the  Generalized  Model 

Three  separate  parameters  have  been  introduced  which  define  the 
gas_surface  interaction  in  a  generalized  manner.  The  average  velocity  of 
the  reflected  molecules  is  related  to  the  velocity  of  the  incident  molecules 
by  the  parameter  where  the  magnitude  of  IL  is  given  by 


U  =  U/l-or. 
J  J 


2.12 


The  velocity  is  a  vector  having  direction  defined  by  the  angle  of  reflec¬ 
tion  .  As  a  first  approximation,  0^.  is  related  to  the  angle  of  incidence 
0  by  the  parameter  P^  given  by 

0.  =  ?  P.  +  (1-P.)  0  2.13 

j  2  j  j 

The  number  flux  of  molecules  reflected  in  the  beam  which  is  symmetrical 
about  the  0^.  axis  is  related  to  the  incident  number  flux  by  the  parameter 
a.  piven  by 


where 


mN.  =  a .  N  =  a.  p  U*n 
J  J  J  K 


S  a.  =  1 


j=]  j 

and  where  J  represents  the  number  of  symmetric  beams, 


2»,2.  Subclasses  of  the  Generalized  Model 

For  each  reflectad  beam,  the  parameters  <j \ >  a ^ ,  and  P  must  bo 
specified  to  determine  the  force  on  the  surface.  By  proper  choice  of  these 
parameters,  the  generalized  model  can  be  reduced  to  more  specific  gas-surface 


1 


} 
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interaction  models.  This  procedure  is  illustrated  in  the  following  three 
sections  for  the  models  of  Maxwell,  Schamberg,  and  Nocilla, 


4  * 


*  m 

’i 

.1 


f 


I 

«6 


2,2,1.  Reduction  to  Maxwell  Model 

In  the  Maxwell  reflection  model,  the  reflection  is  divided  into 

two  components,  specular  and  diffuse.  Define 

=  fraction  of  incident  molecules  2.14 

which  are  reflected  diffusely 

(1-or^)  =  fraction  of  incident  molecules  2.15 

u  which  are  reflected  specularly 

The  velocity  of  molecules  reflected  specularly  is  defined  as  being  equal  to 
the  incidetit  velocity  and  angle  of  reflection  is  equal  to  the  angle  of 
incidence  (elastic  collision  with  the  surface)  (see  Figure  2.2).  For  the 
-diffusely  reflected  component  of  the  reflection  consider  the  velocity  of 
reflection  to  be  related  to  the  incident  velocity  by  the  thermal  accommoda¬ 
tion  coefficient  where 


T.  -  T 
_  t  r 

aT  T  -  T 
i  w 


2.16 


where  T.,  T  ,  and  T  are  defined  as  in  Equation  2.2. 
i  r  w  n 

The  Maxwell  model  is  obtained  from  the  generalized  model,  as 


follows . 


First,  consider  the  specular  component  of  reflection  and  let  this 
be  beam  j  =  1.  Then,  let 

flj  =  0  or  =  U 

=  (l-ffd-)  2.17 

6,  =  0  or  ?.  =  0 

1  i 


La t;  the  diffusely  reflected  component  of  reflection  be  beam  j  =  2. 


Then ,  le t 


a2  = 


-  Ttj 

=  or  U,  =  0;  f- «  1 


a2 


e 


2 


or 


2.18 


The  force  on  a  surface,  due  to’  the  impingement  and  subsequent 
reflections  of  the  Maxwell  type,  is  then 

F  -  r  U-h  (U  -  -  a2U2) 

=  p  U4n  (U  -  (l-ad)  \  -  adU2)  2.19 

where  =  U  having  direction  0  and  U2  =  U  having  direction 


2.2.2.  Reduction  to  Schamberg  Model 

The  Schamberg  model  already  has  much  in  common  with  the  proposed 
generalized  model  in  that  the  reflection  is  considered  as  being  in  a  beam 
which  is  axially  symmetric  about  an  angle  which  is  not  necessarily  in  the 
specular  or  diffuse  directio- .  The  principle  differences  between  the  two 
models  is  the  manner  in  which  the  reflected  velocity  and  the  angle  of  reflec¬ 
tion  are  defined. 

In  the  Schamberg  model,  a  specific  form  for  the  distribution  of 
velocity  in  the  reflected  beam  is  given.  The  velocity  of  reflected  mole¬ 
cules  are  assumed  to  be  equal  in  magnitude  but  distributed  in  direction  and 
number  according  to  a  cosine  law  within  a  beam  width  of  angle  0q  (see  Figure 
2.3)  given  by 
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r,r(0)  =  K  cos  Q-  •  D 


where  is  the  number  of  re-emitted  molecules  per  unit  time  whose  direction 
of  re-emission  (relative  to  the  axis  of  the  beam)  lies  between  0  and  (0+d0) . 
The  constant  K  is  related  to  the  number  flux  of  reflected  molecules  and  is 
dependent  upon  whether  the  reflected  beam  is  two-dimensional  (wedge  shape) 
or  three-dimensional  (conical  shape). 

To  reduce  the  generalized  model  to  the  Schamberg  model,  the 
average  velocity  of  the  reflected  beam  must  be  found.  For  a  three-dimen¬ 
sional  conical  beam  the  average  reflected  velocity  can  be  determined  from 
expressions  derived  by  Schamberg, 

BJ  *  WVr  2M 

where  $3(0^  is  defined  in  Schamberg' s  1959  paper,  and  Vr  is  the  magnitude 

of  the  constant  velocity  of  individual  molecules  in  the  reflected  beam. 

The  quantity  ^Wg)  has  a  maximum  value  of  one  for  0q  =  0  and  a  minimum 

value  of  2/3  for  0q  =  The  velocity  V  is  related  to  the  incident  velocity 

(or  U  in  the  notation  being  used  for  the  generalized  model)  by  a  thermal 

accommodation  coefficient , a,  which  has  been  defined  in  Equation  2.2,  i.e. 

V  =  JUa  V.  2.22 

r  ^  1 

Substituting  Equation  2.22  into  the  expression  for  Uj  in  Equation  2.21,  the 

following  expression  relating  U.  to  U  is  obtained. 

U.  »  §- (0  )  U  2.23 

j  Jo 

Therefore,  in  order  to  reduce  the  generalized  model  to  the  Schamberg  model, 
the  parameter  a ^  must  be  defined  as 
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a/1-ct  «  Jl-a  $3  (0q)  2.24 

The  two  parameters,  a  and  0  ,  of  the  Schamberg  model  are  then  reduced  to  one 
parameter  a  . 

The  angle-of-ref lection  law  proposed  by  Schamberg  is  given  by 

cos  9^  ==  (cos  8)U  2.25 


and  was  discussed  briefly  in  section  2.1.2.  A  plot  of  Equation  2.25  for 
various  values  of  u  is  given  in  Figure  2.4.  A  rough  approximation  to 
variation  in  9^  as  a  function  of  8  for  the  Schamberg  model  can  be  made  by 
discontinuous  linear  relationships.  For  example, 


8 .  =  a .+b .6,  for  9  <  9 
J  J  J 

-  c^d.  8,  for  8  >  8 
j  J  " 


2.26 


For  the  Schamberg  model  Equation  2.25,  the  6^  vs  8  curves  all  pass  through 
the  (0,0)  and  C^>^)  points.  Using  this  information,  a^  =  0,  c^  =  ^  (1-d.). 

A 

Also,  using  the  fact  that  at  9  =  8  the  two  lines  intersect,  Equations  2.26 
reduce  to 


9.  = 
J 


fa 

L  9 


* 


for  9  <  8 


=  ^  (i-d.)  +  d.e 


k 


for  9  £  0 


2.27 


A  one  parameter  family  of  discontinuous  linear  curves  can  be  developed  from 

“it 

Equations  2.27  for  a  choice  of  a  relationship  between  9  and  8.  This  is 

k 

illustrated  in  Figure  2.5  where  8  was  chosen  to  occur  along  the  line  from 
(0,  ~)  to  (^,0).  The  approximation  illustrated  in  Figure  2.5  retains  the 
essential  characteristics  of  the  Schamberg  cosine  variation. 


Figure  2.4.  Plot  of  Schataberg's  angle-of-ref lection  law. 
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2.2,3.  Reduction  to  Nocilla  Model 
14 

Nocilla  postulated  that  the  distribution  of  re-emitted  molecular 

velocities  be  a  drifting  Maxwellian  having  a  drift  velocity  of  in  the 

direction  0  .  The  distribution  function  of  reflected  velocities  is  written 
r 


f  (C)  = 
r 


- Z  exp 

(2TrRirr 


r-ifiL-i 

L  2RT  J 


2.28 


where  C  is  the  molecular  velocity.  Since  number  flux  must  be  conserved  at 

the  surface,  n  can  be  related  to  the  incident  number  flux.  There  are,  then, 
r 

three  parameters  remaining  to  describe  the  reflection,  U  ,  0  ,  and  T  . 

/ 

Nocilla  has  shown  that  the  model  can  be  made  to  closely  match  the 

16 

distribution  obtained  experimentally  by  Hurlbut  ,  for  proper  choice  of  the 

quantities  U^,  0^  and  T  .  However,  Nocilla  doesn't  propose  an  angle-of- 

ref lection  law  or  a  relation  between  the  incoming  and  reflected  velocities. 

Therefore,  in  order  to  develop  an  interaction  model  using  the  distribution 

function  proposed  by  Nocilla,  these  relationships  must  be  provided. 

The  Nocilla  distribution  function  for  reflected  molecules  has  been 

applied  to  the  calculation  of  forces  on  a  solid  surface  in  free  molecular 

6 

flow  by  Hilrlbut.  and  Sherman.  Their  results  can  be  used  to  show  that  the 
Nocilla  model  is  a  subclass  of  the  proposed  generalized  model. 

The  force  on  a  surface  due  to  a  reflection  of  the  Nocilla  type  is 
divided  by  Hurlbut  and  Sherman  into  components  in  the  direction  cf  the 
incoming  beam,  D^,  and  perpendicular,  L^.  (See  Figure  2.6). 
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n^niU^  , 

Dr  =  \  V~2  (ar  +  ~0  T1‘  (l+erf  a  )}  sin  0 

r  2tt*S  2  L  r  r  2  r 

-  sr  cos  er  [x(or)}  cos  bJ 

i*  iuu  —  .  * 

y»  J  1  3L 

Lr  =  2  UCrx(a7  +  2  "  (l+erf  or)}  cos  9 

+  sr  cos  er  {x(cr)}  sin  bJ 

S  =  U  / (2RT 


where 


2.29 


n  raU 
t  r 

2-rr^S 


2.30 


2.31 


a  =  S  sin  0 
r  r  2r 

O’  i 

X(ar)  =  e  r  +  TT^ar(l+erf  ar) 

If  the  term  x(ar)  is  factored  out  of  Equations  ?.29  and  2.30,  and 

if  the  expression  for  a  in  Equations  2.31  is  used,  D  and  L  become. 

r  r  r 


nrraU.2  r  (1+erf  a  ) 

Dr  =  7172  X(ar}  K  CCS  (6  +  «r>  +  ^ 


2tt  S 


n  raU 


2  X(ar) 


sin  0 


2.32 


“r‘“v'r  T  "  (^"er£  o  ) 

Lt  "  ~T~I  x(ar)  LSr  Sin  (Q  +  6r)  + - 77777  ~  ~  cos  0  > 


2tt  S_ 


2X(or) 


From  these  equations,  it  is  -apparent-  that  the  vector  force  on  the  surface 
due  to  a  reflection  of  the  Nocilla  type  can  be  divided  into  components 
normal  to  the  surface  and  in  the  direction  of  the  velocity  U  .  From 
Equations  2.32 


n  raU 

?  _  r  r 

reflection  •  h  2 
£TT  S 


r  U  h  tt2  (1+erf  a  ) n- 

K  T  + - 


2.33 


ft". 


The  outgoing  number  flux  Nr  is  given  by 

,V  T  -  j. 


■r(s 


X(ar) 


Substituting  Equation  2.33  into  2.34  and  using  the  conservation  of  mass 


flux  condition 


mN  -  mN  =  pU*n 
r 


Equation  2.33  reduces  to 


F  ..  .  .  ~  9  U’ft  I  -U  +  6(a  )  n  ! 

reflection  L  r  r  J 


where 


(2LT  y  r  L 

Pa+er'‘°r]  2-37 

In  terms  of  the  proposed  generalized  model,  a  reflection  of  the 

Nocilla  type  can  be  considered  to  be  composed  of  two  beams,  one  normal  to 

surface  and  one  in  the  direction  of  Ur .  The  parameter  Oj  is  not  needed 

since  conservation  of  mass  flux  is  automatically  satisfied  by  the  two  beams, 

Let  beam  i  =  1  be  in  the  U  direction.  Since  one  is  free  to  choose 
J  r 

relationships  for  the  and  6r»  let 

Ur  =  U.  end  @r  -  *  Pr  +  <1_Pr)  9  2,33 

Consider  beam  j  =  2  to  be  normal  to  the  surface.  Then 

U,  =  G(c;r)  =  §(ar)cr  2.39 


where 


cr  =  (2RTr)% 
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then  from  Equation  2.37,  define 


‘«V  =  f 


k  1+erf  a 


X(ar) 


2.40 


Assuming  the  the  velocity  is  related  to  the  incident  velocity  u  by  the 
thermal  accommodation  coefficient  a, 


c  =  Jl-a  U 
r 


2.41 


we  can  define  the  parameter  as 


=  Jl-d  §(or)  2.42 

The  reduction  of  the  generalized  model  to  a  model  using  the 
Nocilla  reflected  distribution  function  is  now  complete.  The  model  obtained 
still  contains  three  parameters;  ,  Pr ,  and  c^.  This  form  of  model  serves 
to  illustrate  some  interesting  characteristics  of  the  Nocilla  distribution 
function.  The  forces  caused  by  a  reflection  of  the  Nocilla  type  are  seen  to 
have  the  character  of  a  Maxwell  reflection  model  except  that,  instead  of  a 
specular  component,  a  component  is  employed.  Also,  similar  to  the  Max¬ 
well  model,  the  magnitude  of  the  component  of  reflection  is  seen  to 
determine,  except  for  one  parameter,  the  magnitude  of  the  normal  component 
of  reflection.  The  normal  component  of  reflection  is  analogous  to  the 
diffusion  component  in  the  Maxwell  model.  Figure  2,7  shows  a  plot  of  the 
coefficient  of  the  normal  component  of  the  Nocilla  distribution  function, 
§(<Tr),  as  a  function  ar>  The  term  a  is  equal  to  the  product  uf  and 
sin  9r  divided  by  .  From  Figure  2.7,  as  Ur  becomes  large  compared  to  c^, 
the  magnitude  of  the  normal  component  becomes  small.  At  0  =  0,  the  normal 


Figure  2.7.  Plot  of  the  coefficient  of  the  normal  component 
of  force  for  the  Nocilla  model  as  a  function  of 
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component  has  a  maximum  value  of  .89  c  .  Due  to  its  dependence  on  sin  6r> 

the  quantity  cr  can  not  in  general  be  assumed  to  always  be  large.  However, 

under  the  special  conditions  of  Ur  being  much  greater  than  and  greater 

than  zero  (this  condition  on  P  insures  that  sin  0  and  therefore  a  will 

r  r  r 

not  be  zero),  the  quantity  §(0^)  may  be  assumed  to  be  small  in  comparison 
to  the  component  of  reflection.  Under  these  conditions,  the  Nocilla 
distribution  function  resembles  the  type  of  reflection  described  by  the 
generalized  model  with  the  average  velocity  replaced  by  the  velocity 
of  the  Nocilla  distribution. 

2,3.  Incorporation  of  Experimental  Results 

2.3.1.  Current  Status  of  Molecular  Beam  Experiments 

Molecular  beam  experiments  have  not  as  yet  been  able  to  obtain 
results  which  could  be  directly  applicable  to  the  calculation  of  forces  and 
torques  on  a  satellite.  The  major  limitation  in  molecular  beam  experiments 
is  the  inability  to  produce  a  neutral  molecular  beam  of  sufficient  intensity 
and  at  a  velocity  which  corresponds  to  the  velocity  of  impingement  of 
atmospheric  molecules  on  a  satellite  in  near  earth  orbit.  In  terms  of  the 
energy  of  interaction,  the  range  of  1  to  10  ev  corresponds  to  that  which 
occurs  in  a  near  earth  orbit.  Interaction  energies  below  1  ev  have  been 
obtained  by  aerodynamic. methods  such  as  expansion  of  a  high  pressure  gas 
through  a  nozzle.  On  the  high  energy  side  of  the  1  to  10  ev  range  (10  ev 
and  above)  neutral  molecular  beams  have  been  obtained  by  the  method  of  ion 
acceleration  with  subsequent  neutralization  by  charge  exchange  techniques 
(see  for  example  Reference  2). 
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Although  the  energy  of  interaction  is  of  primary  importance  in 
determining  the  character  of  the  gas-surface  interaction,  other  factors 
such  as  the  satellite  surface  conditions  and  composition  of  gases  incident 
on  the  surface  are  also  believed  to  influence  the  interac* i cn.  The  major 
limitation  in  duplicating  these  factors  in  molecular  bea a  experiments  has 
been  that  satellite  surface  condition  and  atmospheric  composition  remain 
uncertain. 

As  can  be  concluded  from  the  above  discussion,  available  ex¬ 
perimental  results  are  probably  not  applicable  for  determining  accurately 
the  values  of  the  gas-surface  interaction  parameters  for  the  calculation 
of  aerodynamics  forces  on  satellites.  However,  molecular  beam  experiments 
outside  the  1  to  10  ev  range  could  for  example  suggest  the  form  of  the  9^ 
vs  0  relationship,  or,  indicate  trends  in  the  character  of  the  interaction 
which  could  be  extrapolated  to  the  1  to  10  ev  range. 

2.3.2.  Incorporation  of  Intensity  Distribution 

The  majority  of  experimental  work  is  aimed  toward  obtaining 
information  on  the  distribution  function  of  reflected  molecules  because 
all  other  flow  properties  can  be  found  from  the  distribution  function. 

For  application  to  calculating  forces  on  convex  satellite  shapes,  however, 
less  detailed  information  can  be  used  since  the  actual  form  of  the  distri¬ 
bution  function  is  unimportant.  The  generalized  gas-surface  interaction 
model  developed  in  section  2.1  suggests  the  type  of  experimental  data 
which  would  be  most  useful  for  satellite  application.  For  example, 
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distribution  of  reflected  intensity  provides  information  on  the  number  of 
beam  components  (single  lobed,  J  =  1;  bi-lobed,  J  =  2;  etc).  Also,  inten¬ 
sity  distributions  obtained  for  various  angles  of  attack  to  the  incident 
beam  reveals  characteristics  of  the  0^  vs  0  relation.  As  an.  example,  the 
experimental  results  of  0^  vs  0  for  high  velocity  argon  beams  on  heated 
platimum  obtained  by  Moran,  Wachman,  and  Trilling^  are  plotted  in  Figure 
2.8.  These  results  show  that  as  the  beam  velocity  increases,  the  para¬ 
meter  Pj  (using  the  first  approximation  for  the  0  relationship)  approaches 
zero.  These  results  show  a  slight  departure  from  a  linear  relation  in  @r 
vs  0  at  low  incidence  angles.  The  departure  is  opposite  to  that  postulated 
by  Schamberg’s  (see  Figure  2.4)  cosine  powered  relationship. 

Overspecular  (0^  <  0)  and  backscatter  (0^  >  -j)  results  can  not 
be  reproduced  in  the  Schamberg  or  the  Maxwell  model.  Such  results  are, 
however,  easily  incorporated  into  the  generalized  model.  For  example, 
backscattering  results  could  be  approximated  by  a  linear  relation  by 
allowing  P^  to  take  on  values  between  1  and  2.  Overspecular  results  and 
over-backscatter  (0^  >  tt  -  0)  could  be  approximated  by  2nd  or  3rd  degree 
polynomials  or  by  discontinuous  linear  relations.  These  regions  are 
indicated  in  Figure  2.9  with  examples  of  the  possible  angular  relationships. 


2.3.3.  Incorporation  of  Force  Measurements 

Force  measurements  made  on  flat  surfaces  at  angles  of  attack  to  a 
molecular  beam  can  yield  considerable  information  on  the  gas-surface  inter¬ 
action  and  the  results  are  particularly  suited  for  analysis  in  terms  of 
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iralized  model. 
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Figure  2.9.  The  full  range  of  possible  angle-of-ref lection 
laws  in  the  generalized  model . 
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the  generalized  model.  Since  the  force  on  the  surface  due  to  the  impinge¬ 
ment  can  be  obtained  from  measured  beam  velocity  and  intensity,  the  force 
due  to  the  reflection  of  molecules  can  be  directly  obtained.  Therefore, 
the  magnitude  and  direction  of  the  reflection  force  can  be  determined  to 


yield  information  on  the  0  vs  9  relationship  and  also  the  value  of  cr^/l-Q'.. 

In  order  to  properly  analyze  force  measurements  in  this  manner,  however, 

it  is  necessary  to  know  the  number  of  beam  components  which  make  up  the 

reflection.  Experiments  set  up  to  measure  forces  are  not  usually  set  up 

to  measure  the  distribution  of  intensity  which  could  reveal  the  number  of 

beam  components.  In  the  absence  of  such  information,  special  measuring 

techniques  can  be  employed  to  obtain  the  desired  information  from  the  force 

18 

measurements.  For  example,  the  technique  employed  by  Boring  and  Humphris 
to  obtain  reflected  force  in  the  directions  normal  and  tangential  to  the 
surface  can  be  extended  to  obtain  measured  values  of  the  reflected  force 
at  angles  between  those  two  limits.  The  results  would  reveal  the  lobal 
maxima  of  the  beam  components  and  thereby  yield  information  on  the  9. 


relationship  and  the  quantities  aj*/l-a\ 


2.4.  Significance  of  the  Proposed  Model 

The  generalized  gas-surface  interaction  model  developed  in  this 
chapter  was  shown  to  be  reducible  to  three  currently  accepted  gas-surface 
interaction  models.  In  addition,  the  generalized  model  is  found  to  be 
applicable  to  approximating  a  wide  range  of  possible  gas-surface  interaction 
by  including  the  possibilities  of  over-specular  and  backscatter  reflection. 
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The  generalized  model  is  found  to  also  be  useful  in  the  interpretation  of 
laboratory  experiment  results,  as  illustrated  in  Figure  2.8,  in  order  to 
indicate  trends  and  similarities  in  laboratory  data. 

The  generalized  model  described  in  this  chapter  has  incorporated 

1)  the  effect  of  the  angle  of  incidence  on  the  angle  of  reflection  and 

2)  the  effect  of  the  incident  velocity  on  the  reflected  velocity.  These 

are  undoubtedly  the  major  factors  iL  encing  the  gas-surface  interaction. 
Additional  factors,  however,  can  readily  be  incorporated  into  the  gen¬ 
eralized  model.  For  example,  the  velocity  of  reflection  could  be  considered 

to  be  a  function  of  the  angle  of  incidence  as  well  as  a  function  of  the 

incident  velocity.  This  effect  is,  in  fact,  incorporated  in  the  results 
obtained  in  the  reduction  of  the  generalized  model  to  the  Nccilla  model. 

The  generalized  model  could  also  be  extended  to  include  the  possibility 

that  ct,  and  a.  are  both  functions  of  U  and  0.  Although  refinements  such  as 
J  J 

these  will  not  be  included  in  the  application  of  the  generalized  model  in 
the  work  which  follows,  they  may  be  incorporated  whenever  warrented. 

The  generalized  model  employs  three  parameters  Gy  ay  and  P^ 
for  each  beam  of  the  reflection.  The  two  parameters  Oj  and  determine 
the  magnitude  of  the  force  due  to  the  reflection  and  the  parameter  P^ 
determines  the  direction  of  that  force.  When  using  the  generalized  model, 


as  will  be  seen  in  the  following  chapters,  the  quantity  a^Ji-a^  can  be  con¬ 
sidered  as  a  single  parameter  in  place  of  both  and  a ^  .  That  is,  speci¬ 


fication  of  the  two  quantities  a  and  P^  are  sufficient  for  determin¬ 

ing  the  force  and  torque  acting  on  a  convex  body  in  a  free  molecular  flow. 
In  the  analysis  of  subsequent  chapters,  the  generalized  model 
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will  be  employed  in  the  development  of  equations  which  express  the  aero¬ 
dynamic  properties  of  bodies  in  a  free  molecular  flow. 
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3 .  BASIC  AERODYNAMICS  EQUATIONS 

The  generalized  gas-surface  interaction  model  is  applied  in  this 
chapter  to  the  development  of  the  equations  for  expressing  force  and  torque 
acting  on  spinning  bodies  in  free  molecular  flow.  The  resulting  equations 
have  wide  application  since  they  can  be  interpreted  in  terms  of  many  specific 
models . 

The  equations  for  drag,  lift,  and  torque  are  obtained  in  differ¬ 
ential  form  for  four  body  shapes  in  particular  (disk,  cylinder,  cone,  and 
sphere).  The  results  will  be  used  in  the  next  chapter  to  study  the  effects 
of  spin  and  the  gas~suriace  interaction  on  the  aerodynamic  properties  of 
bodies  of  these  shapes. 

3.1.  Basic  Equations  of  Force  and  Torque 

Consider  a  surface  element,  dA,  on  a  spinning  body.  If  the  posi- 
tive  normal  of  that  element  of  surface  is  n,  the  mass  flux  impinging  on  the 
surface  is 

-  p  U  •  n  dA 

where  p  is  the  density  of  the  free  stream  gas  and  U  is  the  vector  velocity 
of  the  incident  flow  relative  to  the  surface.  The  velocity  U  is  assumed  to 
be  much  higher  than  the  thermal  motion  of  the  free  stream  gas  so  that  ran¬ 
dom  fluctuations  in  the  velocity  can  be.neglected  (hypothermal  assumption). 

The  mass  flux  is  a  positive  number  since  the  quantity  (U-n)  must  always  be 
negative  j n  order  for  molecules  to  hit  the  surface.  The  vector  force  on  the 
surface  due  to  the  impingement  of  molecules  is  given  by 
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(dF).  . ,  =  U(-p  U-n  dA)  3.2 

'incident  v  K 

■ 

where  the  notation  dF  is  used  to  denote  the  force  on  the  surface  area  dA. 

Using  the  notation  introduced  in  the  preceding  chapter  for  the 
generalized  model,  the  vector  force  due  to  the  reflection  of  molecules  is 
given  by 


(dF)  „  =  -a.  U.  (-pU-n  dA) 

reflection  j  j 


3.3 


where  the  minus  sign  is  required  because  the  vector  velocity  must  be  in 
the  direction  away  from  the  surface.  For  cases  in  which  the  reflection  is 
composed  of  more  than  one  beam  component,  the  force  due  to  the  reflection 
is  given  by 

^reflection-"  +  +-"+  “  3'4 

where  J  is  the  number  of  beams  making  up  the  reflection.  Since  the  re¬ 
peated  subscripts  can  be  used  to  imply  the  summation,  Equation  3.3  will  be 
used  in  the  development  of  equations  which  follow. 

The  total  vector  force  on  the  element  of  surface  is  the  sum  of 
force  due  to  impingement  and  force  due  to  reflection  given  by 


dF  =  -(U  ~  a.U.)  p  U-n  dA  3.5 

J  J 

Consider  the  element  of  surface  dA  to  be  at  a  point  defined  by 
the  radius  vector,  R,  from  the  center  of  a  coordinate  system  attached  to 
the  body.  The  vector  torque,  dT,  about  the  center  of  the  coordinate  system 
caused  by  the  impingement  and  reflection  of  molecules  on  the  clement  of  sur¬ 
face  dA  is  then 


v,. 


/ 
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dT  =  -Rx  (U  -  GjUj)  p  U-n  dA 


3.6 


Equations  3.5  and  3.6  form  the  basic  equations  necessary  for 
determining  the  force  and  torque  on  a  convex  body  in  free  molecular  flow. 
These  equations  are  not  applicable  to  concave  body  shapes  which  allow 
molecules  reflected  from  one  part  of  the  body  surface  to  impinge  again  on 
the  body. 


3.2.  Aerodynamics  Equations  for  a  Class  of  Spinning  Body  Shapes 

For  spinning  bodies,  the  velocity  of  impinging  molecules  with 
respect  to  the  surface  elements  depends  upon  the  spin  rate  and  position  of 
the  surface  element  with  respect  to  the  spin  axis.  Let  Um  be  the  free 
stream  velocity  and  let  the  spin  rate  of  the  body  be  defined  by  a  spin 
vector,  Q.  The  vector  velocity  of  impinging  molecules  with  respect  to  any 
point  on  the  surface  exposed  to  the  flow  is  then  given  by 


U  =  U  -  fiXF. 


3.7 


where  R  is  the  radius  vector  from  the  center  of  a  coordinate  system  actached 
to  the  body  (see  Figure  3.1). 

The  reflected  velocity  vector,  IK  ,  is  also  dependent  upon  position 

— ♦ 

on  the  surface  and  spin  rate  since  the  magnitude  and  direction  of  IK  is 
related  to  the  magnitude  and  direction  of  the  incident  velocity  U  through  . 
the  parameters  introduced  in  the  generalized  model.  For  a  given  body  shape 
(within  the  class  of  convex  shapes)  the  velocity  vector  U.  can  be  written  in 

"4 

terms  of  the  vectors  U,  n, Q,  and  R  and  the  appropriate  gas-surface  inter¬ 
action  parameters.  For  example,  consider  the  class  of  body  shapes  which 
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have  circular  cross  section  at  any  point  along  the  spin-axis.  For  this 
class  of  body  shapes,  a  i^cal  unit  tangent  vector,  t,  can  be  defined  as 

3 

!qxr| 

where  |QXR|.  represents  the  magnitude  of  the  vector  quantity  QXR. 

•4 

Since  n  is  the  local  unit  normal  vector  to  the  surface  element 
dA,  a  third  unit  vector,  m,  can  be  defined  to  form  a  triad. 


m  =  n 


x  imi 

|qxr| 


3.9 


The  unit  vector  m  is  also  locally  tangent  to  the  surface.  Consider  now  that 

—* 

the  incident  velocity  vector,  U,  has  component  in  the  directions  o£ 

— ♦  —♦  —5 

n,t.m,  respectively.  That  is,  lec 


U  =  U  t  +  iJ  m+U  n 
t  m  n 


Also,  let  U.  be  defined  in  this  coordinate  system  as 


3.10 


U.  =  (U.)_  t  +  (U.)  m  +  (U.)  n 
J  yt  ,  J  m  J  n 


3.11 


If  the  angle  9  represents  the  angle  of  incidence  with  respect  to  the  surface 
plane  (see  Figure  3.2),  the  components  of  U  are  given  by 

U  =  U  cos  0  cos  y 


U  =  U  cos  8  sin  v 
m  ' 


U  =  -U  sin  0 
n 


3.12 


where  U  is  the  magnitude  of  the  incident  velocity  and  y  is  the  azimuth  angle 
of  incidence  defined  in  Figure  3.2. 
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Since  the  reflection  is  in  the  plane  of  U  and  n,  the  components  of 
IK  can  be  written  in  terms  of  and  y. 

(U.)^  =  -U.  cos  0.  cos  y 

J  t  j  j  Y 

(U.)  =  -U.  cos  0.  sin  v  3.13 

J  m  j  J  Y 

(U.)  =  -U.  sin  0. 

J  n  j  J 

where  is  the  magnitude  of  the  reflected  velocity,  U..  The  quantity  U.  is 
related  to  the  quantity  U  by  the  expression  given  in  chapter  2.  That  is, 

U.  =  Jl-a.  U  3.14 

“4 

The  angle  y  is  found  in  terms  of  the  vector  U,  by  the  following  identities 


U  cos  y  = 


U-t 
cos  0 


TT  _  U-M 

U  v  =  r 
sin  '  cos  0 


3.15 


U  - 


sin  0 

Using  Equations  3.15  and  3.14  in  Equations  3.13,  Equation  3.11  becomes 


-COS0 . 


cos0. 


sin0 . 


J  .  OHIO  .  - 

“j  =  ^  c  +  ^  "  ■'  li^e1  <tU>  nJ  3- 


16 


Using  the  definitions  of  t  given  in  Equation  3.8  and  m  given  in  Equation  3.9, 
Equation  3.16  reduces  to  the  following  form 

,-COS0 


U.  =  yi-«j  {-^r1  [u-nxR]  5<r  (|qxr|)‘ 


C0S9-____ 

+  Tosft  [^-nxenxR)]  nX(QXR)  (jfiXRj) 

Sin6i  -  -  -O 

-  — - *•  [U-n]  ni- 

sm0  L  J 


-2 


3.17 
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Equation  3.17  is  valid  for  any  body  shape  which  is  both  convex  and 

has  circular  cross  sections  at  any  point  along  the  spin  axis.  For  shapes  of 

this  class,  the  vector  fiXR  is  always  perpendicular  to  the  unit  normal  vector 
*— ♦ 

n.  Therefore,  using  Equation  3.7,  the  following  can  be  obtained 


U-n  =  (U  -  QXR) •  n  =  U  -n  3. 18 

00  00 

Substituting  Equations  3.18,  3.17,  and  3.7  into  Equation  3.5,  the 
vector  force  on  an  element  of  surface  of  a  spinning  body  having  circular 
cross-sections  about  the  spin  axis  becomes 


dF  =  -p(Uos4n)dA  ju^  -  QXR  +  C^jnxRj)"2  [U^CQXR)  -  (|qxR|)2]  QXR 

.  .3.19 

+  C.(|QXR|)"2  [U  .nX(QXR)]  nX(QXR)  -  S  [U-n]  n) 

J  80  J  " 

where 


COS0  . 

- 1  3.20 

COS0 
_  sin0. 

and  S.  =  o.Jl-a.  — — “•  3.21 

J  3  3  sine 

The  equation  for  torque  can  be  obtained  directly  from  Equation  3.19  by  tak- 
ing  the  cross  product  with  the  vector  R. 

In  the  above  analysis  the  force  and  torque  acting  on  a  spinning 
body  are  expressed  in  terms  of  the  parameters  of  the  generalized  gas-surface 
interaction  model.  For  a  given  body  shape  in  the  class  of  shapes  for  which 
Equation  3.19  is  valid,  the  vectors  U^,  Q,  R,  and  n  are  known.  In  the 
following  sections,  the  force  and  torque  equations  will  be  developed  for 
four  basic  body  shapes  (the  disk,  cylinder,  cona,  and  sphere)  which  will 


C.  =  ojl-a. 
3  3  3 
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illustrate  the  application  of  Equation  3.19.  The  results  will  he  analyzed  in 
the  following  chapter  for  specific  values  of  the  generalized  gas-surface  inter¬ 
action  parameters. 

3.2.1.  Coordinate  Systems 

Two  coordinate  systems  will  be  used  in  the  development  of  equations 
for  the  disk,  cylinder,  cone,  and  sphere.  The  two  coordinate  syscems  which 
are  the  same  for  each  case  are  defined  as  follows. 

Consider  first  a  coordinate  system  (x,y,z)  which  is  attached  to  the 

•*■4 

free  stream  velocity  vector,  U  .  The  positive  y  axis  will  always  be  defined 

to  be  in  the  positive  direction.  Consider  now  a  body  having  a  spin  vector 

— ♦ 

CJ.  Let  a  coordinate  system  (x  ,y  ,z  )  be  defined  in  which  the  z  axis  is 

s  s  s  s 

always  in  the  direction  of  the  spin  vector  Q.  Furthermore,  consider  the  zg 
axis  to  be  initially  inclined  at  an  angle  8g  with  respect  to  the  z  axis  of 
the  x,y,z  coordinate  systems  and  in  the  x-z  plane.  Also  consider  the  xg  axis 
to  be  initially  collinear  with  the  x  axis.  The  coordinate  systems  are  shown 

in  Figure  3.3, 

o 

The  xs,ys,  and  zg  coordinate  system  associated  with  the  body  is  not 
a  co-rotating  system  with  the  body.  Only  the  zg  axis  is  rigidly  attached  to 
the  spinning  body  by  always  being  in  the  direction  of  the  spin  rate  vector  Q. 
The  xg  and  y^  axis  remain  fixed  in  inertial  space  as  the  body  rotates,  if  no 
torques  act  on  the  body.  If  the  body  is  acted  on  by  external  torques,  the 
xg  and  y^  axis  will  rotate  as  a  function  of  the  precession  rate  caused  by 
the  external  torque.  For  cases  in  which  the  free  stream  velocity  vector  is 
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fixed  in  inertial  space,  the  (x,y,z)  system  is  an  inertial  system  and  the 
two  coordinate  systems  are  similar  to  chose  which  would  be  set  up  for  the 

19  -* 

classical  top  problem.  The  case  for  whicn  U  rotates  in  inertial  space 
(such  as  for  a  body  in  orbit)  is  taken  into  consideration  in  a  later  chap¬ 
ter  by  referring  the  results  obtained  using  the  coordinated  systems 
defined  above, by  a  coordinate  transformation,  to  a  true  inertial  reference. 

The  component  of  force  in  the  y  direction  will  be  referred  to  as 
the  drag 'Component  while  the  x  and  z  components  of  force  will  be  termed 
(lift)  and  (lift)  .  In  terms  of  the  components  of  force  on  ar.  element  of 

X  2 

surface,  the  vector  force  is  given  by 


dF  =  d(lift)  i  +  d(drag)  j  +  d(lift)  k  3.22 

X  z 

where  i,j,k  are  the  unit^vectors  along  the  positive  x,y,z  axis.  The  dif¬ 
ferential  notation  is  used  to  denote  the  vector  force  on  the  element  of 
area,  dA. 


The  vector  torque  is  divided  into  components  along  the  xo,ys, 

z  axis.  The  component  of  torque  along  the  z  axis  is  termed  slow-down 
s  s 

torque,  Tg.  The  x^  and  yg  components  of  torque  are  perpendicular  to  the 

•4 

spin  axis  and  are  denoted  by  T.  and  t,  .  The  vector  torque  dT  is  then 

1s  Js 

given  by 


dT  =  d(T.  )  i  +  d(T  )  j  +  d(Tk  )  k 
s  s  s 


3.23 


where  i^.jgjk^  are  the  unit  vectors  along  the  positive  xs>ys>zs  axis. 

The  surface  element  dA  of  the  body  is  at  a  point  defined  by  the 


radius  vector  R  from  the  center  of  the  coordinate  systems.  In  the 


48 


subsequent  development,  cylindrical  coordinates  about  the  z  axis  are  em- 

s 

ployed  in  which  an  angle  §  is  defined  ac  measured  from  the  positive  x  axis 

s 

— ♦ 

to  the  projection  of  R  onto  the  x  -y  plane  (see  Figure  3.4). 

s  s 


3.2.2.  Spinning  Disk  at  Angles  of  Attack 

A  disk  spinning  about  the  surface  normal  is  in  the  class  of  body 
shapes  for  which  Equation  3.19  is  valid.  Consider  a  circular  disk  of  radius 
r^  with  spin  vector  Q  normal  to  the  surface  at  the  center  of  the  disk.  Fig¬ 
ure  3.5  shows  the  disk  with  the  centers  of  the  coordinate  systems,  defined 
above,  at  the  center  of  the  disk.  The  following  vectors  are  required  to 
find  the  force  and  torque  components  acting  on  the  disk.  Referring  to  Fig¬ 
ures  3.3,  3.4,  and  3.5 


— ♦ 

U 

00 


cos  0 


sJs 


U  sin  6  k 
co  s  s 


0  =  a  k 

s 

"4 

R  =  r  cos  |  i  +  r  sin  §  j 
s  s 


3.24 


n  =  k 

s 

where  r  is  the  distance  from  the  center  of  the  disk  to  the  surface  element 
dA  who.  i  dA  is  given  by 

dLA  =  rdrd^  3.25 

The  angle  9  must  now  be  determined  for  use  in  given  by  Equation  3.20 
and  Sj  given  by  Equation  3.21.  Since 


U-n  =  U  sin  9 


3.26 
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-1 


_U  •n~ 


.  .  -  fu-nl  .  -1  f  CD  i 

0  ’ Sln  ltJ  ‘ sin  LirJ 

Using  Equation  3.24,  0  becomes 


3.27 


_  r 


9  =  sin  L  j  ■ 


U  sin  0 

ao _  _ S 


2  2  2  % 
(U  +  r  0  -  2U  rQ  cos  g  cos  0  y 


3.28 


Define  a  non-dimensional  quantity,  K^,  as  the  ratio  of  the  peripheral  speed 
of  the  disk  to  the  free  stream  velocity,  U  .  That  is,  let 


r  ,n 
"a 

CO 


3.29 


where  r^  is  the  radius  of  the  disk.  Define,  also,  a  non-dimensional  variable 


r'  as  the  ratio  of  the  variable  r  to  r^,  That  is,  let 


/  _  r_ 

r 


3.30 


Using  Equations  3.29  and  3.30  in  Equation  3.28,  6  becomes 

sin  0 


0  =  sin 


-i  r. 


L(l+Kd2(r/)2  -  2Kd  r  '  cos  g  cos  g)^‘ 


3.31 


Substituting  Equations  3.24,  3.25,  and  3.30  into  Equation  3.19  for  d?  and 
taking  components  along  the  x,y,z  axis,  the  following  results  are  obtained 

d(drag)  =  Dd  dr  '  dg  j^r'Cl-C^)  -  Kd(r')2  cos  0g  cos  g  (1-C.) 

+  r  '  sin2  0  (C .  +  S  . )"] 
s  J  J  -i 

d(lift  )  =  D i  dr  '  dg  f"-K  (r')2  cos  g  sin  0  (1-C.) 

z  u  L  a  s  j 

-  r  '  sin  0  cc  s  0  (C.  +  S.)  1 
s  s  J  J  J 


3.32 


3.33 
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d ( lif t  )  =  D  dr7  d|  [K,(r7)2  sin  £  (1-C  )] 


where 


2  2 
D,  =  p.U  sin  0  r 
d  v  oo  s  d 


_  cos  0. 

C.  “  - T1 

J  jv  j  cos  9 

_  sin  0. 

S.  =  a  Vl-a .  — : — r1 
J  .1  sin  0 


and  where  0  is  defined  by  Equation  3.31 

Similarly,  the  components  of  torque  are  obtained,  given  by 

d(T.  )  =  D  r  ■  dr 7  d|  [-(r7)2  sin  ^  sin  9  (l+S.)]  ,  3.35 

3  da  s  1 

s 

d(T .  )  =  D  r ,  dr7  dl|  [(r ')2  cos  5  cos  0  (l+S.)]  -  ■  3.36 

Js  s  J 

d(Tw  )  =  D(  rd  dr7  d%  [(r7)2  cos  §  cos  0g  (1-C  )-Kd(r 7)3(1-C . )]  3.37 


where  D,,  C.,  S.,  and  0  are  as  defined  above, 
d’  j’  J5 


Equations  3.31  through  3.37  are  the  basic  equations  expressing  the 
force  and  torque  or.  a  spinning  or  non-spinning  disk,  Tn  chapter  4,  these 
equations  will  be  evaluated  for  specific  cases  of  the  gas-surfa«-e  interaction 
model . 


3.2.3.  Spinning  Cylinder  at  Angles  of  Attack 

The  orientation  of  the  spinning  cylinder  with  respect  to  the  x,y,z 

and  x  ,y  ,z  coordinate  systems  is  shown  in  Figure  3.6.  The  cylinder  of 
s  s  s 

length  L  is  spinning  about  the  axis  of  the  cylinder.  The  center  of  the 
coordinate  system  is  placed  at  the  geometric  center  of  the  cylinder.  The 
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following  vectors  are  defined 
*-♦ 

U  =  U  cos  0  j  -  U  sin  0  k 
oo  oo  s  s  co  s  s 

S  =  Qk 

s 

3.38 

R  =  r  cos  E  i  +  r  sin  |  j  +  j£k 
cy  *  s  cy  3  s  s 

n  =  cos  §  i  +  sin  §  j 
s  s 

where  r  is  the  radius  of  the  cylinder  and  Z  is  the  variable  o£  integra¬ 
tion  along  the  axis  of  the  cylinder.  The  element  of  surface  dA  is  given 
by 

dA  =  r  d£  dZ  3.39 

cy 

From  Equations  3.38,  the  angle  0  is  found  to  be 


0  =  sin 


■cos  0  sin  £ 
s 


2  h 

(lvK  -  2K  cos  F  cos  0  y 
cy  cy  ’  s 


3.40 


where  the  non-dimensional  spin-rate  parameter  K  is  defined  as 


K 

cy 


r  Q 

c.Y 

U 


3.41 


The  variable  of  integration  Z  appears  only  in  Equation  3.39  above. 
The  equations  for  the  cylinder  can  then  immediately  be  integrated  with 
respect  to  Z  from  -  j  to  +  The  variable  %  still  remains,  however. 

Letting 

2 

D  =  pU  r  L  cos  0  , 
cy  v  05  cy  a 

The  resulting  basic  equations  for  the  cylinder  are  found  to  be 
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d(drag)  =  D  d?!  -sin  E  (1-C.)  +  K  cos  8  sin  E  cos  |  (1-C.) 

cy  J  cy  s  j 

-  cos20  sin2?  (C.  +  S.)l 

s  *  J  J  J 

d(liff)  =  D  d?  I  K  sin  0  cos  ?  sin  ?  (1-C.) 
z  cy  5  L  cy  s  *  J 

-  cos  0  sin  0  sin2?  (C .  +  S.)~l 

s  s  j  j  J 

sin2^  (1-C  ) 

9  J 

ns  0  cos  §  sin^E  (C.  +  S.)"| 

s  J  .1  _j 

sin  0  sin2?  (l-C.)l 
s  3  j  j 

-\ 

sin-  8g  cos  ?  sin  ?  (1  -  C.)  j 

K  sin  ?  (1  -  C.) 

-  cy  *  r 

+  cos  8g  cos  |  sin  ?  (1-C.)J 


=  D 

M 

-K 

.  £■. 
sm  1 

cy 

-  cy 

-  cos  8  i 

s 

f 

- 

=  D  . 
cy 

r 

.cy 

dE  ^sin  6 

=  D 

r 

sin-  0 

cy 

cy 

=  D 

r 

1 K  s 

cy 

cy 

L  cy 

3.42 


3.43 


3.44 

3.45 

3.46 

3.47 


where  and  are  as  defined  previously. 

The  surface  of  the  cylinder  exposed  to  the  flow  is  from  ?  =  tt  to 
2rr  at  all  angles  of  attack  except  8  =  ?  and  -  ?  where  none  of  the  surface 
is  in  the  flow.  Equations  3.42  -  3.47  are  the  basic  expressions  for  force 
and  torque  on  a  spinning  or  non-spinning  cylinder.  These  equations  will  be 
evaluated  in  chapter  4  for  specific  cases  of  the  reflection  model. 


3.2..4,.  .Spinning  Cone  at  Angles  of  Attack 

A  cone  having  half  angle  6,  height  Hc,  and  base  radius  r^,  is 

shown  oriented  in  the  x,y,z  and  x  ,y  ,z  ,  coordinate  systems  in  Figure  3.7. 

s  s  s 

The  center  of  the  coordinate  syscem  is  at  the  center  of  the  base  of  the 
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cone.  The  cpne  is  spinning  about  the  axis  of  the  cone.  The  following  vec¬ 
tors  are  defined 


— ► 

U 

CO 


U  cos  0 

CO  ~  s 


n  - 


k 

s 


U  sin  0  k 
oo  s  s 


R  =  f(A)  cos  §  i  +  f(j 0  sin  |  j  +  S,  k 
.  s  s  s 


n  =  cos  E  cos  6  i  +  sin  §  cos  6  j  +  sin  6  k  3.48 

3  s  s  s 

where  A  is  the  variable  of  integration  along  the  axis  of  the  cone  and  f(A) 

is  the  radial  distance  of  the  surface  dA  from  the  axis  of  the  cone.  The 

distance  f(jj)  is  given  by 


r'U) 


r 

c 


The  element  of  surface  dA  is  given  by 


dA  -  d§  dj e 

COS  0 


The  angle  0  is  found,  using  Equation  3.48,  to  be  given  by 


sin  0  sin  6  -  cos  0  sin  £  cos  6 

sin  0  - - r - « - - - r  3.49 

(1  +  (A  'j  Kc  -  2j&'  Kc  cos  |  cos  0  V5 


where  Kc  is  the  non-dimensional  spin  rate  parameter  defined  as 


CO 


and  A '  is  a  non-dimensional  variable  defined  as 


A 

H 


i'  =  1 
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3 


il 


I  "? 


3 

-* 

;! 

“f 


f  H 


3 


Letting 


_  _  .2  Ic _ 

Dc  P  oo  sin  6 


The  basic  equations  for  the  cone  become 

d(drag)  =  D  dl|  d J),'  "  C.)  +  XK  (<&7)2  cos  g  cos  0  (1  -  C.) 

c  L  J  c  s  j 

-X3  i'  (C.  +  S.)]  3.50 

d(lift)  =  D  d?  d i,'  [xK  U7)2  cos  g  sin  9  (1-C.) 

z  c  L  c  s  j 


2 

~X  (sin  6  cos  0 


+  sin  g  cos  6  sin  9  )(C.  +  S.)  1 
*  s  j  j  J 

d(li£t)x  =  Dc  dg  dz'  [-XKC  W)2  sin  %  (J.  - 

-X2  cos  g  cos  6  Jt(  (C.  +  S  )] 

d(T.  );=  Dcr  d?  [x(-e/)2  sin  es  sin  5  (1_c-) 

O  C  ’ —  J 


3.51 


3.52 


+  y,Z '  (1-J&7)  cot  6  cos  0  (1-C.) 

s  J 

-  XK  (Z1)2  (l-Z  )  cot  6  cos  g  (1-C  ) 

C  J 

+  X2{“(X7)2  sin  6  sin  g 

+  Z'(.l-Z')  cot  6  sin§  cos  6}(C..  +  S^j  3.53 


dC 


[  )  =  Dcrc  dg  d£'  I  -xKc  a,)2(l-je')  cot  6  sin  5  (1-C.) 

Jn  L  J 


-XU')2  sin  0  cos  g  (1-C.) 

S  J 

2 

~X  £7(i  “  £')  cot  6  cos  §  cos  6  (C.  +  S.) 

3  J 

+X.  (.Si'y  sin  6  cos  §  (Cj  +  S4)  J 


3.54 
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d(Tk  )  =  Dcrc  dg  di'  [(XKc  a')3  ~  XU')2  cos  0s  cos  g)  (1-C  )J  3.55 


where 


X  =  cos  8  sin  g  cos  6  -  sin  0  sin  6 
s  s 


3.56 


The  surface  of  the  cone  exposed  to  the  flow  is  a  function  of  the 
angles  0  and  6.  At  8  =  ?,  for  example,  the  entire  conical  surface  is  ex- 

S  S  4 

posed  to  the  flow  (g  =  0  ->  2rr)  .  In  fact,  for  all  values  of  0g  between  ^  and 
^  -  6,  the  entire  conical  surface  is  exposed  to  the  flow.  At  angles  of  0g 
between  ^  +  6^),  decreasing  amounts  of  the  conical  surface  is  exposed  to 
the  flow.  The  shadow  boundary  for  the  flow  on  the  surface  is  always  a 
straight  line  and  can  be  found  in  terms  of  the  variable,  g,  by  finding  where 
U*n  =  0.  Define  the  angle,  p,  to  be  the.  value  of  g  at  the  shadow  boundary. 
Then, 

P  =  sin  ^  (tan  0g  tan  6)  3.57 

EJrom  Equation  3.57,  the  following  three  categories  of  surface  exposure  are 
defined. 

Case  I:  tan  0  tan  8=1 
■s 

Entire  conital  surface  (|  =  0  -*  2rr)  is  exposed  to  the 
flow. 

Case  II:  -1  <•  tan  0  tan  6  <  1 
—  s  — 

'  Only  part  of  the  conical  surface  - 

g  =  (tt  -  P)  -*  (2tt  +  P) 

is  exposed  to  the  flow. 

.Case  III:  tan  0  tan  6  <  -  1 
t  s  — 


.Case  III: 
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No  part  of  the  conical  surface  is  exposed  to  the  flow. 
(The  base  of  the  cone  having  not  been  included  in  the 
basic  equations  for  the  cone). 


3.2.5.  Spinning  Sphere  at  Angles  of  Attack 

A  sphere  of  radius  rg  is  shown  oriented  with  respect  to  the  x,y,z 

and  x  ,y  ,z  coordinate  systems  in  Figure  3.8.  The  spin  axis  is  along  a 
s  s  s 

diameter  of  the  sphere  at  an  angle  (“-0  )  with  respect  to  the  free  stream 

z  s 

velocity.  Cylindrical  coordinates  are  used  to  describe  the  position  of  the 
surface  element  dA.  With  i,  denoting  the  variable  length  along  the  spin  axis, 
dA  is  given  by 

dA  »  r  d§  U 
s 


The  following  vectors  are  defined 


U  =  U  cos  9  j  -  U  sin  9  k 

eo  os  s  s  CO  S  S 


o  =  a  k 


R  =  r(£)  cos  §  i  +  r(X)  sin  £,  +  f.k  3.58 

s  Js  s 

-*  r(4)  „  .  .  r(£)  .  „  .  .  S,  , 

n  =  cos  £  l  +  — 1 u  sm  c  j  -1 - k 

r  ^  s  r  OJsrs 

s  s  s 

wnere  r(j g)  is  the  radial  distance  frcm  the  spin  axis  at  the  point  &  to  the 
element  of  surface  dA 


r(jg)  =  (x  1  -  sF) 


v  s  y 


3-.  55 


The  angle  9  is  obtained  from  Equations  3.58  as 

i,'  sin  9  -  St,  sin  £  cos  0 

C  C  J  < 


sin  0  = 


(1+X  2  K  2  -  2K  &  cos  £  cos  0  F 
s  s  s  s  s 


3.60 
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I 

1 

I 

I 


where  Kg  is  tf'e  non-dimensional  spin-rate  parameter  defined  as 

r  Q 

K  = 
s  U 

CO 

and  Z'  is  a  non-dimensional  variable  of  integration  defined  as 


s 


I 


I 


and 


Letting 


A’  =  [1-U')2]% 


2  2 

D  =  p  u  r 
S  y  co  s 


The  basic  equation  for  the  sphere  becomes 


d(drag)  =  D  d§  &Z  '  [x  (-1  +  K  Z  cos.|  cos  0  )(1-C.) 

s  s  s  s  S  J 


xs  (Cj+Sj)] 


d(lif t  ) 
z 


D  d5  [X.kJ.  cos  §  sin  0„  (1-0.) 

S  S  b  S  ©  J 


X  1  (Z*  cos  5  +  ft  sin  fc*  ;o.  0  )  (C.+S, 

s  s  ?  s  j  ; 


d(lift  ) 
x 


D£  d?  d Z'  C-Xs  \  &s  sin  |  (1-C<) 


I 


d(T.  ) 
s 


-  X  *  cos  %  (C.+S,)] 
^s  s  jj 

Ds  rs  iA'  &S  sin  0s  sln  5 


-  K  z'z  cos  I  +  z'  cos  0  1  X  (1-C 

S  S  S  •  s 


3.61 


)]  3.62 


3.63 


.)  3.64 
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d(T.  )  =  D  r  dE  dZ  [-Z  sio  0  cos  E  -  K  a  SL  sin  El 'X  ( 1— C . )  3.65 

J  SS.  US  S  5  S  S  S  I 

Js  . 

d(T  )  =  D  r  d|  dz'  [K  i  2  -  l  cos  E  cos  0  ]  Xc  (1-C.)  3.66 

k  ss  uss  s  s  s  i 

s  J 

where 

y.  -  Z  cos  0  sin  E  -  Z'  sin  0 
s  s  s  s 

The  shadow  line  of  the  flow  on  the  surface  of  the  sphere  is  a 

curve  in  the  Z'  -  %  plane.  At  a  constant  Z' ,  the  shadow  limits  in  terms  of 

E  can  be  found  by  finding  where  U*n  =  Q.  Define  the  angle  to  be  the  value 

of  E  at  the  shadow  boundary.  Then,  3  is  given  by 

s 

Ps  =  Sin"1  (tan  0g  j~  ) 

s 

For  arbitrary  values  of  0  ,  the  shadow  determines  three  regions  of  exposure 

s 

along  the  spin  axis.  Referring  to  Figure  3.9 

o  ' 

Region  I:  tan  9  „  >  1 

s  Zs 

In  this  region  the  surface  from  E  =  0  to  E  =  2tt  is  exposed  to  the 


flow 

Region  II:.  -1  <  tan  0  <  \ 

— -  "  '  s  ” 

In  this  region  the  shadow  limits  are,  at  any  point,  z' ,  from 
E  =  tt  -  Pg.to  E  =  2tt  +  pg. 

Z ' 

Region  III:  tan  0^  ~~  <  -1 

"  S  Xj 


This  region  of  the  sphere lis  not  exposed  to  the  flow. 


>  —*> 


mb 
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Figure  3.9.  Regions  of  sphere  surface  exposed  to  the  flow. 
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3.3.  Applications 

The  equations  which  have  been  obtained  for  the  four  basic  shapes 
of  the  disk,  cylinder,  cone,  and  sphere  are  applicable  to  a  wide  variety  of 
problems.  In  the  next  chapters  these  equations  will  be  used  to  evaluate  the 
effects  of  spin  rate  and  the  gas  surface  interaction  on  the  dynamic  pro¬ 
perties  of  these  shapes  in  free  molecular  flow.  The  application  of  the 
equations  to  other  problems  is  discussed  as  follows. 

Many  satellite  shapes  can  be  simulated  by  a  combination  of  the 
four  basic  shapes  of  the  disk,  cylinder,  cone,  and  sphere.  If  the  satellite 
is  of  convex  shape,  the  basic  equations  developed  in  this  chapter  can  be 
applied  directly.  Since  the  equations,  are  in  differential  form,  they  are 
equally  applicable  to  bodies  composed  of  segments  of  the  sphere  (such  as 
spherical  caps),  cone  frustrums  and  segments  of  the  cylinder  and  disk. 

Surface  properties  affect  the  character  of  the  gas-surface  inter¬ 
action.  Many  satellites  are  composed  of  surfaces  which  have  widely  differing 
surface  properties  such  as  solar  cells  versus  painted  surfaces  and  varying 
surface  temperature  or  roughness.  Problems  of  this  type  can  be  studied 
using  the  equations  developed  in  this  chapter  by  assigning  different  values 
to  the  parameters  of  the  gas-surface  interaction  model  for  a  specific  region 


of  the  surface. 
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4.  AERODYNAMIC  PROPERTIES  OF  SELECTED  SHAPES 

In  this  chapter  aerodynamic  properties  of  selected  shapes  are 
evaluated  from  the  expressions  developed  in  chapter  3.  For  the  special  case 
of  the  Maxwell  model,  results  are  obtained  to  illustrate  the  effects  of  spin 
on  the  aerodynamic  properties  of  the  four  basic  shapes  of  disk,  cylinder, 
cone,  and  sphere.  Analytical  results  are  also  obtained  in  terms  of  the 
generalized  model  parameters  for  drag  of  the  four  body  shapes  with  zero  spin. 

In  general,  the  aerodynamic  properties  must  be  evaluated  by  using 
numerical  methods.  The  numerical  techniques  employed  in  this  study  are  out¬ 
lined  in  this  chapter  and  results  obtained  for  various  values  of  the  gen¬ 
eralized  model  parameters  are  presented.  These  results  are  compared  with 
those  obtained  for  the  special  case  of  the  Maxwell  model. 

In  the  next  chapter,  the  results  obtained  for  the  cone  and  disk 
are  applied  to  the  problem  of  a  satellite  in  near  earth  orbit. 


4.1.  Special  Case  of  the  Maxwell  Model 

The  reduction  of  the  generalized  gas-surface  interaction  model 

parameters  to  the  parameters  of  the  Maxwell  model  was  demonstrated  in 

section  2.3.1.  Those  results  can  be  applied  to  determining  in  terms  of 

the  Maxwell  model  parameters  the  quantities  C.  and  S.  which  appear  in  the 

J  J 

■rat 

basic  aerodynamic  equations  developed  in  chapter  3.  For  specular  reflec¬ 


tion, 


C. 

J 


_  cos  0. 

jjl-a.  - r1 

J*  J  cos  e 


S.  =  CTV’i  -e.'. 

J  J  J 


sin  0. 


“d 
1  -  ad 


4.1 
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For  diffuse  reflection. 


C.  =  0 
J 


S.  -  0.^1. 


■a, 


T  sin  6 


Combining  diffuse  and  specular  reflection,  the  values  of  and  S.  in  terms 
of  the  parameters  of  the  Maxwell  model  become 
(C,)  »  1  -  Of, 


L.  2 


(3  )  -  1  •  ad  +  Vl-«T  7i“l 

m 


In  the  basic  aerodynamic  equations,  the  follow.' ng  functions  of  and 
occur  repeatedly.  In  terras  of  the  Maxwe.’l  model  parameters,  these  functions 
are  given  as 

i  -  (v  -  «4 

m 


4.3 


(C  )  +  (S  )  -  2(l-„d)  +  STi 

m  J  m 

1  +  (S.)  -  2  -  or ,  +  or.  JbZT,  --1  t 

v  y  d  d  v  T  sin  0 

m 

The  equations  for  the  disk,  cylinder  and  cone  could  be  integrated 

.  i 

analytically  except  for  the  terms  containing  (sin  0)  fhich  appear  in  Equa¬ 
tions  4.3. 


The  angle  of  incidence  of  the  flow,  0,  is  a  function  of  the  spin- 
rate  and  position  of  the  surface  element  dA  as  discussed  in  chapter  3.  For 
the  four  body  shapes  studies,  the  expressions  found  in  chapter  3  for  sin  0 
for  each  shape  are  all  of  the  following  form 
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sin  0  =  A  [1  +  K2d2  -  2Kd  cos  §  cos  0  ]  ^  4.4 

where  Table  4.1  gives  the  values  for  A,  K,  and  d  for  the  four  shapes  studied 
in  chapter  3.  The  terms  K  and  d  are  non-dimensional  where  d  represents, 
except  of  the  case  of  the  cylinder,  a  variable  of  integration  which  may  have 
a  maximum  value  of  one  or  less.  The  spin  rate  parameter  K  is  the  ratio  of 
the  maximum  peripheral  velocity  of  the  body  to  the  free  stream  velocity  U^. 
For  satellite  applications,  K  is  always  less  than  one  as  can  be  determined 
from  the  values  for  K  given  in  Table  4.?  for  a  number  of  past  satellites. 
Since  K  <  1  and  d  <  1  for  satellite  applications, 

(K2d2  -  2Kd  cos  Gs  cos  §)2  <  1  4.5 

and  (sin  0)  ^  can  be  expanded  in  a  binomial  series  as  an  approximation  to 
the  dependence  of  9  on  spin  rate  and  surface  position.  The  result  of  the 
expansion  is 


(sin  0)  ^  [1  -  Kd  cos  0  cos  §  +  ~  (Kd)2  (1  -  cos20  cos2  §) 

A  s  Z.  S 

+  O(Kd)3] 


4.6 


3 

where  terms  of  order  (Kd)'  and  higher  are  not  to  be  retained.  Substitut¬ 


ing  Equation  4.6  into  Equation  4.2  for  and  S^  of  the  Maxwell  model,  the 
following  results  are  obtained 


<Cj>  =  1  '  *d 
m 

(S.)  =  1  -  o'  4  aJl-a  A  ^[1-Kd  cos  0  cos  §  4.7 

j  a.  u  X  s 

m 

+  %  (Kd)2  (1-COS2  0g  COS2  5)] 


Equations  4.7  will  now  be  substituted  into  the  basic  aerodynamics  equations 
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Table  4.1. 

Definitions  of  A,  K,  and  d  obtained  from  the 
results  of  chapter  3 


Body 

Shape _ A _ K _ d 


disk 

cylinder 

cone 

sphere 


sin  0  K, 

s  d 

"  cos  8  sin  f  K 

s  *  cy 

sin  8  sin  6  -  cos  8  sin  £  cos  6  K 
s  s  c 

Sj'  sin  8  -  j &  sin  §  cos  8  K 

s  s  s  s 


rd^ 

U 


r  Q 

cv 


u 


r  O 
c 

U 


r  Q 
s 

U 


t' 


■  a2  ^ 
(i- 

r 

s 


Table  4.2. 


Values  of  the  spin  rate  parameter,  K,  for 
a  number  of  past  satellites 


Satellite 


Orbit  Altitude 


20 


Spin  rate' 


20 


K 


perigee 

- 

apogee  (mi) 

(rpm) 

(10“3) 

ESRO  2A 

215 

680 

35 

.350  -  .519 

ESRO  2B 

205 

- 

677 

35 

.151  -  .223 

LDSCS  19-26 

20,940 

- 

21,068 

150 

2.29 

■ERS  16 

110 

- 

2,260 

10 

.01-97  -  .0301 

ERS  18  and  0V5-1 

5,357 

- 

69,316 

10 

.0306  -  .238 

OSO  3 

336 

- 

354 

30 

.230 

Explorer  32 

173 

- 

1,629 

30 

.169  -  .227 

ESSA  1 

432 

- 

521 

10 

.074  -  .079 

ESSA  2 

843 

- 

885 

10.9 

.084  -  .086 

ESSA  4 

822 

- 

894 

10 

.078 

ESSA  6 

876 

- 

925 

10 

.078 

obtained  in  chapter  3.  For  the  disk,  cylinder,  and  cone,  analytical  results 
can  be  obtained.  The  equations  of  chapter  3  expressing  force  and  torque  on 
a  spinning  sphere  are  not  integrabie  in  closed  form  when  Equations  4.7  are 
employed;  however,  the  torque  on  a  spinning  sphere  has  been  obtained  analyt¬ 
ically  in  terms  of  the  Maxwell  model  parameters  by  using  Hie  techniques 
employed  in  Reference  9,  The  results  from  Reference  9  for  the  sphere  will 
be  given  along  with  the  resulcs  for  drag  of  a  non-spinning  sphere  which  can 
be  bttined  analytically  from  the  equations  given  in  chapter  3. 

The  results  will  be  presented  in  terms  of  coefficients  of  drag, 
lift,  and  torque,  which  are  defined  as  follows. 


Drag  coefficient  = 

Lift  coefficient  = 


%PU  2a 
00 

Lift 

%pU  2a 

CO 


Torque  coefficient  =  C„, 


where  A  is  a  reference  area  associated  with  the  particular  shape  being  .  .. 
studied.  In  addition,  a  reference  radius  r  is  introduced  in  the  expressions 
for  coefficients  of  torque. 


4«1,1«  Spinning  Disk  Properties  in  Terms  of  The  Maxwell  Model  Parameters 


Substituting  Equations  4.7  for  the  Maxwell  model  into  Equations 


3.32  through  3.37  for  the  spinning  disk  at  angles  of  attack  and  integrating 
over  r'  and  5  (0  to  1  for  r'  and  0  to  2rr  for  |)  the  following  results  are 


71 


2  - 


obtained  for  the  case  of  the  spinning  disk,  (A  =  tt  ,  r  =  rd)) 

,2  „  \ 


z  , 


{2 

.  2 
sxn 

e 

s 

+  Of. 
c 

L  (1 

.  - 

2 

s: 

K 

+ 

v1- 

•Of T 

sin 

9 

s 

[1 

+ 

( 

l 

[-2 

sin 

9 

s 

cos 

e„ 

cr 

+ 

2  0fd 

K 

- 

v1- 

-t" 

cos 

9 

s 

[1 

+ 

( 

t 

.  2 


.  2 


(ad7l-aT  Kd  sin  es  cos  eg)/2 

Of,  K,  sin  9 
d  d  s 


4.8 


4.9 

4.10 

4.11 

4.12 

4.13 


A  spinning  disk  is  found  to  have  drag  and  (lift)z  coefficients 

which  are  higher  than  a  disk  with  zero  spin  rate.  This  increase  in  drag  is 
2 

proportional  to  which  would  be  small  for  satellites. 

Besides  the  expected  slowdown  torque,  a  spin  induced  precession 
torque  about  the  y  -axis  is  obtained  which  is  directly  proportioned  to  K^. 
The  precession  torque  arises  because  diffusely  reflected  modules  are  re¬ 
flected  at  a  velocity  relative  to  the  surface  which  is  higher  on  one  half 

TT 

of  the  disk  when  the  disk  has  angles  of  attack  other  than  ^  and  zero. 


4,1.2  Spinning  Cylinder  Properties  in  Terms  of  Maxwell  Model  Parameters 
The  coefficients  of  drag,  lift  and  torque  for  the  spinning 


cylinder  at  angles  of  attack  are  obtained  by  first  substituting  Equations 
4.7  for  the  Maxwell  model  into  Equation  3.42  through  3.47  for  the  cylinder. 
These  equations  are  then  integrated  over  the  surface  area  exposed  to  the  flow 
with  respect  to  §  from  n  to  2n.  Letting  A  =  2  r  L  and  r  =  r  the  coeffi¬ 
cient  of  drag,  lift,  and  torque  for  the  spinning  cylinder  at  ingles  of  attack 
are  found  to  be 

C_  =  cos  0  {-f  cos2  0  +  O' ,  ■§  (3/4  -  cos2  0  ) 

D  s  s  a  3  s 

_  K  2 

+  ?  ad/^~aT  COS  ®s  ^  ^  ”  C0S^  4.14 

8  8 

CT  =  cos  0  cos  0  sin  0  -  a ,  ~  cos  0  sin  0 

L  s  3  s  s  a  3  s  s 

Z  _  K  2 

+  ja^/l-aT  sin  0s  [1  +  (4  -  cos2  0g)]}  4. 15 

C  =  cos  0  K  [-  a ,  +  \  O'  cos  0  ]  4.16 

L  scyL2d3d^  T  sJ 


C  =  -x  a ,  sm  0  cos  0  4.17 

T.  2d  s  s 

i 

s 


CT  =0  4.18 

i 

s 


C_  =  -  2  K  a,  cos  0  4.19 

Tk  cy  d  s 

s 


The  results  obtained  for  the  cylinder  exhibit  the  same  influence  of 
spin  rate  on  drag,  lift  and  slowdown  torque  as  those  of  the  disk.  Unlike 
the  spinning  disk,  the  spinning  cylinder  experiences  a  spin  induced  lift  in 
the  negative  x  direction  which  is  directly  proportional  to  the  spin  rate  • 


parameter  K 
r  cy 


This  spin  induced  lateral  force  on  the  cylinder  is  analogous 
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to  a  similar  effect  called  the  Magnus  force  which  occurs  on  spinning  cyl¬ 
inders  in  continuous  flow  (see  Reference  21) »  In  free  molecular  flew,  the 
lateral  force  is  found,  however,  to  be  in  the  direction  opposite  that  of 
the  analogous  Magnus  force. 

The  finite  value  of  torque  perpendicular  to  the  spin  axis  obtained 
in  Equation  4.17  is  due  entirely  to  geometry  and  is  not  spin  induced  as  was 
found  for  the  disk.  The  center  of  the  surface  area  exposed  to  the  flow  is 
not  at  the  geometric  center  of  the  cylinder  but  is  a  function  of  the  angle 
of  attack.  If  end  contributions  were  included  in  the  development  of  equa¬ 
tions,  CT  would  be  zero. 

'i 

s 

4.1.3.  Spinning  Gone  Properties  in  Terms  of  Maxwell  Model  Parameters 

The  coef'  :ients  of  drag,  lift,  and  torque  for  the  spinning  cone 
at  angles  of  attac  .  .  obtained  by  substituting  Equations  4.7  for  thp  Max¬ 
well  model  into  Equations  3.30  through  3.56  for  the  cone.  These  equations 
can  be  integrated  in  general  form  by  incorporating  the  angle  f3  and  integrat¬ 
ing  over  the  surface  area  exposed  to  the  flow  with  respect  to  X  '  from  0  to  1 

-  2  — 

'and  5  from  (tt  -  {3)  to  (2tt  +  (3)  .  Letting  A  =  tt  r  and  r  =  r  the  coefficient 

c  c 

of  drag,  lift,  and  torque  for  the  spinning  cone  at  angles  of  attack  are  found 
to  be 

=  R„  +  a,  \  —  cos  0  cot  6  cos  (3  +  sin  0  -  R.  j 

D  2  d  Ltt  s  it  s  2  J 

_  r  K  2  2 

+  V1_0rT  l3  +  “4"  (R3  “  R8  COS  Vj 


4.20 
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CT  -  -(1  -  a.j  (2R  sin  6  cos  8  +  2R,  cot  6  sm  8  ) 

ii  a  3  si  s 

z 

+  [~  R^  cos  8s  +  R^  cot  6  sin  9 


Kc  f 

+  ~r~  cR,  cos  9  +  R_  cos  6  sin  0 

4  L  4  s  5  s. 


2  A 
cos  8 


scr7 + vi: 


-Kc  ^  (~ad  (cos  8g  cot  6 


3-sm2S  .  4 

- *”  +  —  sm  9  cos  S) 

TT  TT  s 


+  2  ar. ./I. -or-  B.-  ;js  9  cot  6 
a'  i  /  s  w 


>  1-3  sin^S 

1  ,  .  2e 
3  sm  6 


+  I  sin2  6  cos  f}  -  —  cos  8  sin  8  cot  6  sin  26 
j  Lit  s  ttss 

-  -  cot2  6  cos2  8  cos  6  -  2R  ~6'  1 

tt  s  1  .  2,  ■ 

sm  6 

-  [s  Qtt& 


,  .  2. 
3  sm  6 


+  i-  ®3  -  ~2  v  C-f^ir)] 


c  r  .  2  e 

Y2  [^(cot  6  cos  8£ 


5  sin  d 


+  —  cot  6  sin  8  cos  p) 

TT  S 


+  ar/l-aT  2R?  cos  eg  -1'4  *in  6j 

sin  6 

-  K  or  (-  cos  0  cot  6  cos  6  +  II^  sin  8  ) 
c  d  tt  s  2tt  s 


The  functions  R^,  ,  -  Rg  used  in  Equations  4.20  through  4.25 
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are  defined  as  follows. 


t—  L. 

7“  cos 

8  cos 

3tt 

s 

2  .  2 

_  .  2 

~  sm 

0  sm 

TT 

s 

sin  0 

cos  0 

TT 


4  3  2  2 

R„  =  “  cos  8  cos  6  cot  6  cos  8  (sin  B  +  2) 
2  3tt  s 

,  „  .  3  .  .  2  .  tt+26 

+  2  sin  0  sm  6  - ^ 

s  TT 

.  .  .  r  2  2.  tt+2B  -  sin2B 

+  3  sin  0  cos  0  cos  6 

s  s  TT 

12  2 

H - sin  9  cos  0  sin  6  cos  6  cos  B 

tt  s  s 


„  2  „  „  .  „  tt+2B  -  sm2B 

R_  =  cos  8  cos  6  cot  6  - - K 

3  s  2tt 

+  —  sin  8  cos  0  cos  6  cos  B 
tt  s  s 

.  2  .  .  tt+28 

+  sm  0  sin  o  - 

s  TT 


2.  tt+28 

R,  =  -  —  cos  9  cos  6  cos  B  -  sin  0  sin  6  — — 
4  tt  s  s  TT 


„  .  .  tt+2B  -  sm2B  .  2  . 

R_  =  cos  0  cos  6  -  n -  H —  sin  0  sm  o  cos  p 

5  s.  2tt  tt  s 


R.  =  cos 
6  s 


„  2tt+48  -  sin48  .2  .  .  3  „ 

0  cos  6  - - H  —  sm  0  sm  6  cos  B 

16tt  3tt  s 


2  _  ,  o  „  tt+2B  +  sin2B 

R_  ®  -  —  cos  0  cos  o  cos  B  -  sm  0  sm  6  -  ^ - 

7  3tt  s  s  2tt 


_  2  „  r  2tt+4B  -  sm  4B 

R.  =  cos  0  cot  6  cos  5  - - 

8  s  16tt 

+  ~  sin  0  cos  6  cos  6  cos"^  B 
3tt  s  s 

.  9  .  .  tt+2B  +  sin28 

+  sin*1  0  sm  6  - K 

s  2it 
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In  Equations  4.20  through  4.25,  the  value  of  the  angle  g  depends 
upon  the  angles  9g  and  6  in  accordance  with  the  three  cases  of  surface  ex¬ 
posure  defined  in  the  last  chapter  for  the  cone.  Using  those  definitions 
the  following  ranges  of  g  are  defined. 

Case  I:  tan  6  tan  6  >  1 

s  “ 


e  - 


TT 


Case  II:  -1  <  tan  0  tan  6  <  1 
_  s  — 

g  =  sin  (tan  0^  tan  6) 


4.26 


Case  III:  tan  0  tan  6  <  -  1 
s  ~ 

o  =  .  S 
P  2 

For  case  III,  the  coefficients  of  drag,  lift,  and  torque  all  become  zero 
because  the  conical  surface  is  shaded  from  the  flow. 

The  cone  equations  are  found  to  be  similar  in  form  to  those  ob¬ 
tained  for  the  cylinder  and  disk.  As  in  the  case  of  the  spinning  cylinder, 
a  spin  induced  lateral  force  is  experienced  by  the  spinning  cone  at  angles 
of  attack.  The  spin  induced  lateral  force  is  found  to  be  directly  propor¬ 
tional  to  the  spin  rate  parameter  K£. 

The  cone  is  found  to  also  experience  a  spin  induced  torque,  per¬ 
pendicular  to  the  spin-axis,  about  the  y  axis  of  the  cone.  The  component 

s 

of  torque  about  the  x  axis  is  not  spin  induced  but  is  instead  due  to  the,. 

s  V 

moment  of  the  drag  and  (lift)  forces  about  the  center  of  the  coordinate 

z 

systems. 


4.1.4.  Sp-'nning  Sphere  Properties  in  Teims  of  the  Maxwell  Model  Parameters 
The  sphere  equations  given  in  the  last  chapter  must  in  general  be 


0 
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evaluated  by  numerical  techniques  for  all  cases  of  the  gas^surface  inter¬ 
action  model.  The  equations  given  in  chapter  3  are  particularly  useful 

it,  ">• 
x 

when  applied  to  shapes  such  as  spinning  spherical  caps,  and  hemispheres. 

For  the  special  case  of  a  complete  sphere,  analytical  solutions  for  the 

torque  components  are  given  in  Reference  9.  The  results  of  Reference  9 

will  be  given  here  and  are  the  same  is  would  be  obtained  by  numerically 

solving  Equations  3.61  through  3.66. 

~  2  - 

'  Le’:t‘.ng  A  =  rt  r  and  r  =  r  ,  the  coefficients  of  torque  on  the 

•  S  S  « 

.sphere  are  as  obtained  in  Reference  9, 


s 


CT  =  (Kg  Qfd  sin  0s  cos  9g)/2 


C_,  =  -  K  o'  (2  +  cos2  6  )/2 

T.  s  d  s 


4.27 

4.28 

4.29 


>  s 

The  coefficient  of  drag  for  a  non- spinning  sphere,  in  terms  of  the  Max¬ 
well  model  parameters,  can  be  obtained  analytically  from  the  equations 

-  2 

given  in  chapter  3.  Letting  A  =  rr  . 


°d  =  2  +  I  V--“r  Ks  =  °  4-30 

The  coefficients  of  drag  and  lift  for  the  spinning  sphere  must  be  obtained 
by  numerical  integration  techniques  and  are  therefore  not.  available  in 
analytic  form. 


4.2.  Solutions  in  Terms  of  the  Generalized  Model  Parameters 
4.2.1.  Analytical  Results  (Zero  Spin) 
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For  the  special  case  of  zero  spin  rate,  the  equations  of  chapter  3 
can  be  integrated  in  closed  form  for  certain  cases  of  angle  of  attack. 

These  results  are  given  in  the  following  sections. 


4. 2. 1.1.  Drag  and  Lift  of  Flat  Plate  at  Angles  of  Attack 


The  angle  0  for  the  disk  (Equation  3.31)  becomes  0g  for  the  case 


when  =  0.  Therefore, 


C.  «  O'  .71 -or . 
J  J  1 


-  (t*i +  (i-v 0 


cos  0 


_  sin  (5  P.  +  (1-P.)  9  ") 

S.  =  o  .Jl-a.  - -  1— 

-  1  ’  sm  0 


J  J  J 


which  are  independent  of  x'  and  §. 

Substituting  Equations  4.31  and  4.32  into  Equations  3.32  through 

3 . 37  fur  the  disk  with  =  0,  the  following  results  are  obtained  after 

_  2 

integration  over  x'  from  0  to  1  and  §  from  0  to  2tt.  Letting  A  =  it  r^ 

C  =  2  sin  0  -  2aVl-Q'.  sin  0  cos  [~5  P.  +  (2-P.)  9  ;  4.33 

D  s  y  J  s  L2  J  J  sj 

CT  =  -2  a  J\-a .  sin  G  sin  I  —  P.  +  (2-P.)  0  "1  4.34 

K  J  J  s  L2  j  y  sj 


q  -  ci  =  s  '  °T  - 0 

X  S  ks 


These  equati.  ns  can  be  shown  to  reduce  to  Equations  4.8  through  4.13  ob¬ 
tained  for  the  disk  in  terms  of  the  Maxwell  model  parameters  by  letting 
=  0  and  substituting  the  appropriate  terms  for  a ^ ,  a ^ ,  and  P ^ .  Since 
the  repeated  subscripts  imply  the  summation  over  the  J  beam  components, 
Equation  4.34  for  example  could  be  written 
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J  _ 

C  =  -2  sin  9  S  {<3  .  sin  (6.  +  0  )} 

z  s  j=l  3  J  J  S 

J 

where  .S  CT.  -  1. 

j=l  J 


In  Equations  4,33  and  4.34,  the  first  approximation  to  0^ 


nas  oeen 


utilizied . 


4. 2. 1.2  Drag  of  Cylinder  with  Axis  Perpendicular  to  Flow 

In  the  cylinder  equations  of  chapter  3,  the  angle  0,  given  by 

Equations  3.40,  becomes  if  =  0, 

sin  0  =  -  cos  0  sin  E  4.36 

s 

For  tie  case  when  the  axis  of  the  cylinder  is  perpendicular  to  the  flow 

(0  =  0) ,  sin  0  becomes 

s 

sin  0  =  -  sin  5 

For  §  between  tt  and  3tt/2. 

cos  9  =  -  cos  5 


Therefore,  and  become,  for  §  between  tt  and  3rr/2, 


cos 


G.  =  -a.  </l-ot . 
J  J  j 


f  P.  +  (1-P.K-tt  +  .§)] 
cos  § 


sin  [|  P.  +  (l-P.X-TT  +  |)j 


S.  =  -a  Jl-ci. 
J  J  J 


sin  ^ 


4.37 

4.38 


The  coefficient  of  drag  for  this  special  case  of  the  cylinder  is  found  by 
taking  twice  the  value  obtained  by  integrating  over  5  from  tt  to  3tt/2. 
Letting  A  =  2  r  L,  the  drag  coefficient  of  a  cylinder  is  given  by 

cy 


Pj  +  1,  3 


4,39 
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for  P  =  1,  the  expression  in  the  bracket  should  be  A  value  of  P^  =  3  can 
not  be  meaningfully  applied  to  the  cylinder  problem  unless  a  discontinuous 
linear  relationship  is  chosen  for  the  angle-of-ref lection  law.  Equation  4.39 
is  valid  only  for  values  of  P.  between  0  and  2.  Equation  4.39  can  be  shown 
to  reduce  to  Equation  4.14  for  the  Maxwell  model  by  substituting  the  appro¬ 
priate  values  for  a,,  a.,  and  P..  The  other  coefficients  are  all  zero  for 
J  J  J 

this  special  case  of  a  non-spinning  cylinder  at  9g  =  0. 


4. 2. 1.3.  Drat 


For  a  non-spinning  cone,  Kc  =  0,  the  expression  for  sin  9  reduces  to 


sin  9  =  sin  0  sin  6  -  cos  0  sin  §  cos  6 
s  s 

If  0g  =  ■j,  the  axis  of  the  cone  is  parallel  to  the  flow  and  the  above  equa 
tion  for  sin  9  reduces  to 


sin  0  =  sin  0  sin  6  =  sin  6 
s 


9  =  6 

Therefore,  9  is  found  to  be  independent  of  the  variables  of  integration  £' 
and  §.  For  0g  =  ^  the  angle  (3  is  equal  to  a  constant  of  Since  the  term  x 
in  Equation  3.50  is  also  independent  of  and  |,  Equations  3.50  can  be  inte¬ 
grated  in  closed  form  over  from  0  to  1  and  §  from  0  to  2tt.  Letting 
2 

A  =  tt  rc  the  following  equation  for  the  coefficient  of  drag  of  the  cone  is 
obtained, 

C  =  2  -  2  ajl-a.  cos  P?  P.  +  (2-P.)  d]  4.40 

D  .1  J  12  J  j  J 

The  cone  equation  is  fou>*d  to  have  a  form  similar  to  that  obtained  for  the 

flat  plate  at  angle  of  attack. 
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4. 2. 1.4.  Sphere  Drag 

The  drag  coefficient  of  a  non-spinning  sphere  can  be  obtained  from 
Equation  3.6.1  by  letting  Kg  =  0  and  9^  =  ~.  Then,  sin  0  becomes  . 

sin  9  =  $j' 

the  next  step  is  to  now  change  variables  of  integration  so  that 

sin  y  = 
or  y  =  0 

Substituting  this  into  the  expressions  for  C.  and  and  then  integrating 

TT  _  2 

Equation  3.61  over  5  from  0  to  2rr  and  y  from  0  to  7;  and  letting  A  =  rf  rg  , 
the  following  result  is  obtained 

A  /  1  —  f  p  ^ 

S  ■  2  -  VS  [  p,(4-rj2)  1  ]•  *  °*  4 

where  for  P.  “  0  the  bracketed  term  should  be  zero.  As  in  the  results 
J 

obtained  for  the  cylinder,  Equation  4.41  for  the  sphere  drag  is  valid  only 
for  values  of  P^.  between  0  and  2.  A  value  of  P^  of  4  could  only  occur  if  a 
discontinuous  linear  relationship  were  chosen  for  the  angle-of-ref lection 
law  and  in  that  case  the  equation  obtained  would  not  be  the  same  as  Equa¬ 
tion  4,41. 

4.2.2.  Numerical  Methods 

For  cases  other  than  those  given  in  previous  sections,  the  ex¬ 
pressions  derived  in  chapter  3  for  aerodynamic  properties  contain  integrals 
which  must  be  evaluated  numerically.  In  choosing  a  numerical  technique, 
practical  consideration  must  be  given  to  the  computation  tine  required  to 


t. 


obtain  results  of  suitable  accuracy  which  in  turn  is  dependent  upon  the  com¬ 
puter  being  used  to  perform  the  computations.  The  numerical  evaluation  of 
the  expressions  of  chapter  3  were  made  on  a  Control  Data  160'-!  computer  with 
an  on-line  video  display  system  from  which  graphical  results  were  obtained. 
The  numerical  techniques  employed  are  described  in  the  following. 


4. 2. 2.1.  Single  Variable  of  Integration 

Results  requiring  integration  over  one  variable,  such  as  those  of 
the  cylinder,  were  obtained  by  using  the  Gaussian  integration  formula  over 
and  arbitrary  interval  which  is  given  by 


|  f(y)  dy  =  ~  2  u>.  f(y. ) 

a  ^  1  1 


where 


,b-as  .  b+a 
yi  =  (“)  xi  +  ~T 


The  weights,  u^,  and  abscissans,  x^,  were  obtained  from  Reference  2.2,  The 
twenty  point  formula,  n  =  20,  was  used  in  all  cases.  From  tests,  it  was 
found  that  at  laast  five  place  accuracy  was  obtained  by  using  this  formula, 


4. 2. 2. 2.  Double  Integration 

Results  for  the  disk  and  cone  require  integration  over  two  vari¬ 
ables.  For  these  cases,  the  region  of  integration  was  divided  into  squares 
of  equal  dimensions,  h.  (This  is  possible  for  the  disk  and  cone  since  the 
region  of  integration  is  always  rectangular.)  For  each  square,  a  nine  point 

double  integration  formula  was  used  given  by 
,  9 

— 9  J  ■ J  f(x>y)  <*xdy  =  2  u>.  f(x  ,y  )  +  R 
4li  S  i=l  1  1  1 
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From  Reference  22,  the  abscissas  (x  ,y.)  and  weights  co.  are  given  as 

1  i  i 


-yi} 

u>. 

i 

(0,0) 

16/81 

<±/ 

5  h'±/F> 

25/324 

(0, 

±fT~h) 

10/81 

(±\ 

/f  h,  0) 

10/81 

where  R  =  O(h^) . 

Tests  showed  that  results  of  at  least  5  place  accuracy  could  be 
obtained  in  using  this  formula  by  dividing  the  region  of  integration  into 
about  100  squares. 

The  above  double  integration  formula  was  not  applied  to  the 
sphere  because  the  region  of  integration  defined  by  the  shadow  boundary 
in  the  sphere  is  curved  and  cannot  in  general  be  divided  accurately  into 
squares.  Numerical  resu1 .s  for  the  sphere  were  not  obtained  in  this  study. 

4. 2. 2. 3.  Graphical  Display  of  Results 

In  order  to  evaluate  the  expressions  given  in  Chapter  3  at  various 

values  of  the  quanities  P, ,  9  ,  K  and,  for  the  cone,  6,  a  numerical  inte- 

J  s 

gration  must  be  performed  for  each  case  separately.  Graphical  results  of 
the  variation  in  an  aerodynamic  property  as  a  function  of  these  quantities 
were  obtained  by  evaluating  the  equations  at  six  or  more  separate  values  of 
the  quantity  of  interest  and  then  using  a  six  point  Lagrange  interpolation 
formula  to  plot  the  curve.  For  the  case  of  the  cylinder,  for  example,  288 


84 


separate  numerical  integrations  were  evaluated  and  stored  on  magnetic  tape 
corresponding  to  six  values  for  each  of  the  quantities  P.  and  0g  and  eight 
values  of  K  .  Graphs  of  any  of  the  aerodynamic  properties  as  a  function 
of  any  of  the  three  quantities  can  be  generated  by  obtaining  the  appropriate 
values  from  storage  and  using  the  interpolation  formula  to  obtain  values 
along  the  curve.  Graphs  are  displayed  on  a  X.V.  screen  from  which  pictures 
may  be  obtained.  Numerical  results  are  also  obtained. 

4.3.  Discussion  of  Results 

The  analytical  solutions  obtained  in  this  chapter  illustrate  that 
the  gas-surface  interaction  strongly  influences  the  aerodynamic  properties 
of  both  non-spinning  and  spinning  bodies.  For  non-spinning  bodies,  the 
drag  coefficient  is  the  principal  aerodynamic  property  of  interest.  Fig¬ 
ures  4.1  through  4.8  are  plots  of  some  of  the  analytical  results  obtained 
on  drag  coefficient  as  a  function  of  the  gas  surface  interaction.  ^ 

Figures  4.1  and  4.2  are  plots  of  given  in  Equations  4.8  and 
4.33  for  the  flat  plate  versus  angle  of  attack  for  various  values  of  the 
parameter  or^  of  the  Maxwell  model  and  of  P^  of  the  generalized  model  re¬ 
spectfully.  The  parameters  ot^,  and  a  were  set  equal  to  .5  and  a.  =  1,  and 
the  spin  rate  parameter  is  zero.  These  two  plots  illustrate  the  effect 
of  the  gas-surface  interaction  on  determining  aerodynamic  properties.  The 
effect  of  backscatter  reflections  is  illustrated  in  Figure  4.2  by  the  curves 
for  values  of  P.  between  1  and  2  and  is  seen  to  increase  the  drag  coefficient 
values  at  all  angles  of  attack. 

Figures  4.3  and  4.4  are  plots  of  C  given  in  Equations  4.14  and 


Figure  4.1.  Drag  coefficient  of  a  non-spinning  disk  at  ang 
attack  varying  the  Maxwell  model  parameter  or^. 
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Eicient  of  a  non-spinning  cylinder,  with  axis  of 
perpendicular  to  the  flow,  as  a  function  of  the 
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Figure  4.5.  Drag  coefficient  of  a  non-spinning  cone,  with  axis  of 
the  cone  parallel  to  the  flow,  as  a  function  of  the 
Maxwell  model  parameters. 


4.7. 


.fficient  of  a  aon-spinw"®  sPhete 


function 


Figure 


93 


4.39  for  tha  non-spinning  cylinder  (K  =  0)  with  the  axis  of  the  cylinder 

perpendicular  to  the  flow.  Figure  4.3  gives  the  drag  coefficient,  using 

the  Maxwell  model,  plotted  versus  c^,  for  various  values  of  c^.  Figure  4.4 

is  a  plot  of  C  versus  O',  for  various  values  of  P.  which  are  the  parameters 
D  J  J 

of  the  generalized  model.  For  values  of  P^  between  zero  and  one,  the  range 
in  drag  coefficient  is  the  same  for  both  models.  For  values  of  P^  greater 
than  one,  however,  the  generalized  model  gives  higher  values  of  than  does 
the  Maxwall  model. 

Figures  4.5  and  4.6  are  plots  of  given  in  Equations  4.20  and 

4.40  for  the  non-spinning  cone  (Kc  =  0)  with  the  axis  of  the  cone  parallel 

to  tne  flow.  The  half  angle  of  the  cone  for  these  plots  if  15°.  Figure 
4.5  gives  the  drag  coefficient  versus  for  various  values  of  <*d  which  are 

the  parameters  of  the  Maxwell  model.  Figure  4.6  is  a  plot,  also  of  the 
drag  coefficient,  versus  for  various  values  of  P^  which  are  the  para¬ 
meters  of  the  generalized  model.  In  Figure  4.6,  the  effect  of  backscatter 
and  over-specular  reflections  are  both  illustrated.  The  curve  for  P^  =  -.2 
shows  that  over-specular  reflections  give  values  of  CQ  less  than  2.  In 
general,  over-specular  reflections  decrease  the  drag  coefficient  while  back¬ 
scatter  increases  the  drag  coefficient. 

Figures  4.7  and  4.8  are  plots  of  given  in  Equations  4.30  and 

4.41  for  the  non-spinning  sphere.  Unlik-.  the  previous  plots,  the  two  plots 

of  sphere  drag  coefficients  are  very  similar  for  the  two  cases  of  gas-surface 

interaction  models.  Figure  4.8  illustrates  again  that  backscatter  causes 

increased  values  of  drag  coefficient.  The  results  shown  in  Figure  4.7  us- 

23 

ing  the  Maxwell  model  are  different  than  results  given  by  Cook  which  are 
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similar  to  those  given  in  Figure  4.7  except  that  the  magnitudes  of  drag 

coefficients  are  different  in  the  two  plots  for  the  same  values  of  o;^  and 
0'^,.  For  example,  at  =  1  and  a.j,  =  0,  Reference  23  gives  a  value  of 
for  the  sphere  of  about  2.89.  Figure  4.7  (Equation  4.30)  gives,  however, 
a  value  of  10/3  for  of  the  sphere  at  =  1  and  =  0.  The  difference 
between  the  two  results  is  because  the  results  in  Reference  23  were  ob¬ 
tained  from  Reference  15  in  which  Schamberg  obtained  the  equation  of  for 
the  Maxwell  model  for  3  certain  choice  of  parameters  in  his  model.  In 
reducing  the  Schamberg  model  to  the  Maxwell  case,  the  distribution  of 
velocities  proposed  by  Schamberg  were  still  retained  which,  as  discussed 
in  chapter  2,  introduces  a  factor  of  2/3  in  the  coefficient  of  ^l-o^. 

The  effects  of  spin  on  the  aerodynamic  properties,  is  well 
illustrated,  for  the  Maxwell  model,  by  the  analytical  results  given  in 
section  4.1  of  this  chapter.  In  general,  the  results  show,  for  the  case 
of  the  Maxwell  model 

1.  As  the  spin  rate  increases,  drag  and  (lift)^  increases. 

2 

The  increase  is  proportional  to  K  . 

2.  A  spin  induced  lift  force  in  the  lateral  direction, 

(lift)  , occurs  for  elongated  bodies  at  angles  of  attack. 

The  lateral  force  is  directly  proportional  to  the  spin 
rate  parameter,  K. 

3.  Spin  induced  torques  perpendicular  to  the  spin  axis  are 
experienced  by  all  the  shapes  at  angles  of  attack,  except 
the  cylinder. 
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4.  The  .slowdown,  torque  on  a  spinning  body  is  directly 
proportional  to  K.  Also,  for  the  four  shapes  studied, 
the  slowdown  torque  is  directly  proportional  to 
and  independent  of 

The  effects  of  spin  described  above  were  also  found,  in  general, 
for  the  case  of  the  generalized  model.  The  results  for  the  case  of  the  gen¬ 
eralized  model  were  obtained  by  numerically  evaluating  the  expressions  given 
in  chapter  3. 

The  numerical  results  showed  that  drag  and  (lift)  increased  with 

z 

2 

increased  values  of  K  and  the  increase  was  roughly  proportional  to  K  .  The 
numerical  results  also  showed  that  a  spin- induced  lateral  force  was  experi¬ 


enced  by  the  spinning  cylinder  and  cone  which  is  roughly  proportional  to  K. 
As  for  the, case  of  the  Maxwell  model,  no  lateral  force  was  found  to  occur 
£or  the  spinning  disk.  The  numerical  results  for  torque  perpendicular  to 
the  spin  axis  also  showed  the  same  dependence  on  spin  as  indicated  by  the 
analytical  results  obtained  in  terms  of  the  Maxwell  model  parameters. 

Numerical  results  obtained  for  slowdown  torque  on  a  cylinder  in 
terms  of  generalized  model  parameters  are  given  in  Figure  4.9.  These  re¬ 
sults  show  that  the  slowdown  torque  is  a  function  of  both  and  P  .  For 
the  case  of  the  Maxwell  model  the.  slowdown  torque  is  dependent  on  only 
one  of  the  model  parameters,  The  numerical  results  for  the  slowdown 

torque  of  a  spinning  disk  were  similar  to  those  obtained  for  the  cylinder. 


4.4.  Conclusions 

Two  conclusions  are  evident  from  the  results  obtained  in  this 
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chapter.  First,  the  aerodynamic  properties  of  both  non-spinning  and  spinning 
bodies  is  strongly  affected  by  the  gas  surface  interaction.  Secondly,  dif¬ 
ferent  shapes  have  functionally  different  dependence  on  the  gas-surface 
interaction  parameters  of  a  particular  model. 

These  two  conclusions  suggest  that  gas-surface  interaction  experi¬ 
ments  could  be  performed  in  which  measured  aerodynamic  properties  of  various 
shapes  could  be  utilized  to  obtain  information  on  the  gas-surface  interaction. 
Since  the  spin  induced  torque  and  (lift)^  properties  are  also  strongly  de¬ 
pendent  upon  the  gas-surface  interaction,  considerable  information  on  the 
gas-surface  interaction  could  be  obtained  by  making  measurements  of  these 
aerodynamic  properties  of  spinning  bodies  as  a  function  of  angle  of  attack. 

The  possible  experiments  suggested  by  the  results  of  this  chapter 
would  be  difficult,  if  not  impossible,  using  current  laboratory  methods. 

It  is,  therefore,  proposed  that  tha  aerodynamic  properties  of  spinning 
satellites  be  utilized  to  obtain  information  on  the  gas-surface  interaction. 
The  high  velocity,  largely  neutral,  free  molecular  gas  flow  generated  by  a 
satellite's  motion  through  the  atmosphere  are  almost  ideal  experimental  con¬ 
ditions  for  performing  gas-surface  interaction  experiments.  Tne  remainder 
of  this  study  is  directed  toward  determining  the  feasibility  of  performing 

the  satellite  experiment  suggested  by  the  results  of  this  chapter.  Toward 

% 

this  end,  the  aerodynamic  properties  of  spinning  satellites  are  obtained  in 
the  next  chapter. 
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5.  AERODYNAMIC  PROPERTIES  OF  SATELLITES 

In  chapter  4  the  aerodynamic  properties  of  bodies  in  a  free 
molecular  flow  were  determined  for  the  case  when  the  flow  is  stationary 

with  respect  to  the  body.  In  this  chapter  the  aerodynamic  properties  of 

<y- 

satellites  is  investigated  for  when  the  gas  flow  is  no  longer  fixed  in 
direction . 

The  instantaneous  aerodynamic  properties  of  a  satellite  can  be 
found  at  any  part  in  its  orbit  by  a  suitable  coordinate  transformation 
of  the  equations  obtained  in  the  preceding  chapter.  Since  measurements 
of  satellite  aerodynamic  properties  are  not  usually  made  on  an  instanta¬ 
neous  basis,  the  average  aerodynamic  properties  of  a  satellite  over  one 
orbit  is  investigated  to  determine  how  these  average  properties  depend 
upon  the  gas-surface  interaction. 

The  results  of  this  chapter  will  be  used  in  chapter  6  to  study 
the  feasibility  of  performing  a  satellite  experiment  to  measure  gas-surface 
interaction  parameters.  The  equations  and  procedures  developed  in  this 
chapter  can  also  be  applied  to  the  problem  of  a  tumbling  non-spinning 
satellite.  This  application  is  illustrated  by  obtaining  the  average  drag 
coefficient  for  a  tumbling  non-spinning  disk  as  a  function  of  the  parameters 
of  the  Maxwell  gas-surface  interaction  model. 

5.1.  Coordinate  Transformation 

In  the  equations  developed  in  the  preceding  chapters,  the  aero-' 
dynamic  properties  of  the  various  shapes  were  referred  to  a  coordinate 
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system  which  was  assumed  fixed  in  inertial  space  and  attached  to  the  free 
stream  velocity,  U  .  For  a  satellite  in  orbit  about  the  earth,  the  coor- 

CO 

dinate  system  attached  to  U  rotates  in  inertial  space  as  the  satellite 
travels  in  its  orbit.  Since  the  torques  acting  on  the  satellite  must  be 
referred  to  a  non- rotating  coordinate  system,  a  new  reference  frame  which 
is  attached  to  the  orbital  plane  of  the  satellite  is  chosen.  For  the  pur¬ 
poses  of  illustrating  the  aerodynamic  properties  of  satellites  and  the 

feasibility  of  the  proposed  satellite  experiment,  the  orbit  is  assumed  to 

■ k 

be  circular  and  have  fixed  orientation  in  inertial  space  .  The  atmospheric 
density  is  also  assumed  to  be  constant  over  the  orbit. 

The  initial  orientation  of  the  satellite  spin  vector  in  the  or¬ 
bital  reference  frame  (x0>yo>20)>  *-s  shown  in  Figure  5.1.  The  x0-yo  plane 

is  in  the  orbital  plane  of  the  satellite  and  z  is  the  normal  to  the  or- 

o 

w 4 

bital  plane.  The  free  stream  velocity  vector,  U  ,  rotates  in  the  x  -y 

ca  O  O 

plane  at  a  constant  rate,  a,  equal  to  the  angular  velocity  of  the  satellite's 

—4 

orbit.  For  convenience,  the  satellite  spin  vector,  Q,  is  chosen  to  be 
initially  in  the  xo-zq  plane  at  an  angle  X  from  the  normal  to  the  orbit,  zq. 

With  the  above  definitions,  spherical  trigometry  can  be  used  to 
find  the  angle  9g  in  terms  of  the  two  angles  X  and  a.  This  expression  is 
given  by 

0g  =  sin  ^  (cos  O'  sin  X)  5.1 

where  9  is  the  angle  between  z  and  fi.  (The  notation  and  definition  of  9 
s  s 

k 

In  general,  the  orbital  plane  is  not  an  inertial  frame  of  reference  since 
the  non-spherical  distribution  of  the  earth's  mass  can  cause  the  orbital 
plane  of  a  satellite  to  rotate  in  inertial  space.  This  effect  and  other 
perturbing  effects  of  the  space  environment  are  not  included  in  this  study. 
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used  in  Equation  5.1  is  consistent  with  that  used  in  the  equations  of  pre¬ 
ceding  chapters.) 

5.2.  Instantaneous  Aerodynamic  Properties  of  Satellite 

By  substituting  the  expression  for  0,  given  by  Equation  5.1,  into 
the  equations  of  the  preceding  chapters,  the  aerodynamic  properties  of  a 
satellite  can  be  determined  at  any  position,  a,  in  the  orbit.  For  later 
applications,  it  is  desirable  to  divide  the  instantaneous  force  and  torque 
on  the  satellite  into  components  associated  with  the  x  ,y  ,z  coordinate 


5.2.1.  Components  of  Force 

The  force  acting  on  the  satellite  at  any  position  or  in  the  orbit 

is  divided  into  a  drag  component  in  the  direction  of  U  and  two  components 

of  lift,  normal  to  U  .  For  force,  a  coordinate  system  x  ,v  ,z  is  defined 

oo  ’  J  s  s  s 

which  is  associated  with  the  x0>yo>zc  system  as  shown  in  Figure  5.2.  The 
y^-axis  is  in  the  direction  of  positive  and  is  therefore  in  the  x0_yo 
plane.  The  z-axis  is  in  the  same  direction  as  the  zQ-axis  and  the  x-axis 
then  completes  the  triad.  The  force  on  the  satellite  is  then  divided  into 
components  of  (drag)Vin  the  direction  yj;(lift)x^,  which  is  in  the  plane  of 
the  orbit;  and  (lift)x^,  which  is  normal  to  the  orbital  plane. 

The  components  of  force  defined  in  the  preceding  chapters 
(drag)^,  (lif t)^ ,  and  "(lift^  are  referred  to  the  xf,yf,zf  system  by  per¬ 
forming  the  proper  coordinate  transformations  and  using  Equation  5.1. 


The  results  are 
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cos  X  „  sm  a  sin 

- _  (lift)  - - — 

x  cos  0  z  cos  8 


=  j^(lift) 

4-  [(lift)  +  (lift) 

L  'x  cos  8  v 


s 

cos  X 
z  cos  8 


']  h  +  [drag]  1 

•K 


5.2 


where  are  unit  vectors  in  the  x^,y^,z^  system  and  where 

& 


f  2  2  v 

cos  8g  =  \^l-cos  a  sin  X ) 


Aerodynamic  drag  and  lift  forces  are  known  to  affect  the  elements 

of  a  satellite  orbit.  In  terms  of  the  components  of  force  given  in  Equation 

5.2,  the  principal  effect  of  drag  and  (lift)  (component  of  force  in  the  if 

Xf 

direction  and  in  the  plane  of  the  orbit)  is  to  cause  the  semi-major  axis  of 

the  orbit  to  decay  aud  cause  the  eccentricity  of  the  orbit  to  decrease  to 

zero  (see  for  example  References  24,  25,  and  26).  The  third  component  of 

force,  (lift)  js  normal  to  orbital  plane  in  the  direction  kf.  This  force 
Zf 

causes  a  torque  on  the  orbit  and  has  the  effect  of  causing  the  orbital  plane 

to  precess  in  inertial  space  which  is  analogous  to  a  gyroscope  precessing 

under  the  action  of^an  external  torque. 

A  complete  study  of  the  perturbating  effects  of  drag  and  lift  must 

necessarily  include  also  other  perturbating  forcer  of  the  space  environment 

such  as  solar  radiation  pressure  and  gravity  gradient  forces.  Such  a  study 

is  beyond  the  scope  of  the  present  work.  For  the  purpose  of  the  present 

study,  consider  the  torque  on  a  circular  orbit  caused  by  the  (lift)  com- 

zf 

ponent  of  force.  The  radius  vector,  r  ,  from  the  center  of  the  earth  to 


the  satellite  is  given  by 

r  =  r  sin  a  i  -  r  cos  a  j 


5.3 


The  torque  on  the  orbit  is  then  given  by 
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T  -  r  X  (lift)  kc 
oo  f 


5.4 


i  «••-* 


Since  k..  =  k 
f  o 


T  =  -r  cos  O'  (lift)  i  -  r  sm  ot  (l.ft)  j 
o  o  z,,  o  o  zc  o 

f  r 


5.5 


i  m  * 

f 


Define  the  coefficient  of  torque  on-the-orbit  to  be  given  by 

C_  =  /h  U  2  A  r 
T  O  K  00  o 


5.6 


where  A  is  the  reference  area  of  the  satellite. 
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f  -* 


<  „ » 
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Using  Equations  5.5,  5.6,  and  the  k^  component  of  forces  given 
in  Equation  5.2,  the  coefficients  of  torque  on-the-orbit  about  the  xq  and 


yQ  axis,  respectively,  become 


,r  n  -cos  a  sin  a  sin  X  _ 

K T  ;x  cos  8  I 

o  o  s 


cos  O'  cos  a 

cos  0  L 
s  z 


.2 

er  \  -  ~sin  a  Sln  r 

'  T  Jy  cos  8  °I 

o  o  s 


sin  a  cos  \  ^ 

cos  0  L 
s  z 
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5.2.2.  Components  of  Torque  Acting  on  Satellites 

The  torque  acting  on  the  satellite  at  any  position  a  in  its  orbit 

“4 

is  divided  into  a  slowdown  component  in  the  direction  of  fl,  and  tyo  com- 

— ♦ 

ponents  of  torque  perpendicular  to  Q.  For  satellite  torque,  a  coordinate 


:  .  * 


i  -i 

J 

i"? 

f 

j  Jc 


t  i 
i  ± 


system  x‘,  y  ,  k  is  defined  which  is  associated  with  the  x  ,y  ,z  system 
as  shown  in  Figure  5.3.  The  z^-axis  is  in  the  direction  of  Q  and  is  there¬ 
fore  in  the  x  -y  plane.  The  x  '-axis  is  in  the  same  direction  as  the  y  -axis, 
o  o  s  o 

The  y^-axis  then  completes  the  triad.  The  components  of  torque  on  the 

satellite  are  then  denoted  by  a  slowdown  component,  T^/,  in  the  direction 
-  s 

Q;  a  T  /  component  of  torque  in  the  plane  of  the  orbit  and  in  the  direction 

X 

s 
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x  ;  and  a  T./  component  which  is  out-of-the-piane  of  the  orbit  by  the  angle 
S  T 

JS 

— * 

A ,  "lie  components  of  torque  perpendicular  to  Q  cause  precession  which  also 

has  components  out-of-the  plane  and  in- the- plane  of  the  orbit. 

The  components  of  torque  defined  in  the  preceding  chapters  (T.  , 

s 

T.  ,  and  T  )  are  referred  to  the  x' ,  y' ,  z,  system  by  performing  the 
Js  s  s  s  s 

proper  coordinate  transformation  and  using  Equation  5.1.  The  results  are 

cos  a  cos  X 
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-tv 


T. 
..  1 


cos  6 
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sin  a 


-T 


sin  a  ~j  . 
0sJ  1 
COS  Q !  COS  V 


j  cos  e  j  s 

s  s 


A  -r  T.  . 

cos  6  j  cos  9 
s  s  s  s 
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!  T,  1  k  ' 
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5.9 


This  equation  and  the  equation  of  the  previous  section  allow  the  resui_s>  of 
chapter  4  to  be  used  to  find  the  instantaneous  aerodynamic  properties  of  a 
satellite  at  any  point  a  in  its  orbit. 


5.3.  Average  Aerodynamic  Properties  of  Satellites 

The  average  aerodynamic  properties  of  a  body  is  defined  by  the 


equation 


-  1  T' 

C  =  C  dt 


o 

where  C  represents  the  aerodynamic  property  being  averaged  and 
time  interval  over  which  the  average  is  taken.  As  a  satellite 
through  its  orbit,  the  velocity  vector  rotates  in  the  x0_yQ 
The  position  of  the  velocity  vector  is  determined  by  the  angle 
given  time  t.  For  he  purposes  of  the  present  study,  consider 


5.10 

T7  is  the 
travels 
plane . 
a  at  a 
the  case  of 
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a  satellite  in  near-earth  circular  orbit.  If  the  altitude  loss  due  to  the 
action  of  aerodynamic  drag  is  assumed  to  be  small,  the  angular  velocity  of 
the  satellite’s  orbit  can  be  approximated  by  a  constant  determined  by  the 
average  radius  of  the  orbit  for  a  single  pass  about  the  earth.  Similarly, 
an  average  circular  velocity  equal  to  the  free  stream  velocity  U  and  an 

CO 

average  atmospheric  density  p  can  be  approximated  by  constants. 

For  a  constant  angular  velocity  a  =  cu,  Equation  5.4  can  then  be 
changed  to  an  integral  of  the  angle  a  over  2rr  for  one  complete  orbit.  The 
equation  for  the  average  aerodynamic  properties  of  a  satellite  over  one 
complete  orbit  is  then  given  by 

C  =  “  r2lT  C  da  5.11 

2tt  jo 

Equations  5.2,  5.7,  5.8  and  5.S  can  then  be  substituted  into  Equation  5.11, 
with  the  appropriate  equations  from  the  preceding  chapter,  to  determine  the 
average  properties  of  satellites. 

5.3.1.  Spinning  Disk  , 

In  general,  the  evaluation  of  Equation  5.11  involves  the  use  of 

numerical  techniques.  For  the  case  of  the  properties  of  the  spinning  disk 

in  terns  of  the  Maxwell  model  parameters,  however,  Equation  5.11  can  be 

evaluated  analytically.  The  procedure  is  to  first  substitute  Equations  4.8 

through  4.13  for  the  disk  into  the  equations  fot  the  instantaneous  aero  • 

dyar.mic  properties  given  by  Equations  5.2,  5.7,  5.8  and  5.9.  Equation  5.1 

is  used  to  replace  9  in  terms  of  cv  and  X.  These  results  are  then  sub- 

s 

slituled  into  Equation  5.11  and  integrated  to  find  the  average  properties 
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over  one  orbit.  Two  cases  are  considered:  1)  spinning  disk  with  one  side 
exposed  to  the  flow  and  2)  spinning  disk  with  both  sides  exposed  to  the 
flow.  The  results  of  the  first  case  can  be  applied  to  problems  such  as  a 
spinning  cylinder  or  spinning  cone  in  which  the  ends  of  the  cylinder  and 
base  of  the  cone  can  be  represented  e~  spinning  disks. 


5. 3. 1.1.  Spinning  Disk  with  One  Side  Exposed  to  the  Flow 

Following  the  procedure  outlined  above,  the  results  for  the 
average  properter  of  a  spinning  disk  with  one  si.de  exposed  to  flow  are 


7  _  8  .  3  ,  ,2  .  8  .  3  ^ 

Cn  =  —  sin  A  +  x,  {-  sm  X  -  —  sin  X) 
D  3tt  d  tt  3tt 


+  otf/l-aT  [A  sin2  X  +  sin'"'  X  +  ~  sin“  \)]  5.12 


2  ,1  .2 


3  .  4 


CT .  *  =  ° 
1 
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_  1 


i .  i 

h 


.  ajl-am  K,  sin  X  cos  X 
8  dv  T  d 


5.13 


5.14 


C  -  -  2  aA  K,  sin  X 
1,  /  d  d 

k 

s 


vCT  ^ x 
o  o 


8  2 
—  (l-or, )  sin  X  cos  X 
3tt  d 


5.15 


2  ,1 


L  3  3 

-  a,  \/X-a,„  K Cfr  sin  X  cos  X  +  —  sin  X  cos  X)  5.16 
d  T  d  16  64 


(CT  }y  =  0 
o  ^o 


5.17 


5 .3. 1.2.  Spinning  Disk  with  Both  Sides  Exposed  to  the  Flow 

The  procedure  to  be  followed  for  the  case  of  the  disk  with  both 
sides  exposed  to  the  flow  is  the  same  as  above  except  that,  when  the  flow 
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is  on  the  "back-side"  of  the  disk,  appropriate  changes  in  signs  must  be  made 
on  Equations  4.8  through  4.13  to  account  for  the  fact  that  the  spin  rate 
vector  0  is  in  the  opposite  direction.  The  results  are 


P  _  16  .  3  .  .4  .  .  lb  .3,. 

Cn  ^  —  sm  X  +  o' ,  (~  sm  \  —  sm  a) 
o  on  d  rr  3tt 


+  adv//1“aT  Csin2  X.  +  Kd2  (-g  sin2  X  +  sin4  X)]  5.18 


4  ctrd  Kd  sin  X 


5.19 


5.20 

5.21 


(CT  )x  =  -  (1-c^)  sin2  cos  X 

o  o 

-  Kd2  (-|  sin  X  cos  X  +  sin^  X  cos  X)  5.22 

(C,  )  -  0  5.23 

o  yo 

For  the  case  of  both  sides  of  the  disk  exposed,  the  torque  perpen¬ 
dicular  to  the  spin  axis  of  the  satellite  is  found  to  average  to  zero  over 
one  orbit.  The  average  drag  and  torque  on-the-orbit  about  the  x^  axis  is 
found  to  double  as  would  be  expected. 


5.3,2.  Spinning  Spherical.  Satellite 

The  average  coefficients  of  torque  on  a  spinning  spherical  satellite 
can  be  obtained  in  term.?  of  the  parameters  of  the  Maxwell  model  from  Equations 
4.27  through  4.29.  The  results  are 


no 


=0  5.24 

i.  . 

1  / 
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C  =  7  K  a ,  sit:  X  ccs  X  5.25 

1  .  4  s  a 
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CT  =  -  \  Kg  ad  (6  -  sin2  X)  5.26 
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Unlike  for  the  case  of  the  disk,  the  average,  over  one  orbit,  of  torque 

perpendicular  to  the  spin  axis,  C  ,  is  found  to  be  finite  for  the  case 

1  .  / 

Js 

of  a  spherical  satellite. 

5.3.3.  Spinning  Cone-Disk  Composite  Satellite 

Consider  a  satellite  composed  of  a  cone  spinning  about  the  axis 
of  the  cone  and  a  flat  base  represent'd  by  a  spinning  disk.  For  convenience, 
consider  the  center  of  mass  of  the  satellite  to  be  at  the  center  of  the  base 
of  the  cone.  For  purposes  of  illustration  the  gas-surface  interaction  will 
be  assumed  to  be  of  the  Maxwell  type.  The  equations  expressing  the  aero¬ 
dynamics  of  the  spinning  cone  at  angles  of  attack  have  been  obtained  in 
analytical  form  in  terms  of  the  Maxwell  model  parameters  and  are  given  by 
Equations  4.20  through  4.26  in  chapter  4.  Since  a  cone-disk  composite 
satellite  is  a  convex  body,  Equations  5.12  through  5.17  fox-- the -disk  with 
one  side  exposed  tp  the  flow  can  be  applied  directly.  The  si"n  of  C^, 

is 

(Equation  5.14)  is' reversed  for  this  application  since  only  the  "back¬ 
side"  of  the  spinning  disk  is  exposed  to  the  flow. 

The  average  properties  of  the  cone  part  of  the  spinning  cone-disk 
composite  satellite  are  obtained  using  the  numerical  techniques  described 
previously  to  perform  the  integration  over  o'  in  Equation  5.11.  For  example, 
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for  a  cone  of  a  half  angle  of  15°  and  X  =  45°  the  following  results  were 
obtained  numerically. 

CD  =  ~  [15.4  4  .45  ad  4  o^/l-o^  (10.5  4  2.0A  K2)]  5.27 

CT  ;  ~  [-3.12  4  3.12  ad  -  otjfl  -«T  <1.3  4  8.7  X  ID-3  K2)]  5.28 
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-  ctgj l-QfT  (2.3 
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X  10  K  4  8.2  K  )] 


5.29 


5.30 


(ct  K  =  h t-1-22  +  1,22  ad '  (-85  “  -38  f-2)]  5-31 

o  o 

(CT  )  =  ^  [2.76  K  otd  -  (1.07  K)]  5.32 

o  yo  c 

The  average  properties  of  a  cone-disk  satellite  are  a  function  of 
the  cone  half  angle  6  and  the  angle  This  functional  dependence  is  illus 
trated  in  Figure  5.4  where  the  average  drag  coefficient  is  shewn  plotted  as 
a  function  of  the  angle  X  for  various  values  of  the  angle  6. 


5.4.  Average  Aerodynamic  Properties  of  Tumbling  Bodies 

The  average  aerodynamic  properties  of  a  body  which  is  tumbling  in 
a  random  manner  can  be  found  using  the  techniques  developed  in  this  chapter. 
Let  the  orientation  of  the  body  with  respect  to  the  flow  be  determined  by 
the  angles  a  and  X.  If  the  orientation  is  completely  random,  all  values  of 
X  between  0  and  tt  and  of  a  between  0  and  2tt  occur  for  an  equal  amount  of 
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time.  The  average  of  an  aerodynamic  property  for  a  random- tumbling  body 
is  expressed  by 

2rr 
0 

where  C  is  the  aerodynamic  property  being  averaged  (see  Reference  27), 

A  single  example  of  the  application  of  Equation  5.33  is  given  in 
the  following  for  the  case  of  a  non-spinning  disk. 


VT 

J  C  sin  X  dX  da  5.33 

0 


5.4.1.  Average  Drag  Coefficient  of  a  Random-Tumbling  Disk 

The  instantaneous  drag  coefficient  of  a  disk  in  terms  of  the 
parameters  of  the  Maxwell  model  is  given  by  Equation  4.8.  In  terms  of 
the  angles  a  and  >.  the  instantaneous  drag  coefficient  for  zero  spin  rate 
is  given  by 

.  3 


—  3  3  3  3 

(CL)  =  4  cos  a  sin  X  +  a,(2  cos  a  sin  X  -  4  cos  a  sin  X) 
D  . .  .  d 


disk 


+  Q'^yi-a'^,  2  cos  a  sin  X 


5.34 


Due  to  symmetry,  the  integration  of  Equation  5.33  for  the  disk  can  be  ob¬ 
tained  as  eight  times  the  integration  over  a  from  0  to  ^  and  X  from  0  to  ~ . 
The  result  is 

<CD)  =  1  +  |  5.35 

disk 

The  average  drag  coefficient  of  the  random- tumbling  disk  is  found  to  be 
considerably  less  than  the  maximum  possible  drag  coefficient  of  four  for  a 
stationary  disk  normal  to  the  flow.  The  average  drag  coefficient  gltfen  in 
Equation  5.35  is  exactly  half  the  drag  coefficient  of  a  sphere  (see  Equation 
4.30). 

Approximation  techniques  have  been  commonly  employed  to  determine 
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average  drag  coefficients  of  tumbling  bodies.  The  results  of  two  approxi¬ 
mation  techniques  applied  to  the  disk  problem  are  given  for  comparison  to 
the  results  given  in  Equation  5.35.  Using  results  obtained  by  Schamberg 
(Reference  15) 

(CD)  «  3  +  2  5.36 

Using  techniques  developed  by  Sentman  and  Neice  (Reference  27) 

(CQ)  «  .85  -  .42  <*d  +  .42  ad<yi-a~  5.37 

The  approximation  techniques  used  in  Reference  15  (Equation  5.36)  are  seen 
to  give  a  drag  coef f icient  three  times  greater  than  that  given  by  Equation 
5.35.  The  approximation  technique  of  Reference  27,  on  the  other  hand,  gives 
slightly  lower  results  than  that  given  by  Equation  5.35,  but  the  percent 
error  is  considerably  less  than  the  results  of  Reference  15. 

5.5.  Discussion  of  Results 

The  results  obtained  in  this  chapter  have  been  limited  to  the  use 
of  the  Maxwell  model  because  analytical  expressions  could  be  obtained  in 
most  of  the  cases  studied.  The  procedures  used  to  obtain  these  results  are, 

v 

however,  equally  applicable  to  problems  using  the  generalized  model  if 
numerical  methods  are  employed. 

The  results  obtained  have  shown  than  the  average  aerodynamic  pro¬ 
perties  of  a  satellite  are  of  the  same  form  as  those  obtained  in  chapter  4 
for  stationary  flow  conditions.  The  average  aerodynamic  properties  are 
found  also  to  be  strongly  dependent  upon  thf  parameters  of  the  gas  surface 
interaction  model.  In  general,  the  average  aerodynamic  properties  of 
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the  following  form  as  a  function  of  the  para- 

“d  A2  +  “d^T  A3 

®d  B2  +  V7^  B3 

_  5.38 

*d  °2  +  V1-0^  C3 

where  the  coefficients  A,  B,  etc.,  are  functions  of  the  angle  X  and  the  spin 
rate  parameter  K,  and  also  the  half  angle  6  for  the  cone. 

It  is,  therefore,  evident  that, as  for  bodies  in  a  stationary  flow, 
the  possibility  exists  of  utilizing  measured,  average  satellite-aerodynamic 
properties  to  obtain  specific  information  on  the  gas-surface  interaction. 

The  feasibility  is  assessed  in  the  next  chapter  of  performing  a  satellite 
experiment  to  determine  gas-surface  parameters  by  measuring  the  average  aero 
dynamic  properties  of  satellites. 


satellites  can  be  expressed  in 
meters  of  the  Maxwell  model. 
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6.  FEASIBILITY  OF  A  SATELLITE  EXPERIMENT  TO  DETERMINE  THE  GAS-SURFACE 

INTERACTION  PARAMETERS  AND  THE  ATMOSPHERIC  DENSITY 

The  results  of  the  preceding  chapter  have  indicated  clearly  the 
possibility  of  utilizing  measured  aerodynamic  properties  of  satellites  to 
obtain  information  on  the  gas-surface  interaction.  The  feasibility  of 
this  experimental  concept  is  discussed  in  this  chapter. 

The  obvious  advantage  of  designing  a  gas-surface  interaction  ex¬ 
periment  using  satellites  is  that  a  free  molecular  flow  of  high  velocity 
neutral  molecules  is  generated  by  the  motion  of  the  satellite  in  its  orbit. 
As  discussed  earlier,  it  has  not  been  possible  to  duplicate  these  flow 
conditions  in  the  laboratory.  As  with  all  experiments,  even  those  per¬ 
formed  under  supposedly  controlled  conditions,  information  obtained  on  the 
gas-surface  interaction  from  the  analysis  of  satellite  aerodynamic  pro¬ 
perties  would  be  subject  to  certain  uncertainties  associated  with  the  ex¬ 
periment.  For  satellite  experiments  of  the  type  proposed  in  this  study, 
uncertainties  associated  with  the  space  environment  could  influence  the 
interpretation  and  accuracy  of  results.  A  complete  discussion  of  the 
major  uncertainties  and  their  influence  on  the  proposed  experiment  is 
given  ir.  the  latter  part  of  this  chapter.  However,  since  the  uncertainty 
associated  with  the  orbital  gas-density  is  of  major  importance  to  the 
proposed  experiment,  this  subject  is  discussed  in  the  following. 

The  reason  for  the  uncertainty  in  orbital  gas-density  can  be 
traced  to  the  lack  of  knowledge  or.  the  drag  coefficient  of  satellites 
which,  in  turn,  is  related  to  the  uncertainty  in  the  gas-surface  inter¬ 
action.  As  an  example,  consider  the  measurement  of  the  drag  of  a  satellite 
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in  which  the  drag  is  related  to  the  density  by  the  drag  equation  given  by 

Drag  =  \  p  vj'  A  CD 

In  order  to  determine  the  atmospheric  density  p  from  the  drag  measurement, 

the  drag  coefficient,  C^,  of  the  body  must  be  known.  Due  to  the  lack  of 

knowledge  on  the  gas-surface  interaction,  an  assumed  value  of  the  drag  co- 

28 

efficient  of  2.2  is  commonly  used  to  reduce  drag  data.  From  the  results 
of  chapters  4  and  5,  it  is  seen  that  this  assumption  could  lead  to  consider¬ 
able  error,  perhaps  50%  or  greater,  in  the  determination  of  density  depend¬ 
ing  upon  the  actual  values  of  the  gas-surface  interaction  parameters  and  the 
shape  of  the  body. 

It  is  apparent  from  the  above  discussion  that  a  satellite  experi¬ 
ment  designed  to  obtain  information  on  the  gas-surface  interaction  must, 
also,  simultaneously,  obtain  information  on  the  atmospheric  density.  For 
this  reason,  then,  the  satellite  experiments  proposed  in  this  study  ere  de¬ 
signed  to  determine  the  value  of  the  atmospheric  density  in  addition  to 
determining  the  values  of  gas-surface  interaction  parameters. 

6.1,  Gas-Surface  Interaction  Experiments  Utilizing  Satellite  Aerodynamic 
Properties 

Two  criteria  must  be  met  in  order  to  utilize  maasured  aero¬ 
dynamic  properties  of  satellites  to  obtain  information  on  the  gas-surface 
interaction  and  atmospheric  density.  The  first  criteria  is  that  the  number 
of  aerodynamic  properties  measured  must  at  least  equal  the  number  of  un¬ 
known  quantities  to  be  determined.  The  second  criteria  is  that  the  equations 
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expressing  the  measured  aerodynamic  properties  in  terms  of  the  unknowns 
must  form  an  independent  set  of  equations  which  can  then  be  solved  for  the 
unknowns . 

Two  general  classifications  can  be  made  of  possible  methods  of 
utilizing  satellite  data  which  could  satisfy  the  basic  criteria  given  above. 
These  classifications  of  possible  methods  are: 

Method  1:  Utilizing  data  (such  as  drag)  frorr.  a  number  of 
differently  shaped  and/or  differently  oriented 
satellites . 

Method  2:  Utilizing  data  on  a  number  of  aerodynamic  pro¬ 
perties  of  a  properly  designed  satellite. 

The  methods  of  analysis  are  similar  in  both  classifications.  In  Method  1, 
data  on  past  or  existing  satellites  would  be  used  while  in  Method  2  a 
satellite  is  to  be  designed  for  the  specific  purpose  of  obtaining  informa¬ 
tion  on  the  gas-surface  interaction.  It  is  not  suggested,  however,  that 
either  of  the  two  methods  be  used  exclusively.  In  fact,  as  will  be  pointed 
out  later,  a  combination  of  the  two  methods  appears  desirable  in  terms  of 
a  long  range  program  to  obtain  information  on  the  gas— surface  interaction 
as  a  function  of  the  surface  conditions  and  also  to  obtain  information  on 
the  density  and  composition  of  the  orbital  gas  environment. 

The  procedure,  advantages,  and  disadvantages  associated  with 
the  two  methods  are  discussed  in  the  following  two  sections.  Later  in  this 
chapter  the  feasibility  will  be  illustrated  of  using  Method  2  to  determine^ 
the  parameters  of  the  Maxwell  model  and  the  atmospheric  density. 
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6.1.1.  Method  1:  Utilization  of  Data  From  Satellites  of  Different  Shapes 

In  chapter  5  it  was  found  that  the  average  aerodynamic  properties 
of  satellites  are  dependent  upon  the  parameters  of  the  gas-surface  inter¬ 
action.  In  addition,  the  angle  X, which  is  the  angle  giving  orientation  of 
the  satellite  spin  axis  with  respect  to  orbit, was  also  found  to  have  an 
influence  on  the  coefficients  of  the  unknown  parameters.  These  two  factors 
are  to  be  considered  in  utilizing  data  from  satellites  of  different  shapes 
and  orientation  in  making  the  measurements. 

To  begin  the  discussion  of  Method  1,  consider  that,  for  example, 
drag  data  were  available  on  a  spinning  disk,  a  spinning  cone  (with  flat 
base)  and  a  non-spinning  sphere.  Consider  that  each  of  these  satellite 
shapes  has  the  same  surface  properties,  are  all  in  the  same  circular  orbit, 
and  that  the  gas  density  is  constant  over  the  orbit.  Under  these  considera¬ 
tions  and  the  assumption  that  the  gas-surface  interaction  is  of  the  Maxwell 
type,  the  drag  measurements  of  the  three  would  provide  sufficient  informa¬ 
tion  for  the  determination  of  the  parameters  o^,  and  p.  The  fact  that 
the  values  of  c^,  and  p  can  be  determined  can  be  verified  from  the  equa¬ 
tions  given  in  chapter  5  for  the  spinning  disk  with  both  sides  exposed 
(Equation  5.18)  and  the  numerical  results  given  for  the  cone-disk  satellite 
(Equation  5.27).  The  average  drag  for  a  non-spinning  sphere  would  be  the 
same  as  that  given  in  chapter  4  for  stationary  flow,  (Equation  4.30).  These 
three  equations  satisfy  both  criteria  of  the  experiment,  as  can  be  easily 
verified.  Using  the  same  set  of  equations  one  could  also  verify  that  two 
spinning  disks  at  different  angles  X  to  the  orbital  plane  could  also  provide 
a  system  of  independent  equations  when  combined  with  either  the  cone-disk 
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or  sphere  equations.  One  can  also  easily  verify  that  measurements  of  drag 
on,  for  example,  three  different  sized  spheres  would  not  lead  to  an  indepen¬ 
dent  set  of  equations  which  could  be  solved  for  or^,  o^,  and  p. 

The  possible  experiments  discussed  above  concern  the  determination 
of  the  parameters  of  the  Maxwell  model.  The  same  procedure  would  be  used  if 
instead  the  generalized  gas-surface  interaction  model  is  used.  In  order  to 
use  the  generalized  model,  numerical  methods  would  have  to  be  employed 
throughout  the  analysis. 

The  ideal  conditions  considered  above  would  not  occur  in  practice. 
In  the  actual  application  of  Method  1  of  using  past  drag  data,  assumptions, 
which  would  introduce  errors  into  the  analysis,  must  be  made  in  order  to  use 
data  from  different  satellites.  One  assumption  that  may  have  to  be  made  is 
that  the  gas-surface  interaction  is  the  same  on  all  the  satellites  being  used 
in  the  analysis.  Also,  in  order  to  analyze  dat*-  from  satellites  in  non¬ 
circular  orbits,  an  assumed  atmospheric  density  model  must  be  employed.  The 
same  model  would  also  have  to  be  employed  in  order  to  use  any  data  from 
satellites  which  are  in  different  orbits.  In  addition,  the  satellites  would 
have  to  be  assumed  to  be  approximated  by  convex  shapes,  since  there  is  at 
present  no  adequate  method  available  to  analyze  concave  shapes. 

In  light  of  the  many  assumptions  which  must  be  made  in  the  analysis 
of  existing  satellite  data,  it  is  concluded  that  data  from  a  large  number  of 
satellites  would  have  to  be  analyzed  and  correlated  in  order  to  reduce  the 
errors  introduced  by  the  assumption.  The  results  of  such  an  undertaking 
could,  however,  yield  a  considerable  amount  of  information  not  only  on  the 
gas-surface  interaction,  but  also  on  the  validity  of  the  assumptions  used,. 
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such  as  the  atmospheric  density  model. 


6.1.2.  Method  2:  Utilization  of  Data  From  a  Single  Satellite 

The  design,  manufacture  and  launching  of  a  satellite  is  an  expen¬ 
sive  operation  and  therein  is  the  principal  disadvantage  to  Method  2  in  com¬ 
parison  to  the  less  expensive  analysis  of  existing  data  of  Method  1.  There 
are,  however,  many  advantages  to  a  Method  2  analysis  which  would  utilize  data 
from  a  satellite  which  is  specifically  designed  to  obtain  information  on  the 
gas-surface  interaction.  The  principle  advantage  would  be  the  accuracy  of 
results  obtained,  an  accuracy  which  in  all  probability,  could  not  be  obtained 
in  an  analysis  of  past  satellites. 

In  order  to  illustrate  the  basis  of  Method  2,  consider  the  cone- 
disk  satellite  which  was  analyzed  in  chapter  5.  The  average  drag  and  torque 
coefficients  of  the  cone-disk  satellite  were  of  the  following  form  in  terms 
of  the  parameters  of  the  Maxwell  model. 


C_.  =  A,  +  A  .a,  +  A  J  1-0' 

D  1  /.  u  J  a  i 


CT  B1  +  B2ad  +  B3ad  ^1_QfT 
P 

CT  =  adC 
s 

where  is  the  average  drag  coefficient,  C^,  is  the  average  coefficient  of 

P  _ 

torque  which  acts  perpendicular  to  the  satellite  spin  axis  (C^,  will  be 

P 

interpreted  later  as  being  the  coefficient  of  either  the  i1  -  or  j 1  component 

s  s 

of  torque),  and  C„  is  the  average  coefficient  of  slow  down  torque  which  is 
s 

in  the  direction  of  the  spin  axis.  The  coefficients  A^,  A^,  l  B^,  B2,  B^, 
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and  C  are  functions  of  X,  6,  and  K  as  discussed  in  chapter  "> . 

Consider  now,  that  D,  T  ,  and  T  are  measurable  quantities  so  that 

P  s 

the  following  set  of  equations  could  be  obtained 

D*  =  D/~  U2  A  =  /  .  p  4  A.  p  a ,  +  A,  p  o,  -o'  6.4 

2®  '  2rd  i  r  a  T 

T  *  =  T  /-  n2  A  r  K  =  B,  o  +  B  p  a  +  B,  p  a.  Ji 6.5 
p  p^to  i.  t.  a  j  a  i 

T  '  =  T  /~  U2  A  r  K  =  C  p  O' ,  6.6 

s  s  2  co  v  d 

k  k  k 

in  which  D  ,  T  ,  and  T  are  known  from  measurements.  These  three  equations 
P  s 

satisfy  the  two  basic  criteria  established  previously.  Solutions  for  p,  Of^, 

•k 

and  q^j,  can  be  obtained  analytically  in  terms  of  the  measured  quantities  D  , 

*  * 

T  ,  and  T  and  the  coefficients  A, ,  A„,  etc.  These  results  are 
p  s  L  l 

T  *  CA.  -  D*CB_  +  T  *  (A  B  -  A,.B„) 

P° - - -  6.7 


where  E  =  A^B^  -  A^B^. 


* 

+T  E 

_ s _ 

ra3  -  D*CB3  +  <A2B3  '  A3E2)  Ts* 

C&1  +  D*CBj  -  Ts*  (Ajilj  -  AjP.j) 
- 

T  E 
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The  cone-disk  satellite  is  one  example  in  which  the  aerodynamic 
properties  do  provide  equations  which  satisfy  the  basic  criteria.  This  is 
not  generally  true  of  other  satellite  shapes.  For  example,  tht  sphere  equa¬ 
tions  given  in  chapter  5  show  that  C^,  and  for  the  sphere  are  not  inde- 

P  s 

pendent.  In  another  example,  the  equations  for  the  disk  with  both  sides 

exposed  show  that  CT  is  zero  and  therefore  the  first  criteria  is  not 

P 
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satisfied . 

As  in  Method  1,  the  Maxwell  model  is  used  as  a  convenient  example, 
The  technique  of  analysis  is,  however,  also  applicable  to  the  determination 
of  the  parameters  of  the  generalized  model,  if  numerical  methods  ar  -•  employed. 


6.2.  Accuracy  of  a  Satellite  I’xneriment 

The  accuracy  of  determining  the  gas -surface  interaction  parameters 
and  the  atmospheric  density  in  a  satellite  experiment  is,  of  course,  depen¬ 
dent  upon  the  accuracy  of  the  measurements  of  the  average  satellite  aero¬ 
dynamic  properties.  The  accuracy  is  also  dependent,  however,  on  the  shape  or 
shapes  of  satellites  from  which  data  is  obtained.  For  example,  if  d-ag 
measurements  were  made  on  three  satellites  which  differed  only  slightly  from 
that  cf  a  sphere,  a  solution  would  not  be  possible  because  even  minor  errors 
in  drag  measurements  would  be  greatly  amplified.  These  factors  are  illus¬ 
trated  by  taking  the  specific  example  of  the  cone-disk  satellite  used  in  the 
previous  section. 

The  accuracy  of  determining  the  unknown  parameters  p,  o^,  and  ci^ 

* 

using  a  cone-disk  satellite  is  dependent  upon  the  accuracy  of  measur-'  ig  D  , 

*  * 

T  and  T  .  An  estimate  of  the  errors  in  the  determination  of  p,  or, ,  and 
p  s  d 

Gf^  can  be  made  by  taking  the  first  partial  derivatives  of  each  of  the  un- 

*  *  * 

known  parameters  with  respect  to  D  ,  T  ,  and  Tg  .  Taking  the  derivatives 
of  Equation  6.7,  6.8,  and  6.9,  the  following  estimation  of  the  errors  is 
obtained , 


Aa  =  _  c  13  AD. 


D  E 


A  AT  *  (A„B„  -  AJ&_) 

+  r  -2  — 2—  ■  a  - 3_2_ 

S  E  *  +  ad  E 
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AT  * 
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ji'T^  /T  I  and  |ATg  /Tg  |  are  all  less  than  or  equal  to  some  maximum  possible 
valu..  of  |  AM/M |  .  Then,  for  the  worse  case,  when  all  the  errors  are  added,  the 
maximum  possible  errors  in  the  determination  of  p,  a^,  and  would  be  given 


by 
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6.15 


The  quantities  multiplying  |~~|  on  the  right  hand  side  of  the 
above  equations  represent,  for  a  given  | AM/M ( ,  the  maximum  error  which  is 
associated  with  the  shape  and  orientation  of  the  satellite.  That  is,  these 
quantities  vary  as  a  function  of  the  cone  half  angle  6  and  the  orientation 
of  the  satellite  spin  axis  with  respect  to  the  orbit  which  is  determined  by 
the  angle  X.  Using  numerical  methods,  the  term  multiplying  |AM/M|  in  Equa¬ 
tion  6.15  was  evaluated  and  the  results  are  shown  plotted  in  Figure  6.1  as 
a  function  of  X  for  various  values  of  6. 

Figure  6.1  shows  that  some  X-6  combinations  are  definitely  better 
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than  others.  In  fact,  for  some  combinations  the  coefficient  of  } AM/M |  be¬ 
comes  large,  implying  that  an  experiment  employing  those  \-6  combinations 
would  not  yield  acceptable  results.  The  X-b  combinations  which  give  mini¬ 
mum  possible  error  are,  of  course,  the  most  desirable.  Results  similar  to 
those  shown  in  Figure  6.1  were  also  obtained  for  the  error  estimates  of  p 

and  of,.  The  error  in  determining  a  is  directly  related  to  the  value  of  Of  , 
d  T  g* 

as  seen  in  Equation  6.15.  As  ^  approaches  zero,  the  error  in  the  determina 

tion  of  would  become  infinitely  large.  This  is  expected  since  for  the 

Maxwell  model  the  value  of  a,  de  ermines  the  fraction  of  reflected  force 

d 

which  is  dependent  upon  the  parameter  o' .  Thus,  as  of^  becomes  small,  the 
effect  of  on  the  force  and  torque  on  the  satellite  is  greatly  dii>{nished. 
At  cTj  =  0,  the  parameter  becomes  meaningless. 

The  results  shown  in  Figure  6.1  were  obtained  using  the  coeffi¬ 
cients  of  the  j'  component  of  torque.  Results  were  also  obtained  using 
s 

the  i^  component  of  torque.  Results  were  obtained  for  a  wide  range  of  6 
values  for  both  cases.  These  results  showed  the  following  in  general 

1.  Utilization  of  measured  values  of  the  j/  component  of 

s 

torque  give  consistently  lower  values  of  maximum  possible 
error  in  the  determination  of  p,  Qf^,  and  ar^,  as  compared 
to  the  results  obtained  utilizing  measurements  of  the  i^ 
component  of  torque. 

2.  Values  of  X  greater  than  45°  give  lower  error  when  the 
j  'c  component  of  torque  is  utilized  while  values  of  X 
less  than  45°  give  lower  error  when  the  component  of 


torque  is  utilized. 
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3.  The  error  in  determining  is  significantly  different  for 

using  either  the  i'-or  j  '  -component  of  torque.  Use  of  the 
s  s 

i ^component  of  torque  was  found  to  consistently  give  maxi- 
s 

mum  error  values  of  at  least  an  order  of  magnitude  greater 

than  that  obtained  using  the  -component  of  torque. 

The  difference  in  possible  accuracy  between  using  either  the  i'- 

s 

component  or  the  -component  of  torque  is  best  explained  by  referring  to 
the  numerical  results  given  in  Equations  5.27,  5.28,  and  5.29.  These  equa¬ 
tions  show  that  the  drag  and  the  i^-component  of  torque  are  functionally 
similar  in  their  dependence  on  the  unknown  p,  a^,  and  a^.  The  j ^-component 
of  torque  has,  however,  a  functionally  different  dependence  on  the  unknowns 
than  either  the  drag  or  slow  down  torque.  Therefore,  even  though  it  is 

possible  to  use  either  the  i  -  or  j'  -  components  of  torque  in  an  experiment, 

s  s 

use  of  the  j^-component  of  torque  is  more  desirable  in  terms  of  the  accuracy 
of  the  experiment. 

The  procedure  just  outlined  can  be  used  to  find  an  optimum  satel¬ 
lite  design  for  performing  a  gas-surface  interaction  and  atmospheric  density 
experiment.  In  such  an  optimization  study,  factors  such  as  size  and  weight 
requirements,  orbital  regression  effects,  and  others  would  be  considered  in 
addition  to  the  satellite  shape.  Considering  only  6  and  however,  the 
cone-disk  satellite  was  found  to  give  a  coefficient  of  | AM/M |  in  the  range 
of  3-5  for  each  of  the  error  estimates  in  p,  cv^,  and  o^.  The  inherent 
accuracy  of  this  experiment  is  then  very  good  since  the  errors  in  the  measure¬ 
ments  | AM/M |  could  be  made  small  depending  upon  the  magnitude  of  the  mea¬ 
surable  quantities  and  the  measurement  techniques  employed.  In  addition,  if 
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measurements  at  one  altitude  could  be  correlated  with  those  made  at  other 
altitudes,  the  statistical  error  values  could  be  made  even  lower.  The  cor¬ 
relation  of  measurements  made  at  different  altitudes  would  require  the  use 
of  an  assumed  atmospheric  density  model.  The  magnitude  of  measurable  quan- 
ities  in  a  satellite  experiment  is  discussed  in  the  next  section. 


6.3.  Magnitudes  of  Measurable  Quantities  in  a  Satellite  Experiment 

In  a  satellite  experiment  such  as  described  above,  the  quantities 
*  *  * 


D  ,  T  ,  and  T  would  bo  determined  from  measurements  made  on  the  orbital 
P  s 

decay  rate,  the  satellite  precession  rate  and  the  satellite  spin  decay  rate. 
An  estimate  of  the  magnitudes  of  these  rates  can  be  determined  from  the  ex¬ 
ample  results  obtained  for  the  cone-disk  satellite. 

* 


Consider  first  the  determination  of  D  which  contains  the  drag  D, 

For  the  case  of  circular  orbits,  assuming  that  the  altitude  losses  due  to 

drag  effects  are  small,  the  decay  rate  of  the  orbital  radius,  r  ,  can  be 

o 

approximated  by  the  following  expression  given  in  Reference  1, 


Arn  9 

—  =  -  4ttB  p  (r  ) 
rev  *av  o  av 


6.16 


where  B  is  the  ballistic  coefficient  CpA/2m,  m  is  the  mass  of  the  satellite, 
Pav  is  the  average  atmospheric  gas  density  at  the  average  orbital  radius  r  . 
The  mass  m  of  the  satellite  is  directly  proportional  to  the  density  p  of  the 
material  used  tc  construct  the  satellite.  Let  the  proportionality  be  ex¬ 
pressed  as 

*  — 

m  =  p„  r,  A  6.17 


where  A  is  the  reference  area  of  the  satellite  and  r,  has  units  of  length  and 
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is  the  proportionali  factor  determined  by  the  design  of  the  satellite. 
Using  Equation  6.17  a1  he  definition  of  B,  Equation  6.16  becomes,  dropping 
the  "av"  subscripts  on  p  and  r  , 


Ar  *  r  _  pr 

—/rev  -  -D  =  -2tt  CD  — | 

o  psr,  psr, 
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The  convenient  non-dimensional  form  of  Equation  6.18  can  also  be 
obtained  for  the  precession  rate  and  slow  down  rate  of  the  satellite.  For 
high  spin  rates,  the  Luler  equations  for  the  dynamics  of  a  spinning  body 
under  the  action  of  external  torques  reduce  to 

6.19 


T.  /  =  I.Q  to .  / 
i  1  x 
s  s 


T.  /  =  1,0  u>.  / 
J  1  i 
s  s 


V  =  I2(fi+ 


6.20 


6.21 


where  1^  is  the  moment  of  inertia  about  an  axis  perpendicular  to  the  spin 

axis,  I 2  is  the  moment  of  inercia  about  the  spin  axis,  Q  is  the  spin  rate 

and  the  terms  uj. /,  c and  uj  /  represent  the  components  of  precession  rate 

J  L  K 

S  S  S 

of  the  spinning  body.  Assume  for  the  purposes  of  the  present  study  that  the 

satellite  car  be  designed  to  be  iso-inertial,  1^  =  I^,  regardless  of  the 

external  shape.  In  terms  of  the  notation  of  the  preceding  chapters,  the 

pre:ession  rate  /  is  equal  to  X.  Therefore,  Equation  6.20  becomes 
s 


AX/rev  =  2tt  T.  //QI.ui 

i  i  o 
s 


6.22 


where  u>o  is  the  angular  velocity  of  the  orbit.  The  terms  in  Equation  6.22 
can  be  expressed  as  follows 
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Ti'=CT.  2  ?£kt 
s  Js 


Q  =  U  K/r 

CO 


U  =  r  cu 

co  oo 

substituting  for  T./,  uj  ,  and  0  into  Equation  6.22,  the  following  is  obtained 
Js  ° 

C„ 

j  T  2 

AX/ rev  =  2n  — ^  p  r  6.23 

K  O 

Similar  to  what  was  done  for  the  drag  case,  assume  that  the  moment  of  inertia 

is  proportional  to  the  material  density  of  the  satellite,  p  .  Let 

s 

I.  =  |  p  A  72  r  *  6.24 

1  l  S  £■ 

■k 

where  r^  has  units  of  length  and  is  the  proportionality  factor.  The  quan- 

*  * 
tity  x^  may  or  may  not  be  equal  to  the  similar  term  r^  which  was  used  for 

the  drag  equations.  Substituting  Equation  6.24  in  Equation  6.23,  the  follow¬ 
ing  is  obtained 


AX/rev  =  f 


P  Pcr: 


Tjs  pr0 
_  0  _ £  o 

"  "  K  P  r  o* 


In  a  similar  fashion,  the  satellite  spin  decay  rate  is  found  to  be 


AO.  *  rc,  .  K  P'o 

_/reu  T<  —  2*  _  — 


where  the  term  r^*  would  be  equal  to  x^*  for  an  iso-inertial  body. 

Summarizing,  the  rates  of  change  in  the  measurable  quantities  of  a 


-...a 
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satellite  experiment  are  found  expressible  in  the  following  non-dimensional 
forms . 


Ar 


7rev  =  -2ttC. 


p  ro 
'D  psrl* 


6.27 


Tj'  p  r0 

AX/rev  =  2tt  — 

rS  2 


6.28 


T.  / 

an,  .  ks  p  r» 

g'/rev  =  2tt  -j—  — ; 

S  j 


6.29 


The  products  p  r*,  p  r  *,  and  p  r  s*  are  quantities  which  may  be 

controlled  by  the  satellite  design.  In  general,  large  solid  satellites  would 

have  p  r*  values  which  are  large  whereas  light  satellites  would  have  small 
s 

values  of  pgr*.  In  terms  of  orders  of  magnitude,  an  upper  limit  for  a  large 

2 

solid  satellite  would  be  in  the  order  of  100  gm/cm  .  A  lower  limit  for  pgr* 
-2  2 

could  be  10  gm/cm  for  a  thin  walled  hollow  satellite. 

The  product  p  rQ  varies,  of  course,  with  the  orbital  altitude  of 
the  satellite.  Figure  6.2  shows  the  variation  of  p  rQ  for  a  high,  low,  and 
medium  density  at  ^sphere  versus  the  orbital  altitude  h,  where  h=ro-re;  r^  = 
(radius  of  Earth). 

Using  some  representative  values  of  C^,  ,  and  ,  which  were 

j  '  k ' 

obtained  in  the  numerical  studies  of  chapter  5  on  the'’cone-disk  satellite, 
the  plots  given  in  Figures  6.3,  6.4,  and  6.5  show  bow  the  measurable  quan¬ 
tities  vary  with  altitude  h,  for  pgr*  =  1.  Figure  6.3  can  also  be  inter¬ 
preted  in  terms  of  At/t,  the  rate  of  change  in  the  period  of  the  orbit.  The 
rate  of  change  in  period  is  related  to  the  altitude  decay  rate  by  the 


yOr0  (gm /cm2  ) 
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Figure  6.2.  High,  low,  and  medium  atmospheric  density  variation  with 
altitude  (from  tables  in  Reference  29). 


Altitude,  (km) 


Figure  6.5.  Relative  average  spin  decay  of  a  cone  disk  satellite  per 
(circular)  orbit. 
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following  equation, 

,  -  A  r 

^•/rev  =  —  ~ ^/rev  6.30 

o 

The  results  shown  in  Figure  6.4  are  given  in  terms  of  sec  of  arc/rev. 

The  results  given  in  Figures  6.3,  6.4,  and  6.5,  clearly  illustrate 

the  feasibility  of  proposed  experimental  techniques  in  terms  of  the  magnitude 

of  measurable  quantities.  In  fact,  the  measurements  could  all  be  made  using 

an  optical  technique  such  as  that  proposed  by  the  Coordinate!  Science  Labora- 
12, 13 

tory .  1  The  advantage  to  an  optical  readout  system  is  that  the  satellite 
itself  could  be  completely  passive  in  that  sunlight  reflected  off  the  surface 
of  the  satellite  can  be  used  to  determine  the  orientation  of  the  spin  axis 
and  thus  be  able  to  determine  the  precession  rate  of  the  satellite  spin  axis. 
The  technique  proposed  by  CSL  is  fully  described  in  the  References  12  and  13. 
The  results  of  this  extensive  study  made  by  CSL  of  the  feasibility  of  the 
optical  readout  technique  indicate  that  precession  rates  of  the  order  shown 
in  Figure  6.4  would  be  well  within  the  capabilities  of  the  optical  technique. 
Satellite  spin  rate  data  could  be  obtained  directly  from  the  observations  of 
reflected  sunlight  from  the  satellite.  The  optical  technique  is,  of  course, 
a  standard  method  of  obtaining  orbital  drag  data. 

6.4.  Assessment  of  Results  Obtained  From  a  Satellite  Experiment 

In  the  preceding  discussions,  a  satellite  experiment  to  determine 
the  parameters  of  the  Maxwell  model  was  studied  mainly  to  illustrate  its 
feasibility.  However,  the  accuracy  of  the  interpretation  of  the  measurements 
depends  upon  how  well  the  assumed  model  approximates  the  actual  reflection 
phenomena . 
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In  order  to  illustrate  how,  for  a  special  case,  erroneous  result; 
could  be  obtained,  take  an  experiment  in  which  the  measurements  are  inter¬ 
preted  in  terms  of  the  parameters  of  the  Maxwell  model.  As  an  extreme  case 
consider  that  the  actual  reflection  process  is  a  perfect  backscatter.  That 
is,  irrespective  of  the  angle,  0,  at  which  the  molecules  impinge  on  the  sur¬ 
face,  the  molecules  are  always  reflected  back  in  the  same  direction.  Perfect 
backscatter  is  a  type  of  reflection  which  cannot  be  approximated  by  the  Max¬ 
well  model  (or  the  Schamberg  model  either). 

The  procedure  outlined  for  an  experiment  using  a  cone-disk  satel¬ 
lite  would  yield  values  of  Qf^,  0^,  and  p  for  this  extreme  case  of  perfect 
backscatter;  however,  it  can  be  shown  that  or^  and  would  be  equal  to  one 
and  the  gas-density  would  be  twice  its  actual  value. 

The  above  example  serves  to  illustrate  the  importance  of  inter¬ 
preting  the  satellite  measurements  in  terms  of  a  model,  such  as  the  gen¬ 
eralized  model,  which  can  cover  a  wide  range  of  possible  gas-surface  inter¬ 
actions.  In  the  range  of  possible  reflections  between  specular  and  diffuse, 
the  results  of  chapter  4  indicate  that  use  of  either  the  generalized  model 
or  the  Maxwell  model  could  be  expected  to  yield  results  which  could  reason¬ 
ably  approximate  the  actual  reflection  process.  For  the  case  of  backscatter 
or  over  specular  reflections,  however,  use  of  either  the  Maxwell  or  the 
Schamberg  model  could  not  be  expected  to  yield  valid  results  on  the  gas 
surface  interaction  parameters.  The  Nocilla  model  also  could  give  ex roneous 
results  if  the  actual  distribution  function  were  far  from  the  drifting  Max¬ 
wellian  assumed  by  Nocilla. 

Solutions  for  the  unknowns  when  the  generalized  model  is  used 
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must  be  obtained  using  numerical  methods  since  in  general  the  resulting 
equations  are  trigometric  in  the  parameter  P^  .  Trigometric  equations  can, 
however,  lead  to  multiple  solutions  for  a  given  set  of  measurements.  If 
multiple  solutions  are  obtained,  then  additional  experiments  would  have  to 
be  performed  which  are  designed  to  detect  the  correct  solution.  It  must  be 
concluded  that,  even  considering  the  possibility  of  multiple  solutions,  the 
use  of  the  generalized  model  in  the  interpretation  of  satellite  data  is  to 
be  favored  over  the  use  of  a  given  particular  model.  More  ipecific  models 
could  be  employed  after  the  analysis  using  the  generalized  model  has  indi¬ 
cated  the  general  character  of  the  reflection  process. 

In  order  to  illustrate  the  application  of  the  generalized  model, 
consider  a  reflection  which  is  assumed  to  be  composed  of  three  separate 
beams*,  beam  1  (j  =  l)  is  in  the  backscatter  direction  (P^=2) ,  beam  2  (j=2) 
is  normal  to  the  surface  (?£=!) ,  and  beam  3  (j=3)  is  in  the  specular 
direction  (P3=0)  (see  Figure  6.8).  The  effect  on  the  drag  of  a  body  due 
to  a  reflection  of  this  type  can  be  determined  from  the  equations  given 
in  section  4.3.1.  For  the  values  of  P^.  given  above,  the  drag  coefficient 
for  a  flat  plate,  cylinder,  cone  and  sphere  are  given  by, 

(Vplate  =  sin  9s  (2  +  2V*"*1  +  2  Sin  9s  °2^2 

-  2  cos  20g  ct^/1-0'3  6.31 

(V cylinder  =  2  +  2° +  ?  °2^2  +  3  a3^3 

(Vcone  =  2  +  2ai'/l-«1  +  2  sin  6  l-a2 
-  2  cos  26 


6.33 
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'Vsphere  =  2  +  +  !  6'34 
These  results  serve  to  illustrate  the  statement  made  at  the  end  of  chapter  2 
that  the  aerodynamic  properties  are  specified  by  a  choice  of  only  two  quan¬ 


tities  for  each  reflected  beam;  P.  and  o.Jl-a,.  The  above  equations  show 

J  J  j 


that  the  quantities  <y.  and  ,/l-a..  always  appear  as  the  product  in  the  ex¬ 
pressions  for  aerodynamic  properties  and  could  therefore  be  considered  as  a 
single  parameter. 

The  results  given  in  Equations  6.31,  6.32,  6.33,  and  6.34  also 
illustrate  tha£,  without  prior  knowledge  of  either  or  a ^ ,  only  the  quan¬ 
tities  a  sjl-a ^  could  be  determined  in  an  experiment  which  utilizes  measured 
aerodynamic  properties  of  convex  shapes.  In  application  to  satellite  experi¬ 
ments  then,  if  no  assumptions  are  made  on  the  distribution  function  cf  re¬ 
flected  molecules  (such  assumptions  are  made  in  the  Maxwell,  Schamberg  and 
Nocilla  models),  the  results  of  a  satellite  experiment  using  convex  shapes 

are  limited  to  obtaining  information  on  the  number  of  reflected  beam  com- 

'V. 


ponents,  the  magnitude  of  each  component,  and  the  direction  of  each  compor 
nent.  This  information,  even  though  limited,  would  be  sufficient  to  deter¬ 
mine  the  aerodynamic  drag,  lift,  and  torque  properties  of  any  convex  shaped 
body,  spinning  or  non-spinning.  Such  information  would  also  be  utilized  to 
design  future  satellite  experiments  to  obtain  more  refined  information  on 
the  gas  surface  interaction.  Although  such  shapes  are  not  considered  in 

this  study,  it  is  proposed  that  concave  shapes  could  be  utilized  to  ob'tain 

% 

information  on  the  distribution  function  of  reflected  molecules  since  the 
aerodynamic  properties  of  these  shapes  are  dependent  upon  the  properties  of 
reflected  distribution. 
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6.5.  Significance  of  Satellite  Experiments  to  Determine  Gas-Surface  Inter¬ 
action  Parameters  and  the  Atmospheric  Density 

The  examples  given  in  this  chapter  have  employed  idealized  con¬ 
ditions  such  as  circular  orbits  and  constant  density  which  would  not  be  the 
case  in  an  actual  experiment.  As  mentioned  in  the  introduction  to  this 
chapter,  a  number  of  uncertainties  are  associated  with  satellite  experi¬ 
ments  which  can  affect  the  interpretation  of  results.  Two  of  the  major  un¬ 
certainties  have  already  been  adequately  discussed;  the  uncertainty  associ¬ 
ated  with  the  atmospheric  density  was  discussed  in  the  introduction  and  the 
uncertainties  associated  with  the  gas-surface  'r.teraction  model  has  been 
discussed  in  the  preceding  section.  The  significance  of  performing  a  satel¬ 
lite  experiment  in  light  cf  some  additional  uncertainties  will  be  discussed 
here.  To  be  considered  are  uncertainties  concerning  the  condition  of  the 
satellite  surface  (ie.  degree  of  surface  contamination  by  adsorbed  gases, 
composition  cf  surface-adsorbed  gases,  and  roughness  of  surface),  composi¬ 
tion  of  the  atmosphere,  and  variation  of  atmospheric  density  with  altitude. 

6.5.1.  Consideration  of  Surface  Conditions 

The  gas-surface  interaction  is  known  to  depend  upon  the  degree 
ana  composition  of  adsorbed  surfaces  (see,  for  example,  the  experimental 
results  given  in  Reference  30)  which,  in  the  satellite  environment,  is  not 
well  known.  The  satellite  environment  in  some  respects  acts  like  a  cleans¬ 
ing  environment  in  that  the  very  hrgh  vacuum  combined  with  the  effects  of 
high  energy  solar  and  cosmic  rays  tend  to  rid  the  surface  of  trapped  gas 
molecules.  On  the  other  hand,  the  constant  bombardment  of  high  velocity  gas 
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molecules  as  the  satellite  travels  in  its  orbit  causes  the  surface  to  again 
be  contaminated  with  trapped  gas  molecules.  The  cleansing  and  contaminating 
processes  take  place  continuously.  The  equilibrium  conditions  if  reached 
have  not  as  yet  been  determined.  This  uncertainty  enters  the  interpretation 
of  results  when  measurements  made  at  one  orbital  altitude  are  correlated 
to  those  made  at  other  altitudes. 

Instead  of  treating  the  surface  conditions  as  uncertainties  in  the 
experiment,  it  is  suggested  that  a  satellite  experiment  be  designed  to  ob¬ 
tain  information  on  the  uncertainties.  One  possible  approach  is  to  first 
contaminate  the  satellite  surface  with  a  known  contamination  which  would  be 
expected  to  degas  at  a  known  rate  in  the  relative  vacuum  of  the  orbital  en¬ 
vironment.  The  degassing  rate  could  then  be  correlated  with  the  measured 
rates  of  change  in  the  aerodynamic  properties  of  the  satellite  to  determine 
the  effect  of  contamination  on  the  gas-surface  interaction.  That  such  an 
approach  would  be  feasible  is  indicated  first  by  the  results  obtained  in 

this  chapter  on  the  possible  accuracy  of  a  satellite  experiment.  Secondly, 

31 

in  some  preliminary  studies  performed  by  Cohen  ,  it  was  concluded  that 
contaminates  such  as  water  vapor  on  a  metal  surface  could  be  expected  to 
degas  down  Lo  1/e  of  the  initial  surface  coverage  in  a  time  of  about  5  to 
10  weers.  This  rate  of  change  estimate  combined  with  the  expected  accuracy 
of  the  satellite  measurements  indicates  that  such  an  experiment  may  be 
feasible . 

The  effects  of  other  surface  properties  on  the  gas-surface  inter¬ 
action  could  also  be  studied  by  making  accurate  measurements  of  the  aero¬ 
dynamic  properties  as  a  function  of  time  and  then  correlating  the  results 
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with  known  variations  of  the  environmental  conditions.  For  example,  for 
satellites  with  long  lifetimes,  the  satellite  surface  could  be  expected  to 
become  rougher  with  time  due  to  the  bombardment  of  micrometeorites,  dust 
particles  and  high  energy  cosmic  rays  (see  for  example  Reference  32).. 

The  effect  of  the  gas- surface  interaction  on  surface  roughness  could  then  be 
assessed  by  correlating  changes  in  the  satellite  aerodynamic  properties  with 
„hat  is  known  about  the  roughening  effects  of  the  space  environment.  Data 
from  satellites  which  have  changing,  or  controlled  surface  temperatures 
could  be  utilized  in  a  similar  manner  to  obtain  information  on  the  effect 
of  surface  temperature. 

6.5.2.  Consideration  of  the  Composition  of  the  Atmosphere 

The  composition  of  the  earth's  atmosphere  is  known  to  vary  con¬ 
siderably  with  altitude  and  solar  activity  (see  for  example  Reference  29)  . 
The  average  molecular  weight  at  1000  km.  varies  from  1.47  for  a  low-dencity 
atmosphere  to  a  valu.  15.04  for  a  high  density  atmosphere.  At  300  km. 
the  variation  is  fron  16.89  to  22.46  for  the  low  and  high  density  atmo¬ 
spheres.  (Values  of  molecular  weight  obtained  from  Reference  29) .  As  a 
satellite  orbit  decays,  then,  the  species  of  gas  molecules  which  impinge 
on  the  satellite  surface  will  change  in  concentration.  The  change  in  gas- 
specieo  concentration  with  altitude  enters  into  the  interpretation  ,>f  re¬ 
sults  when  measurements  made  at  one  altitude  are  correlated  to  those  made 
at  other  altitudes.  Two  uncertainties  are  involved;  the  gas-surface  inter¬ 
action  as  a  function  of  gas-species  concentration  and  the  concentration  of 
gas-species  as  a  function  of  altitude. 
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As  with  the  uncertainties  in  surface  conditions,  it  is  proposed 
here  that  the  uncertainties  in  the  atmospheric  composition  be  considered  un¬ 
known  and  be  determined  by  satellite  experiments.  One  possible  approach  is 
to  consider  each  species  of  gas  separately.  That  is,  instead  of  the  average 
density  of  the  atmosphere,  p,  consider  that  the  free  stream  density  is  com¬ 
posed  of  a  separate  density,  p^,  for  each  species  where  p  =  Sp^.  In  addi¬ 
tion,  for  each  species  of  gas,  i,  consider  the  gas— surface  interaction  to  be 

determined  by  parameters  of  the  generalised  model  (cr -a/I-g;.  )  .  and  (P.)..  For 

J  j  J 

example,  the  drag  of  a  non-spinning  sphere  would  then  be  given  by 

TT 


Dsphere 
h  U2  A 

CO 


r  —  r4(1'cos  2  (pi)i)Ti 

=  l  pi  l2  +  j  (V1"Qfj).  L  (p^.c^-CPj).)  Jl 
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Expressions  such  as  Equation  6.31  could  also  be  developed  for  other  satellite 

shapes.  Then,  using  the  same  procedures  as  outlined  earlier  in  this  chapter, 

it  is  proposed  that  the  unknowns  p.,(P.).,  and  ( a .Jl-a .)  could  be  deter- 

1  J  J  ^ 

mined  from  measurements  made  on  Lhe  aerodynamic  properties  of  satellites. 
Certainly  for  such  a  large  number  of  unknowns,  a  single  satellite  experiment 
would  not  provide  sufficient  information  for  the  determinations.  However,  a 
number  of  properly  designed  satellites  could  be  utilized  along  with  a  care¬ 
ful  re-evaluation  of  past  drag  data  to  obtain  a  considerable  amount  of  infor¬ 
mation  on  the  unknowns.  The  values  of  p.  could  be  correlated  with  models 

i 

of  the  variation  in  the  composition  of  the  atmosphere  to  determine  which 
model,  if  @ny,  gives  the  most  consistent  results.  Such  an  analysis  would 
then  provide  information  on  both  the  atmospheric  composition  and  the  gas- 
surface  interaction  as  a  function  of  gas-species. 
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6.5.3.  Variation  of  Atmospheric  Density  Models 

In  the  determination  of  the  atmospheric  density,  a  model  of  the 
relative  variation  of  atmospheric  density  with  altitude  must  be  assumed. 
Such  a  model  must  be  utilized  in  the  analysis  of  elliptic  orbit's  and  also 
in  correlating  measurements  made  at  one  altitude  with  those  made  at  other 
altitudes.  An  error  in  the  determination  value  of  p  is  then  introduced  due 
to  the  uncertainty  associated  with  the  assumed  m6del  of  the  atmosphere. 

As  with  the  other  uncertainties  in  a  satellite  experiment,  it  is 
suggested  that  this  uncertainty  could  also  be  removed  by  a  number  of 
satellite  experiments  and  re-evaluation  of  past  drag  data.  The  procedure 
would  be  much  the  same  as  outlined  in  the  previous  section  on  atmospheric 
composition  models. 
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7 .  SUMMARY  AND  CONCLUSIONS 

A  generalized  model  of  the  gas-surface  interaction  was  developed 
to  cover  a  wide  range  of  possible  interactions.  This  model  was  incorporated 
into  the  aerodynamic  equations  of  spinning  convex  bodies  in  a  free  molecule 
flow  in  order  to  study  fully  the  influence  of  the  gas-surface  interaction  on 
the  aerodynamic  properties  of  satellites  in  this  flow  regime. 

Analysis  of  the  aerodynamic  properties  of  four  spinning  bodies, 
(disk,  cylinder,  cone,  and  sphere)  at  angles  of  attack  revealed  the  strong 
influence  of  the  gas-surface  interaction,  especially  on  the  torque  properties. 
The  aerodynamic  torque  acting  on  a  body  in  free  molecular  flow  was  found  to 
be  caused  by  1)  the  moment  of  drag  and  lift  forces  about  the  center  of  mass 
of  the  body  and  2)  forces  tangent  to  the  surface  induced  by  the  spinning  of 
the  body.  Aerodynamic  torques  of  the  first  type  are  experienced  by  both 
spinning  and  non-spinning  satellites  and  are  well  known.  Aerodynamic  torques 
of  the  second  type  which  are  spin  induced  have  components  both  parallel  and 
perpendicular  to  the  spin  axis  of  the  body.  The  component  of  orque  parallel 
to  the  spin  axis  would  cause  the  expected  decay  in  the  spin  'te  of  a  spin¬ 
ning  body.  The  components  of  torque  perpendicular  to  the  spin  axis  would 
cause  a  gyroscopic  precession.  The  spin  induced  torque  on  a  body  was  found 
to  be  more  strongly  dependent  on  the  gas-surface  interaction  than  aerodynamic 
torques  of  the  first  type. 

Spin  induced  effects  were  also  found  in  the  aerodynamic  drag  and 
lift  properties  of  spinning  bodies  at  angles  of  attack.  In  general,  it  was 
found  that  spinning  bodies  experience  higher  values  of  drag  and  lift  than  do 
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non- spinning  bodies.  Also  of  interest  is  a  lateral  force  experienced  by 
spinning  bodies  in  free  molecular  flow  which  is  entirely  spin  induced.  The 
spin  induced  lateral  force  was  found  to  be  also  strongly  dependent  on  the 
gas-surface  interaction,  as  for  the  spin  induced  aerodynamic  torques.  The 
aerodynamic,  spin  induced  lateral  force  is  analogous,  but  opposite  in  direc¬ 
tion,  to  the  Magnus  effect  on  spinning  bodies  in  viscous  flow. 

The  development  of  the  generalized  gas-surface  interaction  model 
and  che  analysis  of  aerodynamic  properties  of  spinning  bodies  formed  the 
basis  for  proposing  satellite  experiments  to  obtain  information  on  the  gas- 
surface  interaction  as  well  as  the  orbital  atmospheric  density.  It  was 
found  that  the  average  aerodynamic  properties  of  spinning  satellites  are 

V 

strongly  dependent  on  the  parameters  of  a  given  gas-surface  interaction  model. 
It  is,  therefore,  proposed  that  the  measured  average  aerodynamic  properties 
of  spinning  satellites  be  utilized  to  determine  precise  values  of  the  gas- 
surface  interaction  parameters  and  the  orbital  gas  density.  The  preliminary 
phase  of  the  study  of  the  feasibility  of  these  satellite  experiments  was 
conducted.  This  phase  of  the  study  covered  1)  the  consideration  of  schemes 
utilizing  the  aerodynamic  properties  of  satellites  of  various  shapes  and 
orientations,  2)  the  assessment  of  the  accuracy  of  determining  the  gas  sur¬ 
face-interaction  parameters  and  the  orbital  gas  density,  3)  the  estimate  of 
the  magnitude  of  measurable  quantitites  in  a  satellite  experiment,  and  4) 
the  investigation  of  the  possible  effects  on  accuracy  introduced  by  uncertain¬ 
ties  in  the  space  environment  and  satellite  surface  conditions.  On  the  basis 
of  these  results,  it  was  established  that  the  proposed  satellite  experiments 
are  feasible  and  could  provide  signal i-cant  information  on  both  the  gas-surface 
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interaction  at  satellite  velocities  and  the  near-earth  atmosphere. 

l'ue  to  the  considerable  uncertainties  associated  with  the  gas-sur- 

© 

face  interaction  at  satellite  velocities,  it  was  found  that  the  generalized 
modul  of  the  gas-surface  interaction  is  necessary  in  the  interpretation  of 
results  of  the  proposed  satellite  experiments.  Since  the  aerodynamic  pro¬ 
perties  of  convex  bodies  in  a  free  mol.  flow  are  not  dependent  on  the  spe¬ 
cific  form  of  the  distribution  of  reflected  molecules,  gas-surface  inter¬ 
action  models  which  incorporate  a  specific  distribution  function  are  not  only 
unnecessary  but  also  undesirable  in  that  considerable  error  can  be  introduced 
in  the  interpretation  of  aerodynamic  measurements  by  using  such  models.  For 
this  reason,  the  use  of  the  generalized  model  developed  in  this  study  is  pre¬ 
ferred  since  no  assumption  was  made  on  the  distribution  of  reflected  molecules 
other  than  the  existance  of  an  average  velocity  and  direction.  It  was  sug¬ 
gested  that  the  generalized  model  could  also  be  applied  to  the  interpretation 
of  results  obtained  from  laboratory  experiments  in  order  to  parameterize  in. 
a  general  manner  the  results  of  molecular  beam  studies.  The  parameterization 
of  these  results  would  facilitate  the  comparison  of  the  various  results  and 
could  serve  as  a  basis  for  suggesting  more  precise  gas-surface  interaction 
models.  Laboratory  experiments  on  gas-surface  interaction  which  measure 
forces  and  torques  are  particularly  suitable  in  using  the  generalized  model 
to  interpret  the  results. 

The  study  of  feasibility  of  the  satellite  experiments  proposed  in 
this  study  has  been  supported  by  the  National  Aeronautics  and  Space  Adminis¬ 
tration  at  the  George  Marshall  Space  Flight  Center.  The  aerodynamic  pro¬ 
perties  of  near-earth  satellites  are  of  major  importance  in  determining  both 
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the  lifetime  of  a  satellite  and  the  motion  of  the  satellite  about  its  center 
of  mass.  Since  NASA  is  planning  in  the  future  to  orbit  satellites  of  larger 
size  (such  as  manned  space  stations)  and  satellites  requiring  greater  orien¬ 
tation  stability  than  satellites  of  the  past,  there  is  an  urgent  need  for 
precise  knowledge  of  the  character  of  the  gas-surface  interaction  and  the 
composition  of  the  atmosphere  in  order  to  properly  and  economically  design 
these  future  satellites.  For  this  reason,  NASA  is  considering  an  extensive 

AO  A  7 

experimental  program  called  Project  ODYSSEY"3  which  is  planned  to  obtain 
information  on  the  gas-surface  interaction  and  the  orbital  environment  by  a 
number  of  satellite  experiments.  The  techniques  of  analysis  developed  and 
the  results  obtained  in  this  study  have  direct  application  to  the  design  of 
satellite  experiments  and  the  interpretation  of  results  of  these  experiments, 
as  well  as  the  interpretation  of  existing  satellite  data,  with  the  objective 
to  obtain  information  on  the  gas-surface  interaction  and  the  atmospheric 
composition . 
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In  rarefied  gas  flow  problems  there  are  two  effects  which  influence  the  flow; 
(1)  collisions  of  gas  molecules  with  other  gas  molecules,  and  (2)  collisions  of  gas 
molecules  with  solid  surfaces  (the  gas-surface  interaction).  This  study  deals  with 
free  molecular  flow  in  which  the  effect  of  gas-gas  collision  can  be  neglected  and 
the  gas- surface  interaction  has  the  dominate  influence  on  the  flow.  The  results 
°f  this  study  have  application  to  satellites  since  free  molecular  flow  conditions 
occur  at  orbital  altitudes  above  120km. 


Knowledge  of  the  gas-surface  interaction  is  required  in  order  to  determine  the 
aerodynamic  properties  of  satellites.  At  satellite  velocities  (7-8km/sec)  the 
interaction  of  neutral  atmospheric  gas  molecules  with  the  satellite  surface  occurs 
at  energies  in  the  1  to  10  ev  range.  It  is  just  this  energy  range  which  has  not 
been  sail* f actor ily  duplicated  in  the  laboratory;  therefore,  at  present,  laboratory 
gas-surface  interaction  data  can  not  be  applied  directly  to  the  determination  of  the 
aerodynamics  properties  of  satellites.  It  is  proposed  in  this  study  that  satellite 
experiments  be  performed  to  obtain  the  needed  information  from  measurements  of  the 
aerodynamics  properties  of  satellites.  In  order  to  interpret  the  satellite  data, 
a  generalized  gas-surface  interaction  model  was  developed  and  used  in  the  analysis 
of  this  study.  --  - 
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